April 17, 2008 15:40 WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



On Gamma-Ray Bursts 

R. Ruffini, 1 ' 2 ' 3 M.G. Bcrnardini, 1 ' 2 C.L. Bianco, 1 ' 2 ' 3 L. Caito, 1 ' 2 P. Chardonnet, 1 ' 3 ' 4 C. 

Cherubini, 1 ' 5 M.G. Dainotti, 1 ' 2 F. Fraschetti, 6 A. Geralico, 1 ' 2 R. Guida, 1 ' 2 B. Patricclli, 1 ' 2 M. 

Rotondo, 1 ' 2 J. A. Rueda Hernandez, 1 ' 2 G. Vereshchagin, 1,2 ' 3 S.-S. Xue 1 ' 3 

1 ICRANet and ICRA, Piazzale della Repubblica 10, 1-65122 Pescara, Raly. 
2 Dip. di Fisica, Universita di Roma "La Sapienza" , Piazzale Aldo Moro 5, 1-00185 Roma, Raly. 
3 ICRANet, Universite de Nice Sophia Antipolis, Grand Chateau, BP 2135, 28, avenue de 

Valrose, 06103 NICE CEDEX 2, France. 
4 Universite de Savoie, LAPTH - LAPP, BP 110, F-74941 Annecy-le-Vieux Cedex, France. 
5 Universita Campus Biomedico, Facolta di Ingegneria, Via Alvaro del Portillo 21, 00128 Rome, 

Italy. 

6 CEA/Saclay, F-91191 Gif-sur-Yvette Cedex, Saclay, France. 
E-mails: rujfini@icra.it, maria.bernardini@icra.it, bianco@icra.it, letizia.caito@icra.it, 
chardon@lapp.in2p3.fr, cherubini@icra.it, dainotti@icra.it, fraschetti@icra.it, geralico@icra.it, 
roberto.guida@icra.it, barbara.patricelli@icranet.org, michael.rotondo@icra.it, 
jorge.rueda@icra.it, veresh@icra.it, xue@icra.it. 

We show by example how the uncoding of Gamma-Ray Bursts (GRBs) offers unprece- 
dented possibilities to foster new knowledge in fundamental physics and in astrophysics. 
After recalling some of the classic work on vacuum polarization in uniform electric fields 
by Klein, Sauter, Heisenberg, Euler and Schwinger, we summarize some of the efforts 
to observe these effects in heavy ions and high energy ion collisions. We then turn to 
the theory of vacuum polarization around a Kerr-Newman black hole, leading to the 
extraction of the blackholic energy, to the concept of dyadosphere and dyadotorus, and 
to the creation of an electron-positron-photon plasma. We then present a new theoret- 
ical approach encompassing the physics of neutron stars and heavy nuclei. It is shown 
that configurations of nuclear matter in bulk with global charge neutrality can exist on 
macroscopic scales and with electric fields close to the critical value near their surfaces. 
These configurations may represent an initial condition for the process of gravitational 
collapse, leading to the creation of an electron-positron-photon plasma: the basic self- 
accelerating system explaining both the energetics and the high energy Lorentz factor 
observed in GRBs. We then turn to recall the two basic interpretational paradigms of 
our GRB model: I) the Relative Space-Time Transformation (RSTT) paradigm and 
2) the Interpretation of the Burst Structure (IBS) paradigm. These paradigms lead to 
a "canonical" GRB light curve formed from two different components: a Proper-GRB 
(P-GRB) and an extended afterglow comprising a raising part, a peak, and a decaying 
tail. When the P-GRB is energetically predominant we have a "genuine" short GRB, 
while when the afterglow is energetically predominant we have a so-called long GRB or 
a "fake" short GRB. We compare and contrast the description of the rclativistic expan- 
sion of the electron-positron plasma within our approach and within the other ones in 
the current literature. We then turn to the special role of the baryon loading in discrimi- 
nating between "genuine" short and long or "fake" short GRBs and to the special role of 
GRB 991216 to illustrate for the first time the "canonical" GRB bolometric light curve. 
We then propose a spectral analysis of GRBs, and proceed to some applications: GRB 
031203, the first spectral analysis, GRB 050315, the first complete light curve fitting, 
GRB 060218, the first evidence for a critical value of the baryon loading, GRB 970228, 
the appearance of "fake" short GRBs. We finally turn to the GRB-Supernova Time Se- 
quence (GSTS) paradigm: the concept of induced gravitational collapse. We illustrate 
this paradigm by the systems GRB 980425 / SN 1998bw, GRB 030329 / SN 2003dh, 
GRB 031203 / SN 20031w, GRB 060218 / SN 2006aj, and we present the enigma of the 
URCA sources. We then present some general conclusions. 
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1. Introduction 

After almost a century of possible observational evidences of general relativistic 
effects, all very weak and almost marginal to the field of physics, the direct obser- 
vation of gravitational collapse and of black hole formation promises to bring the 
field of general relativity into the mainstream of fundamental physics, testing a vast 
arena of unexplored regimes and leading to the explanation of a large number of 
yet unsolved astrophysical problems. 

There are two alternative procedures for observing the process of gravitational 
collapse: either by gravitational waves or by joint observation of electromagnetic 
radiation and high energy particles. The gravitational wave observations may lead 
to the understanding of the global properties of the gravitational collapse process, 
inferred from the time-varying component of the global quadrupolc moment of the 
system. Their observation is also made difficult and at times marginal due to the 
weak coupling of gravitational waves with the detectors. The observation of the 
electromagnetic radiation in the X, 7, optical and radio bands, and of the associated 
high energy particle component, is carried out by an unprecedented observational 
effort in space, ground and underground observatories. This effort is offering the 
possibility of giving for the first time a detailed description of the gravitational 
collapse process all the way to the formation of the black hole horizon. 

The Grossmann Meetings have been dedicated to foster the mathematical and 
physical developments of Einstein theories. They have grown in recent years due to 
remarkable progress both in fundamental physics and in astrophysics. In this sense, I 
will present here some highlights of recent progress on the theoretical understanding 
of Gamma-Ray Bursts (GRBs, for a review see e.g. Ref. 1 and references therein), 
which nicely represents two complementary aspects of the problem: progress in 
probing fundamental theories in yet unexplored regimes, as well as understanding 
the astrophysical scenario underlying novel astrophysical phenomena. 

ft is particularly inspiring that this MG11 takes place in Dahlcm, close to where 
many of the fundamental breakthroughs and discoveries of modern physics have 
indeed occurred. A few hundred meters from here, in Ehrcnbcrgstrasse 33, Albert 
Einstein once lived while introducing and developing the theory of general relativity 
(see Fig. f). 

A few hundred meters from here there is also the Kaiser- Wilhelm-Institut, where 
Lise Meitner, Otto Hahn and Fritz Strassmann (see Fig. 2) continued the work on 
neutron capture by Uranium initiated a few years before by Fermi. 2 This work led 
to the revolutionary evidence for Uranium fission. 3 ' 4 There is no need to stress the 
enormous consequences of this discovery following the work of Fermi, 5 Feynman, 
Metropolis & Teller, 6 Oppcnheimer, 7 Wigncr, 8 etc., and, just to keep symmetry, 
the work of Kurchatov, Sakharov, and Zel'dovich (see e.g. Ref. 9). The message 
from the nuclear and thermonuclear work leads to nuclear reactors and to explosive 
events typically of ~ 10 22 ergs/pulse (see Fig. 2). 

What 1 would like to stress in this lecture is the possible role of a theoretical 
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Fig. 1. Einstein's home in Dahlem 



work, coeval to the work of Otto Halm and Lisa Mcitncr, developed close to Dahlem, 
at the University of Leipzig. Such a theoretical work may very well lead in the 
future to the understanding of yet bigger explosions and have an equally, if not 
more, fundamental role on the existence and the dynamics of our Universe. We 
refer here to the work pioneered by Klein, 10 Sauter, 11 Euler and Heisenberg, 12,13 
Weisskopf 14,15 (see Fig. 3). The work deals with the creation of electron-positron 
pairs out of the vacuum generated by an overcritical electric field. In the following we 
give evidences that this process is indeed essential to the extraction of energy from 
the black hole and occurs in the GRBs. The characteristic energy of these sources 
is typically on the order of 10 49 -10 54 ergs/pulse. I am giving here, as examples of 
these GRBs, the light curves of GRB 980425 (with a total energy of ~ 10 48 ergs) 
and GRB 050315 (with a total energy of - 10 53 ergs). 

So much for fundamental science. From the astrophysical point of view, GRBs 
offers an equally rich scenario, being linked to processes like supernovae, coalescence 
of binary neutron stars, black holes and binaries in globular clusters, and possibly 
intermediate mass black holes. After reviewing some of the fundamental work on 
vacuum polarization, I will focus on a new class of 4 particularly weak GRBs which 
promises to clarify the special connection between GRBs and supernovae. I will then 
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Fig. 2. Left: Lisc Meitner and Otto Hahn in the Kaiser- Wilhclm-Institut in Berlin. Right: 
Picture of a nuclear reactor and an atomic bomb, both resulting from Uranium fission. 



conclude on the very surprising aspect that GRBs can indeed originate from a very 
wide variety of different astrophysical sources, but their features can be explained 
within a unified theoretical model, which applies in the above-mentioned enormous 
range of energies. The reason of this "uniqueness" of the GRBs is strictly linked to 
the late phases of gravitational collapse leading to the formation of the black hole 
and to its theoretically expected "uniqueness" : the black hole is uniquely character- 
ized by mass, charge and angular momentum. 16 The GRB phenomenon originates 
in the late phases of the process of gravitational collapse, when the formation of 
the horizon of the black hole is being reached. The phenomenon is therefore quite 
independent of the different astrophysical settings giving origin to the black hole 
formation. 

2. Vacuum polarization in a uniform electric field 

We recall some early work on pair creation following the introduction by Dirac of 
the relativistic field equation for the electron and leading to the classical results of 
Klein, 10 Hciscnberg & Euler, 13 Schwinger. 17-19 

2.1. Klein and Sauter work 

It is well known that every relativistic wave equation of a free relativistic particle 
of mass m e , momentum p and energy £ , admits symmetrically "positive energy" 
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Fig. 3. Above Left: A historical picture taken at the time of the Heisenberg-Euler work. Sur- 
rounding Hciscnbcrg, counterclockwise, are some of the major interpreters of modern physics: 
Rudolf Pcicrls, George Placzek, Giovanni Gentile, Giancarlo Wick, Felix Bloch, Wicky Wcisskopf 
and Fritz Sautcr. Above Center: H. Euler. Above Right: Julian Schwingcr. Below Left: Ob- 
served light curve of GRB 980425 with our theoretical fit with a total energy of 10 48 ergs. Below 
Right: Observed light curve of GRB 050315 with our theoretical fit with a total energy of 10 53 
ergs. 

and "negative energy" solutions. Namely the wave-function 

V' ± (x,t) ~e*( k - x - £±t ) (1) 

describes a relativistic particle, whose energy, mass and momentum must satisfy, 

4 =m 2 e c A + c 2 |p| 2 ; £± = ±^m 2 e c i + c 2 |p| 2 , (2) 

this gives rise to the familiar positive and negative energy spectrum (£±) of positive 
and negative energy states (^(x, t)) of the relativistic particle, as represented in 
Fig. 4. In such free particle situation (flat space, no external field), all the quantum 
states are stable; that is, there is no possibility of "positive" ("negative") energy 
states decaying into a "negative" ("positive") energy states, since all negative energy 
states are fully filled and there is an energy gap 2m e c 2 separating the negative energy 
spectrum from the positive energy spectrum. This is the view of Dirac theory on 
the spectrum of a relativistic particle. 20 ' 21 

Klein studied a relativistic particle moving in an external constant potential V 
and in this case Eq. (2) is modified as 

[£ - V] 2 = m 2 c A + c 2 |p| 2 , ; £ ± = V ± \/m 2 c 4 + c 2 |p| 2 . (3) 
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He solved this relativistic wave equation by considering an incident free relativistic 
wave of positive energy states scattered by the constant potential V, leading to 
reflected and transmitted waves. He found a paradox that in the case V > £ + m e c 2 , 
the reflected flux is larger than the incident flux j Te { > jj nc , although the total flux 
is conserved, i.e., ji nc = jref + jtran- This was known as the Klein paradox (see 
Ref. 10). This implies that negative energy states have contributions to both the 
transmitted flux j t ran and reflected flux j ro f ■ 



4 £ 



positive continuum £ + > m e c 2 

//////////////////// 



m e c 



-m e c 



77777777777777777777 



negative continuum £_ < m e c 2 



Fig. 4. The mass-gap 2m e c 2 that separates the positive continuum spectrum £+ from the negative 
continuum spectrum £_ . 



Sautcr studied this problem by considering an electric potential of an external 
constant electric field E in the z direction. 11 In this case the energy £ is shifted 
by the amount V(z) = —eEz, where e is the electron charge. In the case of the 
electric field E uniform between z\ and z-i and null outside, Fig. 5 represents the 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



7 



corresponding sketch of allowed states. The key point now, which is the essence of 
the Klein paradox, 10 is that the above mentioned stability of the "positive energy" 
states is lost for sufficiently strong electric fields. The same is true for "negative 
energy" states. Some "positive energy" and "negative energy" states have the same 
energy- levels, i.e., the crossing of energy-levels occurs. Thus, these "negative energy" 
waves incident from the left will be both reflected back by the electric field and partly 
transmitted to the right as a '"positive energy" wave, as shown in Fig. 5. 22 This 
transmission is nothing else but a quantum tunneling of the wave function through 
the electric potential barrier, where classical states are forbidden. This is the same 
as the so-called the Gamow tunneling of the wave function through nuclear potential 
barrier (Gamow- wall). 23 




Fig. 5. In the presence of a strong enough electric field the boundaries of the classically allowed 
states ("positive" or "negative") can also be so tilted that a "negative" is at the same level as 
a "positive" (level crossing). Therefore a "negative" wave-packet from the left will be partially 
transmitted, after an exponential damping due to the tunneling through the classically forbidden 
states, as a "positive" wave-packet outgoing to the right. This figure is reproduced from Fig. II in 
Ref. 22, and = m e c 2 , eV = V(z), uj = £. 



Sauter first solved the relativistic Dirac equation 20 ' 21 in the presence of the 
constant electric field by the ansatz, 



^(x,i) = e^+ fe ^- £±t )x S 3M 



(4) 



Where the spinor function Xs 3 { z ) obeys the following equation (70,71 are Dirac 
matrices) 



ftc7 3 ^ + jo(V(z) - £±) + (m e c 2 + ic^ 2 Py + ic~f\p x ) 



Xs 3 (z)=0, (5) 
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and the solution Xs 3 {z) can be expressed in terms of hypergeometric functions. 11 
Using this wave-function tp s (x,t) (4) and the flux ictplj^ips, Sauter computed the 
transmitted flux of positive energy states, the incident and reflected fluxes of nega- 
tive energy states, as well as exponential decaying flux of classically forbidden states, 
as indicated in Fig. 5. Using continuous conditions of wave functions and fluxes at 
boundaries of the potential, Sauter found that the transmission coefficient \T\ 2 of 
the wave through the electric potential barrier from the negative energy state to 
positive energy states: 

lr . 9 (transmission fluxl ^ m l° 3 ,„s 
1 1 |incidcnt flux| v ; 

This is the probability of negative energy states decaying to positive energy states, 
caused by an external electric field. The method that Sauter adopted to calculate 
the transmission coefficient |T| 2 is the same as the one Gamow used at that time to 
calculate quantum tunneling of the wave function through nuclear potential barrier 
(Gamow- wall), leading to the a-particle emission. 23 

2.2. Heisenberg-Euler-Weisskopf effective theory 

To be able to explain clastic light-light scattering, 12 Heisenberg and Euler 13 and 
Wcisskopf 14 ' 15 proposed a theory that attributes to the vacuum certain nonlinear 
electromagnetic properties, as if it were a dielectric and permeable medium. 13 ' 14 

Let C to be the Lagrangian density of electromagnetic fields E,B; a Legendre 
transformation produces the Hamiltonian density: 

U = E i § i -C. (7) 

In Maxwell's theory, the two densities are given by 

£ M = ^(E 2 -B 2 ), % = i(E 2 + B 2 ). (8) 
07T an 

To quantitatively describe nonlinear electromagnetic properties of the vacuum based 
on the Dirac theory, the above authors introduced the concept of an effective La- 
grangian £ ff of the vacuum state in the presence of electromagnetic fields, and an 
associated Hamiltonian density 

Aff = £m + AC, H cS = H M + AH. (9) 
From these one derives induced fields D, H as the derivatives 

In Maxwell's theory, AC = in the vacuum, so that D = E and H = B. In Dirac's 
theory, however, AC is a complex function of E and B. Correspondingly, the vacuum 
behaves as a dielectric and permeable medium 13 ' 1 in which, 

A = ^ tjkEk, Hi — ^^^ ik B k , (11) 
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where complex e,fc and fiik are the field-dependent dielectric and permeability ten- 
sors of the vacuum. 

The discussions on complex dielectric and permeability tensors (e^ and /ijfc) can 
be found for example in Rcf. 24. The effective Lagrangian and Hamiltonian densities 
in such a medium is given by 

£ eff = ^(E-D-B-H), W off = ^-(E-D + B.H). (12) 
In this medium, the conservation of electromagnetic energy has the form 

- divS =±( E .f + B.§). S-£.x„. (1.) 

where S is the Poynting vector describing the density of electromagnetic energy 
flux. By considering electromagnetic fields complex and monochromatic 

E = E(w) cxp —i{ujt); B = B(w) exp —i(ujt), (14) 

of frequency u>, the dielectric and permeability tensors are frequency-dependent, i.e., 
Cik(u)) and fiik(co). Substituting these fields and tensors into the r.h.s. of Eq. (13), 
one obtains the dissipation of electromagnetic energy per time into the medium, 

Q = {Im [ Cifc (w)] EiE* k + Im [ Mife (w)] B^} . (15) 

This is nonzero if ejfc(w) and iiikiui) contain an imaginary part. The dissipation 
of electromagnetic energy is accompanied by heat production. In light of the third 
thermodynamical law of entropy increase, the energy Q of electromagnetic fields 
lost in the medium is always positive, i.e., Q > 0. As a consequence, Im[ejk(u;)] > 
and Im[/ijfc(o;)] > 0. The real parts of ejfc(w) and ^^(u)) represent an electric and 
magnetic polarizability of the vacuum and lead, for example, to the refraction of 
light in an electromagnetic field, or to the elastic scattering of light from light. 12 The 
7iij(u>) = y/eik(uj)iikj{^) is the reflection index of the medium. The field-dependence 
of at and fiik implies nonlinear electromagnetic properties of the vacuum as a 
dielectric and permeable medium. 

The effective Lagrangian density (9) is a relativistically invariant function of the 
field strengths E and B. Since (E 2 — B 2 ) and (E • B) 2 are relativistic invariants, 
one can formally expand AC in powers of weak field strengths: 

AC = «: 20 (E 2 -B 2 ) 2 + «:o2(E-B) 2 + K 3o(E 2 -B 2 ) 3 + Kl2 (E 2 -B 2 )(E.B) 2 + . . . , (16) 

where fvjj are field-independent constants whose subscripts indicate the powers of 
(E 2 — B 2 ) and E • B, respectively. Note that the invariant E • B appears only in 
even powers since it is odd under parity and electromagnetism is parity invariant. 
The Lagrangian density (16) corresponds, via relation (7), to 

AH = k 2 ,o(E 2 - B 2 )(3E 2 + B 2 ) + Ko , 2 (E • B) 2 

+K 3 , (E 2 - B 2 ) 2 (5E 2 + B 2 ) + ki, 2 (3E 2 - B 2 )(E • B) 2 + . . . . (17) 
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To obtain H c s in Dirac's theory, one has to calculate 

An = E [ a ■ h^v + eA )) + P m - c2 \ ^} ' ( 18 ) 
k 

where an, (3 are Dirac matrices, A is the vector potential, and {tpk{x)} are the wave 
functions of the occupied negative-energy states. When performing the sum, one 
encounters infinities which were removed by Weisskopf, 14 ' 15 Dirac, 25 Heisenberg 26 
by a suitable subtraction. 

Heisenberg 26 expressed the Hamiltonian density in terms of the density ma- 
trix p(x,x') = ^ k ^P k {x)ipk{x') 2b Euler and Kockel, 12 and Heisenberg and Euler 13 
calculated the coefficients Kij. They did so by solving the Dirac equation in the 
presence of parallel electric and magnetic fields E and B in a specific direction, 

M%) - VW,«, S 3 = ei {xp '- £t) My)Xs a (x), n = 0,l,2,... (19) 

where {u n (y)} are the Landau states 27,28 depending on the magnetic field and 
Xs 3 {x) are the spinor functions calculated by Rcf. 11. Heisenberg and Euler used the 
Euler-Maclaurin formula to perform the sum over n, and obtained for the additional 
Lagrangian in (9) the integral representation, 



e 2 

AC 1 



-s ds 

16tt 2 he J s 3 



is 2 EB 



cos(s[E 2 - B 2 + 2i{EB)] 1 / 2 ) + c.c. 
cos(s[E 2 - B 2 + 2i(EB)} 1 / 2 ) - c.c. 



+ (^) , + yafli'-W)]. (») 

where E, B are the dimensionless reduced fields in the unit of the critical field E c , 

*-f- 

Expanding this in powers of a up to a 3 yields the following values for the four 
constants: 

a 2 rr 32na . 13 . . 

«2,0 = -^2 E c > «o : 2 = "k 2 ,o, «3 : o = -gjg- E c , «i )2 = y K 3,o- (22) 

Weisskopf 14 adopted a simpler method. He considered first the special case in which 
E = 0, B 7^ and used the Landau states to find AH of Eq. (17), extracting 
from this k 2 ,o and k 3 ^- Then he added a weak electric field E ^ to calculate 
perturbatively its contributions to AH in the Born approximation (see for example 
Landau and Lifshitz 27 ' 28 ). This led again to the coefficients (22). 

The above results receive higher corrections in QED and are correct only up to 
order a 2 . Up to this order, the field-dependent dielectric and permeability tensors 
eik and (11) have the following real parts for weak fields 

Re(e 4fe ) = Sit + ^ [2(E 2 - B 2 )S lk + 7BiB k ] + 0(a 2 ), 
45 

4o; - — — — 

Re(Mifc) = S ik + — [2(E 2 - B 2 )5 lk + 7E t B k ] + 0(a 2 ). (23) 
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2.3. Imaginary part of the effective Lagrangian 

Heisenberg and Euler 13 were the first to realize that for E ^ the powers series 
expansion (16) is not convergent, due to singularities of the integrand in (20) at 
s = tt/E, 2-k/E, .... They concluded that the powers series expansion (16) does not 
yield all corrections to the Maxwell Lagrangian, calling for a more careful evaluation 
of the integral representation, see Eq. (20). Selecting an integration path that avoids 
these singularities, they found an imaginary term. Motivated by Sauter's work 11 on 
Klein paradox, 10 Heisenberg and Euler estimated the size of the imaginary term in 
the effective Lagrangian as 

_^V(^)V^, (24) 

7T V h / 

and pointed out that it is associated with pair production by the electric field. 
This imaginary term in the effective Lagrangian is related to the imaginary parts 
of field-dependent dielectric e and permeability [i of the vacuum. 

In 1950's, Schwinger 17-19 derived the same formula (20) once more within the 
quantum field theory of Quantum Electromagnetics (QED), 



L aE 2 



V ir 2 h 

n=l 



and its Lorentz- invariant expression in terms of electromagnetic fields E and B, 

f _ ol£ 2 1 mr/3/e ( nnEA 

V - ^ 2-j n 2 tanhn7r/3/ £ 6XP ^ e ) ' W 

n=l ' y ' 



where 



1 1 J = ^V^ E2 - B2 ) 2 + 4(E • w± (E 2 - in 



(27) 



The exponential factor e 71 " »=e in Eqs. (6) and (24) characterizes the transmission 
coefficient of quantum tunneling, Rcf. 13 introduced the critical field strength E c — 

2 3 

^f£- (21). They compared it with the field strength E e of an electron at its classical 
radius, E e = e/r 2 where r e = ah/(m e c) and a = 1/137. They found the field 
strength E e is 137 time larger than the critical field strength E c , i.e., E e = a~ 1 E c . 
At a critical radius r c = a x / 2 h/ {m e c) > r e , the field strength of the electron would 
be equal to the critical field strength E c . 

As shown in Fig. 4, the negative-energy spectrum of solutions of the Dirac equa- 
tion has energies £- < —m e c 2 , and is separated from the positive energy-spectrum 
£ + > m e c 2 by a gap 2m e c 2 w 1.02MeV. The negative-energy states are all filled. 
The energy gap is by a factor 4/a 2 w 10 5 larger than the typical binding energy of 
atoms (~ 13.6eV). In order to create an electron-positron pair, one must spend this 
large amount of energy. The source of this energy can be an external field. 

If an electric field attempts to tear an electron out of the filled state the gap 
energy must be gained over the distance of two electron radii. The virtual particles 
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give an electron a radius of the order of the Compton wavelength A = h/m e c. Thus 
we expect a significant creation of electron-positron pairs if the work done by the 
electric field E over twice the Compton wave length h/m e c is larger than 2m e c 2 

eE (-—] > 2m e c 2 . 
\m e cj 

This condition defines a critical electric field 

2 3 

E c = ~ 1.3- 10 16 V/cm, (28) 

above which pair creation becomes abundant. To have an idea how large this critical 
electric field is, we compare it with the value of the electric field required to ionize 
a hydrogen atom. There the above inequality holds for twice of the Bohr radius and 
the Rydberg energy 

eE ion ( > a 2 meC 2 i 

\am e c ) 

so that E c w ££ on /a 3 is about 10 6 times as large, a value that has so far not been 
reached in a laboratory on Earth. 



3. Pair production in Coulomb potential of nuclei and heavy-ion 
collisions 

By far the major attention to build a critical electric field has occurred in the 
physics of heavy nuclei and in heavy ion collisions. We recall in the following some 
of the basic ideas, calculations, as well as experimental attempts to obtain the pair 
creation process in nuclear physics. 



3.1. The Z = 137 catastrophe 

Soon after the Dirac equation for a relativistic electron was discovered, 29-31 Gor- 
don 32 (for all Z < 137) and Darwin 33 (for Z = 1) found its solution in the point-like 
Coulomb potential V(r) — —Za/r, < r < oo. Solving the differential equations 
for the Dirac wave function, they obtained the well-known Sommerfeld's formula 34 
for the energy-spectrum, 

-1/2 



£(n,j) 



1 



Za 



n-\K\ + {K 2 - Z 2 a 2 y/ 2 



(29) 



Here the principle quantum number n = 1, 2, 3, • • • and 

-(/ + l),if j 
if 



K 



i i 



I > 
/ > 1 



(30) 



-(.7 + 1/2) = 

(i + 1/2) = Z, n 2 , 

where I — 0, 1, 2, ... is the orbital angular momentum corresponding to the upper 
component of Dirac bi-spinor, j is the total angular momentum, and the states with 
K = =Fl,=p2, =f3, ■ • -,T( n — 1) are doubly degenerate, while the state K = —n is 
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a singlet. 32,33 The integer values n and K label bound states whose energies are 
£(n,j) £ (0,m e c 2 ). For the example, in the case of the lowest energy states, one 
has 



S(lSi) = y/1 - (Za)2, (31) 



S(2S h )=S(2P h ) = J 1 + ^-^ 2 , (32) 



f(2Pj) = yi-i(Za)a. (33) 

For all states of the discrete spectrum, the binding energy mc 2 — £(n,j) increases 
as the nuclear charge Z increases, as shown in Fig. 6. When Z — 137, £(lS , 1 / 2 ) = 0, 
£{2S 1/2 ) = £(2P 1/2 ) = 1/V2 and £(2S 3/2 ) = V3/2. Gordon noticed in his pioneer 
paper 32 that no regular solutions with n = l, j = 1/2,1 = 0, and K = —1 (the 
I1S1/2 ground state) are found beyond Z = 137. This phenomenon is the so-called 
" Z =137 catastrophe" and it is associated with the assumption that the nucleus is 
point-like in calculating the electronic energy-spectrum. 




Fig. 6. Atomic binding energies as function of nuclear charge Z. This figure is reproduced from 
Fig. 1 in Ref. 35. 
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3.2. Semi- Classical description 

In order to have further understanding of this phenomenon, we study it in the semi- 
classical scenario. Setting the origin of spherical coordinates (r, 0, <p) at the point-like 
charge, we introduce the vector potential = (A, A$), where A = and A is 
the Coulomb potential. The motion of a relativistic "electron" (scalar particle) with 
mass m and charge e is described by its radial momentum p r , angular momenta p$ 
and the Hamiltonian, 



H ± = ±mc\l + (-PL)2 + (_P*-)2 - V (r), (34) 
V mc mar 

where the potential energy V(r) = eA 0} and ± corresponds for positive and negative 
solutions. The states corresponding to negative energy solutions are fully occupied. 
The angular momentum p^ is conserved, for the Hamiltonian is spherically symmet- 
ric. For a given angular momentum p^, the Hamiltonian (34) describes electron's 
radial motion in the following the effective potential 



E± = ±mc 2 Jl + (^) 2 -V(r). (35) 
V mcr 

The Coulomb potential energy V(r) is given by 

Ze 2 

V{r) = —, (36) 

where Ze 2 = \Qe\. 

In the classical scenario, given different values of angular momenta p<p, the stable 
circulating orbits (states) are determined by the minimum of the effective potential 
E + (r) (35). Using dE + (r)/dr = 0, we obtain the stable orbit location at the radius 
Rl in the unit of the Compton length A = h/mc, 



BlW = Z ^l-(^), (37) 

where a = e 2 /he and p^ > Za. Substituting Eq. (37) into Eq. (35), we find the 
energy of the electron at each stable orbit, 



£(p,) Sm in(E + )=^l-(^L)' 



(38) 



The last stable orbits (minimal energy) are given by 

P4> ^Zah+0+, R L { P4> )^0+, £( P4> )^0+. (39) 

For stable orbits forp^ 3> 1, the radii Rl 3> 1 and energies £ — > mc 2 + 0~; electrons 
in these orbits are critically bound since their banding energy goes to zero. As the 
energy-spectrum (29) (see Eqs. (31,32,33), Eq. (38) shows, only positive or null 
energy solutions (states) exist in the presence of a point-like nucleus. 
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In the semi-classical scenario, the discrete values of angular momentum p^ are 
selected by the Bohr-Sommerfeld quantization rule 

J ' P4> d<P~h{l+ l -) 1 pt{l)~H{l + ±), 1 = 0,1,2,3,... (40) 

describing the semi-classical states of radius and energy 

s ' (,) = W 1 "^) 2 ' <41) 

eu^i-ffi)'. (42) 

Other values of angular momentum p^, radius Rl and energy £ given by Eqs. (37,38) 
in the classical scenario are not allowed. When these semi-classical states are not 
occupied as required by the Pauli Principle, the transition from one state to another 
with different discrete values (Zi , Z2 and Al = I2 — h = ±1) is made by emission 
or absorption of a spin-1 (h) photon. Following the energy and angular-momentum 
conservations, photon emitted or absorbed in the transition have angular momenta 
P<j>(h) — P<j>(h) = — h) = ±fi and energy £{h) — £{h)- As required by the 
Heisenberg indeterminacy principle A0Ap^ ~ A^p^iX) > h, the absolute ground 
state for minimal energy and angular momentum is given by the I = state, ~ 
ft/2, R L ~ ZaX^l - (2Zaf > and £ ~ mc 2 ^/i - {2Za) 2 > for Za < 1/2. 
Thus the stability of all semi-classical states I > is guaranteed by the Pauli 
principle. In contrast for Za > 1/2, there is not an absolute ground state in the 
semi-classical scenario. 

We see now how the lowest energy states are selected by the quantization rule in 
the semi-classical scenario out of the last stable orbits (39) in the classical scenario. 
For the case of Za < 1/2, equating Eq. (39) to p = h(l + 1/2) (40), we find the 
selected state I = is only possible solution so that the ground state I = in 
the semi-classical scenario corresponds to the last stable orbits (39) in the classical 
scenario. On the other hand for the case Za > 1/2, equating Eq. (39) to p^ — 
h(l + l/2) (40), we find the selected state 1 = 1 = (Za—l)/2 > in the semi-classical 
scenario corresponds to the last stable orbits (39) in the classical scenario. This state 
I = I > is not protected by the Heisenberg indeterminacy principle from quantum- 
mechanically decaying in ?i-steps to the states with lower angular momentum and 
energy (correspondingly smaller radius Rl (41)) via photon emissions. This clearly 
shows that the "Z = 137-catastrophe" corresponds to Rl — > 0, falling to the center 
of the Coulomb potential and all semi-classical states (/) are unstable. 

3.3. The critical value of the nuclear charge Z cr — 173. 

A very different situation is encountered when considering the fact the nucleus is 
not point-like and has an extended charge distribution. 36-44 When doing so, the 
Z = 137 catastrophe disappears and the energy-levels £(n,j) of the bound states 
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15*, 2P and 25, ■ ■ ■ smoothly continue to drop toward the negative energy continuum 
as Z increases to values larger than 137, as shown in Fig. 6. The reason is that the 
finite size R of the nucleus charge distribution provides a cutoff for the boundary 
condition at the origin r — > and the energy-levels £(n,j) of the Dirac equation 
are shifted due to this cutoff. In order to determine the critical value Z cr when the 
negative energy continuum (£ < —m e c 2 ) is encountered (see Fig. 6), Zel'dovich and 
Popov 40-44 solved the Dirac equation corresponding to a nucleus of finite extended 
charge distribution, i.e., the Coulomb potential is modified as 

[ r>R 

V(r) = < Z V, N ' (43) 

where R <~ 10~ 12 cm is the size of the nucleus. The form of the cutoff function f(x) 
depends on the distribution of the electric charge over the volume of the nucleus 
(x = r/R,0 < x < 1, with /(l) = 1). Thus, f(x) = (3 - x 2 )/2 corresponds to a 
constant volume density of charge. 

Solving the Dirac equation with the modified Coulomb potential (43) and cal- 
culating the corresponding perturbative shift A£r of the lowest energy level (31) 
one obtains 40 ' 44 



A£ R = m e c 



2 ( Q 2 (2£ e - A ) 2 ^ 
7,(1 + 27.) 



1-27, / f(x)x 2 ^dx 



o 



(44) 



where £ = Za, j z — y/l — £ 2 and A = \n(K/m e cR) >• 1 is a logarithmic parameter 
in the problem under consideration. The asymptotic expressions for the 15*1/2 energy 
that were obtained are 43,44 



v/I - ^cotMAv/i - o < £ < l, 

£(lS 1/2 ) = m e c 2 \ A- 1 , £ = 1, (45) 

VF^Tcot(A v ^T), £>i. 

As a result, the "Z — 137 catastrophe" in Eq. (29) disappears and £ (15*1/2) = 
gives 

Co = l + ^+0(A- 4 ); (46) 

the state 15i/ 2 energy continuously goes down to the negative energy continuum 
since Za > 1, and £(15i/ 2 ) = —1 gives 

«— 1 + 3ArrW +0(A " ) (47) 

as shown in Fig. 6. In Rcf. 40,44 it is found that the critical value ^ c n ^ = Z c a for 
the energy-levels n5i/ 2 and nPi/ 2 to reach the negative energy continuum is equal 
to 

d n) =l + ^+0(A- 3 ). (48) 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



17 

The critical value increases rapidly with increasing n. As a result, it is found that 
Z cr ~ 173 is a critical value at which the lowest energy- level of the bound state lS'1/2 
encounters the negative energy continuum, while other bound states encounter the 
negative energy continuum at Z cr > 173 (see also Ref. 38 for a numerical estimation 
of the same spectrum). We refer the readers to Ref. 40-45 for mathematical and 
numerical details. 

When Z > Z cr = 173, the lowest energy-level of the bound state lS'1/2 enters 
the negative energy continuum. Its energy-level can be estimated as follows 

£(lS 1/2 ) = m e c 2 - ^ < -m e c 2 , (49) 

where f is the average radius of the lS'1/2 state's orbit, and the binding energy 
of this state satisfies Za/f > 2m e c 2 . If this bound state is unoccupied, the bare 
nucleus gains a binding energy Za/f larger than 2m e c 2 , and becomes unstable 
against the production of an electron-positron pair. Assuming this pair-production 
occurs around the radius f, we have energies for the electron (e_) and positron (e + ) 
given by 

e_ = v /| C p_|2 +m 2 c 4_^. e+ = v /|cp + | 2 + m2 c 4 + ^, (50) 

where p± are electron and positron momenta, and p_ = p + . The total energy 
required for the production of a pair is 

e-+=e-+e+ = 2 V / |cp-| 2 +m2 c 4, (51) 

which is independent of the potential V(f). The potential energies ±eV(f) of the 
electron and positron cancel each other out and do not contribute to the total energy 
(51) required for pair production. This energy (51) is acquired from the binding 
energy (Za/f > 2m e c 2 ) by the electron filling into the bound state lSi/ 2 . A part 
of the binding energy becomes the kinetic energy of positron that goes out. This is 
analogous to the familiar case that a proton [Z = 1) catches an electron into the 
ground state 1Si/ 2 , and a photon is emitted with the energy not less than 13.6 eV. 
In the same way, more electron-positron pairs are produced, when Z >• Z cr = 173 
the energy-levels of the next bound states 2P 1 / 2 , 2S 3 / 2 , • • • enter the negative energy 
continuum, provided these bound states of the bare nucleus are unoccupied. 

3.4. Positron production 

Ref. 46,47 proposed that when Z > Z cr , the bare nucleus spontaneously produces 
pairs of electrons and positrons: the two positrons a go off to infinity and the effective 
charge of the bare nucleus decreases by two electrons, which corresponds exactly to 
filling the K-shcll. b A more detailed investigation was made for the solution of the 



a Hypcrfine structure of lS'1/2 state: single and triplet. 

b The supposition was made in Ref. 46,47 that the electron density of lS'1/2 state, as well as the 
vacuum polarization density, is dclocalizcd at Z — » Z cr . Later it was proven to be incorrect. 41,42,44 
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Dirac equation at Z <~ Z cr , when the lowest electron level IS1/2 merges with the 
negative energy continuum, by Ref. 40-43,48. This further clarified the situation, 
showing that at Z > Z cr , an imaginary resonance energy of Dirac equation appears 

e = e -i^, (52) 

where 

e = -m e c 2 - a(Z - Z cr ), (53) 

r ~ 6{Z - Z cr ) cxp (-byJ-^-^J , (54) 

and a, b are constants, depending on the cutoff A (for example, b — 1.73 for Z = 
Z cr = 173, see Ref. 41,42,44). The energy and momentum of the emitted positrons 
are |e | and |p| = y/\e \ - m e c 2 . 

The kinetic energy of the two positrons at infinity is given by 

£ P = je | — m e c 2 = a(Z - Z cr ) + . . . , (55) 

which is proportional to Z — Z cr (as long as (Z — Z cr ) <C Z cr ) and tends to zero 
as Z — > Z cr . The pair-production resonance at the energy (52) is extremely nar- 
row and practically all positrons are emitted with almost same kinetic energy for 
Z <~ Z cr , i.e., nearly monoenergetic spectra (sharp line structure). Apart from a 
pre-exponential factor, T in Eq. (54) coincides with the probability of positron pro- 
duction, i.e., the penetrability of the Coulomb barrier. The related problems of 
vacuum charge density due to electrons filling into the K-shcll and charge renor- 
malization due to the change of wave function of electron states are discussed by 
Ref. 49-53. An extensive and detailed review on this theoretical issue can be found 
in Ref. 35,44,45,54. 

On the other hand, some theoretical work has been done studying the possibility 
that pair production due to bound states encountering the negative energy contin- 
uum is prevented from occurring by higher order processes of quantum field theory, 
such as charge renormalization, electron self-energy and nonlinearities in electrody- 
namics and even the Dirac field itself. 55-61 However, these studies show that various 
effects modify Z cr by a few percent, but have no way to prevent the binding energy 
from increasing to 2m e c 2 as Z increases, without simultaneously contradicting the 
existing precise experimental data on stable atoms. 62 

Following this, special attention has been given to understand the process of 
creating a nucleus with Z > Z cr by collision of two nuclei of charge Z\ and Zi such 
that Z = Z\ + Z2 > Z cr . 46 ~ 48,63 ' 64 To observe the emission of positrons coming 
from pair production occurring near an overcritical nucleus temporarily formed by 
two nuclei, the following necessary conditions have to be fullfillcd: (i) the atomic 
number of an overcritical nucleus must be larger than Z cr = 173; (ii) the lifetime 
(the sticking time of two-nuclear collisions) of the overcritical nucleus must be much 
longer than the characteristic time (h/m e c 2 ) of pair production; (iii) inner shells 
(K-shell) of the overcritical nucleus should be unoccupied. 
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When in the course of a heavy- ion collision the two nuclei come into contact, 
some deep- inelastic reactions have been claimed to exist for a certain time At s . The 
duration At s of this contact (sticking time) is expected to depend on the nuclei 
involved in the reaction and on beam energy. For very heavy nuclei, the Coulomb 
interaction is the dominant force between the nuclei, so that the sticking times At s 
are typically much shorter and on the average probably do not exceed 1 ~ 2 • 1CP 21 
sec. 35 Accordingly the calculations (see Fig. 7) also show that the time when the 
binding energy is overcritical is very short, about 1.2-10 -21 sec. Theoretical attempts 
have been proposed to study the nuclear aspects of heavy-ion collisions at energies 
very close to the Coulomb barrier and search for conditions, which would serve as 
a trigger for prolonged nuclear reaction times, (the sticking time At s ) to enhance 
the amplitude of pair production. 35 ' 62, 65-67 Up to now no conclusive theoretical or 
experimental evidence exists for long nuclear delay times in very heavy collision 
systems. 




Fig. 7. Energy expectation values of the lsa state in a U+U collision at 10 GeV/nucleon. The 
unit of time is h/m e c 2 . This figure is reproduced from Fig. 4 in Ref. 35. 



It is worth noting that several other dynamical processes contribute to the pro- 
duction of positrons in undercritical as well as in overcritical collision systems. 55-58 
Due to the time-energy uncertainty relation (collision broadening), the energy- 
spectrum of such positrons has a rather broad and oscillating structure, consider- 
ably different from a sharp line structure that we would expect from pair-production 
positron emission alone. 
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3.5. Experiments 

As already remarked, if the sticking time At s could be prolonged, the probability 
of pair production in vacuum around the super heavy nucleus would be enhanced. 
As a consequence, the spectrum of emitted positrons is expected to develop a sharp 
line structure, indicating the spontaneous vacuum decay caused by the overcritical 
electric field of a forming super heavy nuclear system with Z > Z cr . If the sticking 
time At s is not long enough and the sharp line of pair production positrons has not 
yet well-developed, in the observed positron spectrum it is difficult to distinguish 
the pair production positrons from positrons created through other different mecha- 
nisms. Prolonging the sticking time and identifying pair production positrons among 
all other particles 68,69 created in the collision process are important experimental 
tasks. 70 " 77 

For nearly 20 years the study of atomic excitation processes and in particular 
of positron creation in heavy-ion collisions has been a major research topic at GSI 
(Darmstadt). 78-82 The Orange and Epos groups at GSI (Darmstadt) discovered 
narrow line structures (see Fig. 8) of unexplained origin, first in the single positron 
energy spectra and later in coincident electron-positron pair emission. Studying 
more collision systems with a wider range of the combined nuclear charge Z = 
Z\ + Z2 they found that narrow line structures are essentially independent of Z . 
This rules out the explanation of a pair-production positron, since the line would 
be expected at the position of the lsa resonance, i.e., at a kinetic energy given 
by Eq. (55), which is strongly Z dependent. Attempts to link this positron line to 
spontaneous pair production have failed. Other attempts to explain this positron 
line in term of atomic physics and new particle scenario were not successful as well. 35 

The anomalous positron line problem has perplexed experimentalists and the- 
orists alike for more than a decade. Moreover, later results obtained by the Apex 
collaboration at Argonne National Laboratory showed no statistically significant 
positron line structures. 83,84 This is in strong contradiction with the former results 
obtained by the Orange and Epos groups. However, the analysis of Apex data was 
challenged in the comment by Ref. 85,86 for the Apex measurement would have 
been less sensitive to extremely narrow positron lines. A new generation of experi- 
ments (Apex at Argonne and the new Epos and Orange setups at GSI) with much 
improved counting statistics has failed to reproduce the earlier results. 35 

To overcome the problem posed by the short time scale of pair production (1CP 21 
sec), hopes rest on the idea to select collision systems in which a nuclear reaction 
with sufficient sticking time occurs. Whether such a situation can be realized still 
is an open question. 35 In addition, the anomalous positron line problem and its 
experimental contradiction overshadow on the field of studying the pair production 
in heavy ion collisions. In summary, clear experimental signals for electron-positron 
pair production in heavy ion collisions are still missing 35 at the present time. 
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Fig. 8. Two typical example of coincident electron-positron spectra measured by the Epose group 
in the system U+Th (left) and by the Orange group in U+Pb collisions (right). When plotted 
as a function of the total energy of the electron and positron, very narrow line structures were 
observed. This figure is reproduced from Fig. 7 in Ref. 35. 

4. The extraction of blackholic energy from a black hole by 
vacuum polarization processes 

We recall here the basic steps leading to the study of a critical electric field in a Kerr- 
Newman black hole. We recall the theoretical framework to apply the Schwinger 
process in general relativity in the field of a Kerr-Newman geometry as well as 
the process of extraction of the "blackholic" energy. We then recall the basic con- 
cepts of dyadosphere and dyadotorus leading to a photon-electron-positron plasma 
surrounding the black hole. 

4.1. The mass-energy formula of black holes 

The same effective potential technique (see Landau and Lifshitz 87 ) which allowed 
the analysis of circular orbits around the black hole was crucial in reaching the 
equally interesting discovery of the reversible and irreversible transformations of 
black holes by Christodoulou and Ruffini, 88 which in turn led to the mass-energy 
formula of the black hole 



with 




where 

S = inpl = 4n(rl + = 16rr (^j M? , (58) 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



22 



is the horizon surface area, M- u is the irreducible mass, r + is the horizon radius 
and p + is the quasi-spheroidal cylindrical coordinate of the horizon evaluated at 
the equatorial plane. Extreme black holes satisfy the equality in Eq. (57). 

From Eq. (56) follows that the total energy of the black hole Ebh can be split 
into three different parts: rest mass, Coulomb energy and rotational energy. In prin- 
ciple both Coulomb energy and rotational energy can be extracted from the black 
hole (Christodoulou and Ruffini 88 ). The maximum extractable rotational energy 
is 29% and the maximum extractable Coulomb energy is 50% of the total energy, 
as clearly follows from the upper limit for the existence of a black hole, given 
by Eq. (57). We refer in the following to both these extractable energies as the 
blackholic energy. We outline how the extraction of the blackholic energy is indeed 
made possible by electron-positron pair creation. We also introduce the concept of a 
"dyadosphere" or "dyadotorus" around a black hole and we will the outline how the 
evolution of the electron-positron plasma will naturally lead to a self-acceleration 
process creating the very high Lorentz gamma factor observed in GRBs. 

4.2. Vacuum polarization in Kerr-Newman geometries 

We already discussed the phenomenon of electron-positron pair production in a 
strong electric field in a flat space-time. We study the same phenomenon occurring 
around a black hole with electromagnetic structure (EMBH). For simplicity and in 
order to give the fundamental energetic estimates, we postulate that the collapse 
has already occurred and has led to the formation of an EMBH. Clearly, this is 
done only in order to give an estimate of the transient phase occurring during the 
gravitational collapse. In reality, an EMBH will never be formed because the vacuum 
polarization process will carry away its electromagnetic energy during the last phase 
of gravitational collapse. Indeed being interested in this transient phenomenon, we 
can estimate its energetics by the conceptual analysis of an already formed EMBH, 
and this is certainly valid for the estimate of the energetics which we will encounter 
in a realistic phase of gravitational collapse. 

The spacetime around the EMBH is described by the Kerr-Newman geome- 
try whose metric we rewrite here for convenience in Boycr-Lindquist coordinates 



where A = r 2 — 2Mr + a 2 + Q 2 and £ = r 2 + a 2 cos 2 8, as before and as usual M is 
the mass, Q the charge and a the angular momentum per unit mass of the EMBH. 
We recall that the Reissner-Nordstr0m geometry is the particular case a = of 
a non-rotating black hole. Natural units G = h = c = 1 will be adopted in this 
section. 

The electromagnetic vector potential around the Kerr-Newman black hole is 
given in Boyer-Lindquist coordinates by 




A = -gS _1 r(rfi - a sin 2 9d<t>). 



(60) 
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The electromagnetic field tensor is then 

F = dA = 2Q£~ 2 [(r 2 - a 2 cos 2 9)dr A dt — 2a 2 r cos 9 sin OdO A dt 

- a sin 2 9(r 2 - a 2 cos 2 6)dr Ad<j> + 2ar(r 2 + a 2 ) cos 9 sin 9d9 A #]. (61) 

After some preliminary work in Refs. 89-91, the occurrence of pair production 
in a Kerr-Newman geometry was addressed by Deruelle. 92 In a Reissner-Nordstrom 
geometry, QED pair production has been studied by Zaumen 93 and Gibbons. 94 The 
corresponding problem of QED pair production in the Kerr-Newman geometry was 
addressed by Damour and Ruffini, 95 who obtained the rate of pair production with 
particular emphasis on: 

• the limitations imposed by pair production on the strength of the electro- 
magnetic field of a black hole; 96 

• the efficiency of extracting rotational and Coulomb energy (the "blackholic" 
energy) from a black hole by pair production; 

• the possibility of having observational consequences of astrophysical inter- 
est. 

The third point was in fact a far-reaching prevision of possible energy sources for 
gamma ray bursts that are now one of the most important phenomena under current 
theoretical and observational study. In the following, we recall the main results of 
the work by Damour and Ruffini. 

In order to study the pair production in the Kerr-Newman geometry, they in- 
troduced at each event (t, r, 9, a local Lorentz frame associated with a stationary 
observer O at the event (t,r,9,cp). A convenient frame is defined by the following 
orthogonal tetrad 97 

u(°' = (A/Z) 1 / 2 (dt-asm 2 9d ( j)), (62) 
w« = (S/A) 1 / 2 ^, (63) 
u/ 2 > = Y}' 2 d9, (64) 
u/ 3 > = sinc?E- 1 / 2 ((r 2 + a 2 )dcp - adt). (65) 

In this Lorentz frame, the electric potential Aq, the electric field E and the magnetic 
field B are given by the following formulas (c.e.g. Ref. 98), 

E a = J° ] F af} , 

We then obtain 



-Qr(EA)- 1 / 2 , 



(66) 
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while the electromagnetic fields E and B are parallel to the direction of a/ 1 ) and 
have strengths given by 

= Q5T 2 (r 2 - a 2 cos 2 9), (67) 
B (1) = Q5T 2 2ar cos 6>, (68) 

respectively. The maximal strength i? max of the electric field is obtained in the case 
a = at the horizon of the EMBH: r = r + . We have 

£ max - Q 2 /r\. (69) 

Equating the maximal electric field strength (69) to the critical value (28), one 
obtains the maximal black hole mass M max ~ 7.2 • 10 6 M Q for pair production to 
occur. For any black hole with mass smaller than M max , the pair production process 
can drastically modify its electromagnetic structure. 

Both the gravitational and the electromagnetic background fields of the Kerr- 
Newman black hole are stationary when considering the quantum field of the elec- 
tron, which has mass m e and charge e. If m e M 3> 1, then the spatial variation scale 
GM/c 2 of the background fields is much larger than the typical wavelength H/m e c 
of the quantum field. As far as purely QED phenomena such as pair production 
are concerned, it is possible to consider the electric and magnetic fields defined by 
Eqs. (67,68) as constants in the neighborhood of a few wavelengths around any 
events (r, 8, <p, t). Thus, the analysis and discussion on the Sauter-Euler-Heisenberg- 
Schwinger process over a flat space-time can be locally applied to the case of the 
curved Kerr-Newman geometry, based on the equivalence principle. 

The rate of pair production around a Kerr-Newman black hole can be obtained 
from the Schwingcr formula (26) for parallel electromagnetic fields e = E^ and 
(3 = as: 

f e 2 E (1) B (1) 1 / mrB (1) \ ( nirE< 



,,^ h UrM-^7> (70) 



V At: 2 

The total number of pairs produced in a region D of the space-time is 

N = Jd<x^^, (71) 

where y/—g = Xsin#. In Ref. 95, it was assumed that for each created pair the 
particle (or antiparticle) with the same sign of charge as the black hole is expelled 
to infinity with charge e, energy w and angular momentum while the antiparticle 
is absorbed by the black hole. This implies the decrease of charge, mass and angular 
momentum of the black hole and a corresponding extraction of all three quantities. 
The rates of change of the three quantities are then determined by the rate of pair 
production (70) and by the conservation laws of total charge, energy and angular 
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momentum 

Q = -Re, 

M = -R(w), (72) 
L = -R(l<t,), 

where R = N is the rate of pair production and (w) and (1$) represent some suitable 
mean values for the energy and angular momentum carried by the pairs. 

Supposing the maximal variation of black hole charge to be AQ = —Q, one 
can estimate the maximal number of pairs created and the maximal mass-energy 
variation. It was concluded in Ref. 95 that the maximal mass-energy variation in 
the pair production process is larger than 10 erg and up to 10 58 erg, depending on 
the black hole mass. They concluded at the time "this work naturally leads to a 
most simple model for the explanation of the recently discovered j-ray bursts". 



4.3. The "Dyadosphere" 

We first recall the three theoretical results which provide the foundation of the 
EMBH theory. 

In 1971 in the article "Introducing the Black Hole", 16 the theorem was advanced 
that the most general black hole is characterized uniquely by three independent 
parameters: the mass-energy M, the angular momentum L and the charge Q making 
it an EMBH. Such an ansatz, which came to be known as the "uniqueness theorem" 
has turned out to be one of the most difficult theorems to be proven in all of physics 
and mathematics. The progress in the proof has been authoritatively summarized 
by Ref. 99. The situation can be considered satisfactory from the point of view 
of the physical and astrophysical considerations. Nevertheless some fundamental 
mathematical and physical issues concerning the most general perturbation analysis 
of an EMBH are still the topic of active scientific discussion. 100 

In 1971 Christodoulou and Ruffini 88 obtained the mass-energy formula of a 
Kerr-Newman black hole, see Eqs. (56-58). As we just recalled, in 1975 there was 
the Damour and Ruffini work 95 on the vacuum polarization of a Kerr-Newman 
geometry. The key point of this work was the possibility that energy on the order 
of 10 54 ergs could be released almost instantaneously by the vacuum polarization 
process of a black hole. At the time, however, nothing was known about the distance 
of GRB sources or their energetics. The number of theories trying to explain GRBs 
abounded, as mentioned by Ruffini in the Kleinert Festschrift. 101 

After the discovery in 1997 of the afterglow of GRBs 102 and the determination of 
the cosmological distance of their sources, we noticed the coincidence between their 
observed energetics and the one theoretically predicted by Damour and Ruffini. 95 
We therefore returned to these theoretical results with renewed interest developing 
some additional basic theoretical concepts 103-107 such as the dyadosphere and, more 
recently, the dyadotorus. 
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As a first simplifying assumption we have developed our considerations in the 
absence of rotation using spherically symmetric models. The space-time is then 
described by the Reissner-Nordstrom geometry, whose spherically symmetric metric 
is given by 

d 2 s = g tt {r)d 2 t + g rr (r)d 2 r + r 2 d 2 9 + r 2 sin 2 dd 2 (f> , 



(73) 



where g u (r) = - 



j _ 2GM , Q 2 G 



= —a 2 (r) and g r r{r) = a 2 {r). 
The first result we obtained is that the pair creation process does not occur at 
the horizon of the EMBH: it extends over the entire region outside the horizon in 
which the electric field exceeds the value E* of the order of magnitude of the critical 
value given by Eq. (28). The pair creation process, as we recalled in the previous 
sections (see e.g. Fig. 5), is a quantum tunnelling between the positive and negative 
energy states, which needs only a level crossing but can occur, of course, also for 
E* < E c if the extent of the field is large enough, although with decreasing intensity. 
Such a process occurs also for E* > E c . In order to give a scale of the phenomenon, 
and for dcfinitencss, in Ref. 105 we first consider the case of E* = E c . Since the 
electric field in the Reissner-Nordstrom geometry has only a radial component given 
by 108 



E M = ^2 ' 

this region extends from the horizon radius 

r+ = 1.47 • 10 5 a*(1 + v 7 ! -C 2 ) cm 

out to an outer radius 103 

h_\ * fGM\ 5 /mp\ 5 /_e \ 3 /_ (} 
mc ) \ c 2 J 



(74) 



(75) 



Vgm 



1.12- 10 8 y^cm, (76) 



where we have introduced the dimensionless mass and charge parameters \x = -j^-, 

« = A * 1; see Fig - 10 - 

The second result has been to realize that the local number density of electron 
and positron pairs created in this region as a function of radius is given by 

Q 



n e + e -(r) 



4-Kr 2 



\ mc ) 



(r*) 



(77) 



and consequently the total number of electron and positron pairs in this region is 



N°+ 



Q-Qc 



1 + 



(r* — r + ) 



h 

mc 



(78) 



where Q c = E c r\. 

The total number of pairs is larger by an enormous factor r* / (h/mc) > 10 18 
than the value Q/e which a naive estimate of the discharge of the EMBH would 
have predicted. Due to this enormous amplification factor in the number of pairs 
created, the region between the horizon and r* is dominated by an essentially high 
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Fig. 9. The energy extracted by the process of vacuum polarization is plotted (solid lines) as 
a function of the mass M in solar mass units for selected values of the charge parameter £ = 
1, 0.1, 0.01 (from top to bottom) for an EMBH, the case £ = 1 reachable only as a limiting process. 
For comparison we have also plotted the maximum energy extractable from an EMBH (dotted 
lines) given by Eq. (56). Details in Ref. 109. 

density neutral plasma of electron-positron pairs. We have defined this region as the 
dyadosphere of the EMBH from the Greek duas, duadsos for pairs. Consequently 
we have called r* the dyadosphere radius r* = rd s . 103-105 The vacuum polarization 
process occurs as if the entire dyadosphere is subdivided into a concentric set of 
shells of capacitors each of thickness h/m e c and each producing a number of e+e~ 
pairs on the order of <~ Q/e (see Fig. 10). The energy density of the electron-positron 
pairs is given by 



(79) 
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+Q -Q 




mc 





plasma (e + e~ y) 

Fig. 10. The dyadosphere of a Reissner-Nordstrom black hole can be represented as equivalent to 
a concentric set of capacitor shells, each one of thickness h/m e c and producing a number of e~*~e — 
pairs of the order of ~ Q/e on a time scale of 10 — 21 s, where Q is the EMBH charge. The shells 
extend in a region of thickness Ar, from the horizon r+ out to the dyadosphere outer radius r^s 
(see text). The system evolves to a thermalised plasma configuration. 



(see Figs. 2-3 of Ref. 104). The total energy of pairs converted from the static 
electric energy and deposited within the dyadosphere is then 

4" 



Edya = ( 1 — 

2r+ V r ds 



rds 



(80) 



As we will see in the following this is one of the two fundamental parameters of 
the EMBH theory (see Fig. 11). In the limit £±- ->■ 0, Eq. (80) leads to E dya -> ±2l, 
which coincides with the energy extractable from EMBHs by reversible processes 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



29 



(Mi r = const.), namely Ebh — M ir = ^f-, see Fig. 9. Due to the very large 
pair density given by Eq. (77) and to the sizes of the cross-sections for the process 
e + e~ <-> 7 + 7, the system has been assumed to thermalize to a plasma configuration 
for which 



n e + = n e 



(81) 



where n° +e _ is the total number density of e+e~-pairs created in the dyado- 
sphere. 104 ' 105 This assumption has been in the meantime rigorously proven by Ak- 
senov, Ruffini and Vereshchagin. 110 

The third result which we have introduced again for simplicity is that for a 
given Edya we have assumed either a constant average energy density over the entire 
dyadosphere volume, or a more compact configuration with energy density equal to 
its peak value. These are the two possible initial conditions for the evolution of the 
dyadosphere (see Fig. 11). 
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Fig. 11. Left) Selected lines corresponding to fixed values of the Edya are given as a function 
of the two parameters fi £, only the solutions below the continuous heavy line are physically rele- 
vant. The configurations above the continuous heavy lines correspond to unphysical solutions with 
r ds < r +- Right) Two different approximations for the energy density profile inside the dyado- 
sphere. The first one (dashed line) fixes the energy density equal to its peak value, and computes 
an "effective" dyadosphere radius accordingly. The second one (dotted line) fixes the dyadosphere 
radius to its correct value, and assumes an uniform energy density over the dyadosphere volume. 
The total energy in the dyadosphere is of course the same in both cases. The solid curve represents 
the real energy density profile. 



The above theoretical results permit a good estimate of the general energetics 
processes originating in the dyadosphere, assuming an already formed EMBH. In 
reality, if the GRB data become accurate enough, the full dynamical description of 
the dyadosphere formation will be needed in order to explain the observational data 
by the general relativistic effects and characteristic time scales of the approach to 
the EMBH horizon. 111 - 114 
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4.4. The "Dyadotorus" 

We turn now to examine how the presence of rotation modifies the geometry of the 
surface containing the region where electron-positron pairs are created as well as 
the conditions for the existence of such a surface. Due to the axial symmetry of the 
problem, we have called this region the "dyadotorus". 115 

Following Damour 116,117 we introduce at each point of the spacetime the orthog- 
onal Carter tetrad 

w(°) = {A/^) 1/2 {dt-asm 2 6dip), 

w « -(s/A) 1 / 2 ^, jv = ^i 2 de, 

= sintfE- 1 / 2 ^ 2 + a 2 )dip - adt] , 

where £ = r 2 + a 2 cos 2 9 and A = r 2 - 2Mr + a 2 + Q 2 , M being the mass, Q the 
electric charge and a the angular momentum per unit mass of the black hole. Thus, 
in the Lorentz frame defined by the above tetrad the Kerr-Newman metric reads 

ds 2 = -(^W) 2 + (c^) 2 + (c^ 2 )) 2 + (c^) 2 , (82) 

and the outer horizon of the black hole is localized at r + — M + y/ M 2 — (a 2 + Q 2 ). 
In this frame the electric and magnetic field obtained from the potential 

A = -gr(SA)- 1 / 2 w<°) , (83) 

are parallel, i.e., 

E (1) = Q£~ 2 (r 2 - a 2 cos 2 6) , (84) 

= 2aQY>- 2 rcos6 . (85) 

The rate R of pair creation using the Schwinger approach is 95 




where 

oo 

2Im£ = (47r)- 1 (S (1) e/7r) 2 ]T n- 2 {nvB {1) / E (1) ) 

n=l 

coth(n7ri?( 1 )/i?( 1 )) exp(-n7rm 2 /_E( 1 )) , (87) 

where e = m 2 / E c . 

We define the dyadotorus by the condition |E| = kE c , where 10 _1 < k < 10. 
Solving for r and introducing the dimensionless quantities £ = Q/M, a = a/M, 
H = M/M Ql £ = kE c M q w 1.873 x 1(T 6 and r = r/M (with M Q w 1.477 x 10 5 
cm) we get 

(—) = - " 2 cos2 0±J - ^l" 2 cos 2 ^ , (88) 

where the ± signs correspond to the two different parts of the surface. 
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The two parts of the surface join at the particular values 8* and tt — 8* of the 
polar angle where 8* — arccos ( 2% /2a \f^s) • ^ nc rc< 4 u i rcrncn t that cos6>* < 1 can 
be solved for instance for the charge parameter £, giving a range of values of £ for 
which the dyadotorus takes one of the shapes (see fig. 12) 

surface I 01111580 " * ^ C * (89) 
y thorus-like if £ < £* 

where £* = 8/i£a 2 . 

In Fig. 12 we show some examples of the dyadotorus geometry for different sets 
of parameters for an extreme Kerr-Newman black hole (M 2 = a 2 + Q 2 ), we can 
see the transition from a toroidal geometry to an ellipsoidal one depending on the 
value of the black hole charge. 

Fig. 13 shows instead the projections of the surfaces corresponding to different 
values of the ratio |E|/_E C = k for the same choice of parameters as in Fig. 12 (b), as 
an example. We see that the region enclosed by such surfaces shrinks for increasing 
values of n. 



5. On the observability of electron-positron pairs created by 
vacuum polarization in Earth-bound experiments and in 
astrophysics 

In summary, from the considerations we have presented in the previous sections 
three different earth-bound experiments and one astrophysical observation have 
been proposed for identifying the polarization of the electronic vacuum due to 
a supercritical electric field postulated by Sauter-Heisenberg-Euler-Schwinger (see 
Ref. 11,13,17,118): 

(1) In collisions of heavy ions near the Coulomb barrier, as first proposed in Ref. 46, 
47 (see also Ref. 44,48,119). Despite some apparently encouraging results (see 
Ref. 120), such efforts have failed so far due to the small contact time of the 
colliding i ons / 9 > 8 °. 8 2,83,i2i Typically the electromagnetic energy involved in the 
collisions of heavy ions with impact parameter Zi ~ 10~ 12 cm is E\ ~ 10~ 6 erg 
and the lifetime of the diatomic system is t\ ~ 10~ 22 s. 

(2) In collisions of an electron beam with optical laser pulses: a signal of positrons 
above background has been observed in collisions of a 46.6 GeV electron beam 
with terawatt pulses of optical laser in an experiment at the Final Focus Test 
Beam at SLAC; 122 it is not clear if this experimental result is an evidence for 
the vacuum polarization phenomenon. The energy of the laser pulses was Ei ~ 
10 7 erg, concentrated in a space-time region of spacial linear extension (focal 
length) I2 ~ 10~ 3 cm and temporal extension (pulse duration) ti <~ 10~ 12 s. 122 

(3) At the focus of an X-ray free electron laser (XFEL) (see Ref. 123-125 and 
references therein). Proposals for this experiment exist at the TESLA collider 
at DESY and at the LCLS facility at SLAC. 123 Typically the electromagnetic 
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Fig. 12. The projection of the dyadotorus on the X — Z plane (X = rsin#, Z = rcosB are 
Cartesian-like coordinates built up simply using the Boyer-Lindquist radial and angular coordi- 
nates) is shown for an extreme Kerr-Newman black hole with fi = 10 and different values of the 
charge parameter £ = [1,1.3,1.49,1.65] X 10~ 4 (from (a) to (d) respectively). The black circle 
represents the black hole horizon. 



energy at the focus of an XFEL can be E3 ~ 10 6 erg, concentrated in a space- 
time region of spacial linear extension (spot radius) I3 ~ 10~ 8 cm and temporal 
extension (coherent spike length) £3 ~ 10 _13 s. 123 

and from astrophysics 

(1) Around an electromagnetic black hole (black hole), 95 ' ' 105 giving rise to the 
observed phenomenon of GRBs. 126-129 The electromagnetic energy of an black 
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Fig. 13. The projections of the surfaces corresponding to different values of the ratio |E|/B C = k 
are shown for the same choice of parameters as in Fig. 12 (b), as an example. The gray shaded 
region is part of the "dyadotorus" corresponding to the case k = 1 as plotted in Fig. 12 (b). The 
region delimited by dashed curves corresponds to K = 0.8, i.e., to a value of the strength of the 
electric field smaller than the critical one, and contains the dyadotorus; the latter in turn contains 
the white region corresponding to k = 1.4, i.e., to a value of the strength of the electric field greater 
than the critical one. 



hole of mass M — 10M© and charge Q ~ 0.1M/VG is E4 ~ 10 54 erg and it is 
deposited in a space-time region of spacial linear extension I4 ~ 10 8 cm 105 ' 112 
and temporal extension (collapse time) t± <~ 10~ 2 s. 130 

As we will see in the following, the creation of an electron-positron plasma is 
indeed essential to explain not only the GRB energetics but also the unprecedentedly 
large Lorentz gamma factor observed in GRBs. 



6. Electrodynamics for nuclear matter in bulk 

We have seen how critical fields may be generated in heavy ion collisions, in colli- 
sions of electron beams with optical laser pulses and in X-ray free electron lasers. 
The explanation of GRBs leads to a theoretical framework postulating the existence 
of critical and overcritical fields in black holes, in order to extract their blackholic 
energy. It becomes natural, then, to ask if there is any mechanism which can lead to 
the existence of a critical field not only on the above mentioned microscopic scale, 
but also on macroscopic scales, possibly to be encountered at the onset of the pro- 
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cess of gravitational collapse. It is clear that such processes should be common to a 
large variety of initial conditions occurring in the gravitational collapse either of a 
neutron star or of two binary neutron stars or again of a binary system formed by 
a neutron star and a white dwarf or, finally, in the general formation of an interme- 
diate mass black hole. It is already clear from the work of Damour and Ruffini that 
the gravitational collapse giving birth to the black holes in active galactic nuclei 
with masses much larger that 1O 6 M cannot give rise to instantaneous vacuum po- 
larization processes leading to GRBs. For this reason we present here an alternative 
treatment of the electrodynamics for nuclear matter in bulk, presenting for the first 
time a unified approach to neutron star physics and nuclear physics, covering the 
range of atomic number from A <~ 10 2 to A <~ 10 57 . 



6.1. The Thomas- Fermi equations for heavy ions 

It is well know that the Thomas-Fermi equation is the exact theory for atoms, 
molecules and solids as Z — > oo. 131 We show in this section that the relativistic 
Thomas-Fermi theory developed to study atoms for heavy nuclei with Z ~ 10 636 - 132 
gives important basic new information about the state of nuclear matter in bulk in 
the limit of N ~ (?npi anc k/m n ) 3 nucleons of mass m n and about its electrodynamic 
properties. 

The analysis of bulk nuclear matter in neutron stars composed of a degenerate 
gas of neutrons, protons and electrons has traditionally been approached by micro- 
scopically implementing the charge neutrality condition by requiring the electron 
density n e (x) to coincide with the proton density n p (x) 

n e (x) = n p (x). (90) 

It is clear, however, that especially when conditions close to gravitational collapse 
occur, there is an ultra-relativistic component of degenerate electrons whose con- 
finement requires the existence of very strong electromagnetic fields, in order to 
guarantee the overall charge neutrality of the neutron star. Under these conditions 
Eq. (90) will necessarily be violated. We will show here that they will develop electric 
fields close to the critical value E c introduced by Ref. 11,13,17-19: 

E , . »£. m 

Special attention to the existence of critical electric fields and the possible 
condition for electron-positron (e + e~) pair creation out of the vacuum in the 
case of heavy bare nuclei with the atomic number Z > 173 has been given by 
Ref. 36,42,44,47,62,133. They analyzed the specific pair creation process of an 
electron-positron pair around both a point-like and extended bare nucleus by direct 
integration of Dirac equation. These considerations have been extrapolated to much 
heavier nuclei Z ^> 1600, implying the creation of a large number of e + e~ pairs by 
using a statistical approach based on the relativistic Thomas-Fermi equation by 
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Rcf. 134,135. Using substantially the same statistical approach based on the rela- 
tivists Thomas-Fermi equation, Ref. 132,136 have analyzed the electron densities 
around an extended nucleus in a neutral atom all the way up to Z ~ 6000. They 
have shown the effect of the penetration of the electron orbitals well inside the nu- 
cleus, leading to a screening of the nuclei positive charge and to the concept of an 
"effective" nuclear charge distribution. 

All of this work assumes for the radius of the extended nucleus the semi-empirical 
formula, 137 

R c w r a A 1/3 , r = 1.2 • 10" 13 cm, (92) 

where the mass number A = N n + N p , N n and N p are the neutron and proton 
numbers. The approximate relation between A and the atomic number Z — N p 

Z~\ (93) 
was adopted in Ref. 134,135, or the empirical formula 

was adopted in Ref. 132,136. 



6.2. Electroweak equilibrium in Nuclear Matter in Bulk 

We outline an alternative approach of the description of nuclear matter in bulk: 
it generalizes the above treatments, already developed and tested for the study of 
heavy nuclei, to the case of N ~ (mp\a, nc ^/m n ) 3 nucleons. This more general ap- 
proach differs in many aspects from the ones in the current literature and reproduces 
the above treatments in the limiting case of A smaller than 10 6 ,. We will look for a 
solution implementing the condition of overall charge neutrality of the star as given 
by 

N e = N p , (95) 

which significantly modifies Eq. (90), since now N e (N p ) is the total number of 
electrons (protons) of the equilibrium configuration. 

We present here only a simplified prototype of this approach, outlining the es- 
sential relative role of the four fundamental interactions present in the neutron star 
physics: the gravitational, weak, strong and electromagnetic interactions. In addi- 
tion, we also implement the fundamental role of Fermi-Dirac statistics and the phase 
space blocking due to the Pauli principle in the degenerate configuration. The new 
results essentially depend from the coordinated action of the above five theoretical 
components and cannot be obtained if any one of them is neglected. 

Let us first recall the role of gravity. In the case of neutron stars, unlike the case 
of nuclei where its effects can be neglected, gravitation has the fundamental role of 
defining the basic parameters of the equilibrium configuration. As pointed out by 
Ref. 23 at a Newtonian level and by Ref. 138 in general relativity, configurations of 
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equilibrium exist at approximately one solar mass and at an average density around 
the nuclear density. This result is obtainable considering only the gravitational 
interaction of a system of Fermi degenerate self-gravitating neutrons, neglecting all 
other particles and interactions. This situation can be formulated within a Thomas- 
Fermi self-gravitating model (sec e.g. Ref. 139). 

In the present case of our simplified prototype model directed at revealing new 
clectrodynamic properties, the role of gravity is taken into account simply by con- 
sidering in line with the generalization of the above results a mass-radius relation 
for the baryonic core 

R NS = R c n H mpia " ck . (96) 

This formula generalizes the one given by Eq. (92) extending its validity to N w 
(wpi anc k/TO„) 3 , leading to a baryonic core radius R c w 10km. We also recall that a 
more detailed analysis of nuclear matter in bulk in neutron stars (see e.g. Ref. 140, 
141) shows that at mass densities larger than the "melting" density of 

p c = 4.34 • 10 13 5 /cm 3 , (97) 

all nuclei disappear. In the description of nuclear matter in bulk we have to consider 
then the three Fermi degenerate gases of neutrons, protons and electrons. In turn 
this naturally leads to considering the role of strong and weak interactions among 
the nucleons. In the nucleus, the role of the strong and weak interaction, with a 
short range of one Fermi, is to bind the nucleons, with a binding energy of 8 MeV, 
in order to balance the Coulomb repulsion of the protons. In the neutron star case 
we have seen that the neutron confinement is due to gravity. We still assume that an 
essential role of the strong interactions is to balance the effective Coulomb repulsion 
due to the protons, partly screened by the electron distribution inside the neutron 
star core. We shall verify, for self-consistency, the validity of this assumption for the 
final equilibrium solution we will obtain. 

We now turn to the essential weak interaction role in establishing the relative 
balance between neutrons, protons and electrons via the direct and inverse /3-decay 

p + e — >n + v e , (98) 
n — >p + e + D e . (99) 

Since neutrinos escape from the star and the Fermi energy of the electrons is zero, 
as we will show below, the only non-vanishing terms in the equilibrium condition 
given by the weak interactions are 

[{P^cf + Mlc^ - M„c 2 = [(P p F cf + - M p c 2 + \e\V? oul , (100) 

where P% and Pp are respectively, the neutron and proton Fermi momenta, and 
Koui 1S t ne Coulomb potential of the protons. At this point, having fixed all these 
physical constraints, the main task is to find the electron distributions satisfying 
not only the Dirac-Fermi statistics but also the electrostatic Maxwell equations. 
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The condition of equilibrium for the Fermi degenerate electrons implies a zero value 
for the Fermi energy 

[(Pf c) 2 + m 2 c 4 ] 1/2 - mc 2 + eV coul (r) = 0, (101) 

where Pf is the electron Fermi momentum and V^ ou i(r) is the Coulomb potential. 

6.3. Relativistic Thomas-Fermi Equation for Nuclear Matter in 
Bulk 

In line with the procedure already followed for heavy atoms 132 ' 136 we adopt here 
the relativistic Thomas-Fermi Equation 

-i 3/2 ' 



1 cP X (x) J 1 



(102) 



where a = e 2 /(he), 9(x — x c ) represents the normalized proton density distribution, 
the variables x and \ are related to the radial coordinate and the electron Coulomb 
potential V con \ by 

—k{%T< ""0-^. (103, 
and the constants x c (r — R c ) and (3 are respectively 

3iV p \ 1/3 n mcR c ( 4tt ^ 1/3 



4* j < Ps -r\w P ) ■ (104) 

The solution has the boundary conditions 

X(0)=0; x(oo) = 0, (105) 

with the continuity of the function x an d its first derivative \' a * the boundary of 
the core R c . The crucial point is the determination of the eigenvalue of the first 
derivative at the center 

X'(0) = const., (106) 

which has to be determined by satisfying the above boundary conditions (105) and 
constraints given by Eq. (100) and Eq. (95). 

The difficulty of the integration of the Thomas-Fermi equations is certainly one 
of the most celebrated chapters in theoretical physics and mathematical physics, still 
challenging a proof of the existence and uniqueness of the solution and strenuously 
avoiding the occurrence of exact analytic solutions. We recall after the original 
papers of Ref. 142,143, the works of Ref. 144-147 all the way to the many hundred 
papers reviewed in the classical articles of Ref. 131,148,149. The situation here is 
more difficult since we are working on the special relativistic generalization of the 
Thomas-Fermi equation. We must therefore proceed by numerical integration in 
this case as well. The difficulty of this numerical task is further enhanced by a 
consistency check in order to satisfy all the various constraints. 
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We start the computations by assuming a total number of protons and a value 
of the core radius R c . We integrate the Thomas- Fermi equation and determine the 
number of neutrons from the Eq. (100). We iterate the procedure until a value of A 
is reached consistent with our choice of the core radius. The paramount difficulty of 
the problem is the numerical determination of the eigenvalue in Eq. (106) which al- 
ready for A w 10 4 had presented remarkable numerical difficulties. 136 In the present 
context we have been faced for a few months with an apparently insurmountable nu- 
merical task: the determination of the eigenvalue seemed to necessitate a significant 
number of decimal places in the first derivative (106) comparable to the number of 
the electrons in the problem! We shall discuss elsewhere the way we overcame this 
difficulty by splitting the problem on the basis of the physical interpretation of the 
solution. 150 The solution is given in Fig. (14) and Fig. (15). 



Ar (cm) 

10" 9 10" 8 10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 10 1 




r/R c 

Fig. 14. The solution \ °f the relativistic Thomas-Fermi equation for A = 10 57 and core radius 
R c = 10km is plotted as a function of radial coordinate. The left red line corresponds to the 
internal solution and it is plotted as a function of radial coordinate in units of R c in logarithmic 
scale. The right blue line corresponds to the solution external to the core and it is plotted as 
function of the distance Ar from the surface in a logarithmic scale in centimeters. 



A relevant quantity for exploring the physical significance of the solution is given 
by the number of electrons within a given radius r 

N e (r) = [ Air(r'fn e {r')dr' . (107) 
Jo 

This allows the determination of the distribution of the electrons inside and outside 
the core for selected values of the A parameter, and follows the progressive penetra- 
tion of the electrons in the core as A increases [ see Fig. (16)]. Then we can evaluate 
the net charge inside the core 

N nct = N p - N e (R c ) <N p , (108) 

generalizing the results in Ref. 132,136 and consequently determine the electric field 
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Fig. 15. The same as Fig. (14): enlargement around the core radius R c showing explicitly the 
continuity of function \ an d its derivative \' from the internal to the external solution. 

at the core surface, as well as inside and outside the core [see Fig. (17)] and evaluate 
as well the Fermi degenerate electron distribution outside the core [see Fig. (18)]. 

It is interesting to explore the solution of the problem under the same conditions 
and constraints imposed by the fundamental interactions and the quantum statistics 
and imposing the corresponding Eq. (95) instead of Eq. (90). Indeed a solution exists 
and is much simpler 

n n (x) = n p (x) = n e (x) — 0, x = 0. (109) 



6.4. The energetic stability of the solution 

Before drawing our conclusions we should check the theoretical consistency of the 
solution. We obtain an overall neutral configuration for the nuclear matter in bulk, 
with a positively charged baryonic core with 

\2/ \ 2 / -, \ 2 

m \ I e \ / 1 



*— 092 1 J I -I • (110) 



and an electric field at the baryonic core surface (see Fig. (17) ) 

^ = 0.92. (Ill) 
The corresponding Coulomb repulsive energy per nucleon is given by 

U£3 = mc 2 ^ 1.78 .10- 6 (MeV), (112) 

2a \ inn J 

well below the nucleon binding energy per nucleon. It is also important to verify 
that this charge core is gravitationally stable. We have in fact 

^_ = a -i/a (HL Y « i.5 6 . io-4. (us) 
VGM \m^J 
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r/R c 

Fig. 16. The electron number (107) in units of the total proton number N p is given as function of 
radial distance in units of the core radius R c , for selected values of A, again in a logarithmic scale. 
It is clear that by increasing the value of A, the penetration of electrons inside the core increases. 
The detail shown in Fig. (17) and Fig. (18) demonstrates how for N ~ (mpi an ck/ m n) 3 a relatively 
small tail of electrons outside the core exists and generates on the baryonic core surface an electric 
field close to the critical value. A significant electron density outside the core is found. 



Distance outside the core surface (units of X c ) 

1 10 100 

1 , , , 




1000 100 10 1 

Depth inside the core surface (units of X c ) 

Fig. 17. The electric field in units of the critical field E c is plotted around the core radius R c . 
The left (right) diagram in the red (blue) refers the region just inside (outside) the core radius 
plotted logarithmically. By increasing the density of the star the field approaches the critical field. 
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10 




Ar (cm) 

Fig. 18. The density of electrons for A = 10 57 in the region outside the core; both scales are 
logarithmic. 

The electric field of the baryonic core is screened to infinity by an electron distri- 
bution given in Fig. (18). 

As has been the case previously, any new solution for a Thomas- Fermi system has 
relevance and finds its justification in the domain of theoretical and mathematical 
physics. We expect that as in the case of other solutions that have appeared in the 
literature of the relativistic Thomas-Fermi equations, this new one presented here 
will find important applications in physics and astrophysics. There are a variety of 
new effects that such a generalized approach naturally leads to: (1) the energetics 
of the global neutrality solution is greatly different from the one obtained from 
the condition of local neutrality; (2) the formation process for a neutron star can 
also have specific new signatures, due to reaching a more tightly bound system; 
(3) we expect important consequences on the initial conditions in the physics of 
gravitational collapse of the baryonic core as soon as the protons and neutrons 
become relativistic and the critical mass for gravitational collapse to a black hole is 
reached. The consequent collapse to a black hole will have very different energetics 
properties, since the initial conditions will imply the existence of a critical electric 
field. Such a field will naturally lead to very strong processes of pair creation during 
the following phases of gravitational collapse. This research is ongoing. 

We now turn to the interpretation of the GRB data within the above theoretical 
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framework and recall some basic interpretational paradigms that we have introduced 
in order to reach a systematic understanding of these sources. 

7. The first paradigm: The Relative Space-Time Transformation 
(RSTT) paradigm 

The ongoing dialogue between our work and that of others who model GRBs still 
rests on some elementary considerations presented by Einstein in his classic arti- 
cle of 1905. 151 These considerations are quite general and even precede Einstein's 
derivation of the Lorentz transformations from first principles. We recall here Ein- 
stein's words: "We might, of course, content ourselves with time values determined 
by an observer stationed together with the watch at the origin of the coordinates, 
and coordinating the corresponding positions of the hands with light signals, given 
out by every event to be timed, and reaching him through empty space. But this 
coordination has the disadvantage that it is not independent of the standpoint of 
the observer with the watch or clock, as we know from experience. 

Einsteins message is simply illustrated in Fig. 19. If we consider in an inertial 
frame a source (solid line) moving with high speed and emitting light signals (dashed 
lines) along the direction of its motion, a far away observer will measure a delay 
At a between the arrival time of two signals respectively emitted at the origin and 
after a time interval At in the laboratory frame, which in our case is the frame 
where the black hole is at rest. The real velocity of the source is given by 

v 

and the apparent velocity is given by: 

As pointed out by Einstein, the adoption of coordinating light signals simply by 
their arrival time as in Eq. (115), without an adequate definition of synchroniza- 
tion, is incorrect and leads to insurmountable difficulties as well as to apparently 
"superluminal velocities as soon as motions close to the speed of light are considered. 

The use of At a as a time coordinate, often tacitly adopted by astronomers, 
should be done, if at all, cautiously. The relation between At a and the correct time 
parameterization in the laboratory frame has to be taken into account 

Ar 1 f to+At 

At a = At = At - - / v (t') dt' . (116) 

C C J to 

In other words, the relation between the arrival time and the laboratory time cannot 
be done without a knowledge of the speed along the entire world line of the source. 
In the case of GRBs, such a world line starts at the moment of gravitational collapse. 
It is of course clear that the parameterization in the laboratory frame has to take 
into account the cosmological redshift z of the source. We then have at the detector 

At d a = (l + z)At a . (117) 



Ar 
At 



Ar 
At? 



(114) 
(115) 
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In the current GRB literature, Eq. (116) has been systematically neglected by 
addressing only the afterglow description and neglecting the previous history of the 
source. Often the integral equation has been approximated by a clearly incorrect 
instantaneous value 



&t a ^— 2 . (118) 



At 

The approach has been adopted to consider the afterglow part of the GRB phe- 
nomenon separately without knowledge of the entire equation of motion of the 
source. 

This point of view has reached its most extreme expression in the work reviewed 
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by Ref. 153,154, where the so-called "prompt radiation" , lasting on the order of 10 2 
s, is considered as a burst emitted by the prolonged activity of an "inner engine. 
In these models, generally referred to as the "internal shock model, the emission 
of the afterglow is assumed to follow the "prompt radiation" phase. 155-159 As we 
outline in the following sections, such an extreme point of view originates from the 
inability to obtain the time scale of the "prompt radiation" from a burst structure. 
These authors consequently appeal to the existence of an "ad hoc" inner engine in 
the GRB source to solve this problem. 

We show in the following sections how this difficulty has been overcome in our 
approach by interpreting the "prompt radiation" as an integral part of the afterglow 
and not as a burst. This explanation can be reached only through a relativistically 
correct theoretical description of the entire afterglow (see next sections). Within 
the framework of special relativity we show that it is not possible to describe a 
GRB phenomenon by disregarding the knowledge of the entire past world line of the 
source. We show that at 10 2 seconds the emission occurs from a region of dimensions 
of approximately 10 16 cm, well within the region of activity of the afterglow. This 
point was not appreciated in the current literature due to the neglect of the apparent 
supcrluminal effects implied by the use of the "pathological" parametrization of the 
GRB phenomenon by the arrival time of light signals. 

We now turn to the first paradigm, the relative space-time transformation 
(RSTT) paradigm, 126 which emphasizes the importance of a global analysis of the 
GRB phenomenon encompassing both the optically thick and the afterglow phases. 
Since all the data are received in terms of the detector arrival time, it is essential 
to know the equations of motion of all relativistic phases of the GRB sources with 
7 > 1 in order to reconstruct the corresponding time coordinate in the laboratory 
frame, see Eq. (116). Contrary to other phenomena in nonrclativistic physics or 
astrophysics, where every phase can be examined separately from the others, in the 
case of GRBs all the phases are inter-related by their signals received in the arrival 
time if. In order to describe the physics of the source, there is the need to derive 
the laboratory time tas a function of the arrival time t~ along the entire past world 
line of the source using Eq. (117). 

An additional difference, also linked to special relativity, between our treatment 
and others in the current literature relates to the assumption of the existence of 
scaling laws in the afterglow phase: the power law dependence of the Lorentz gamma 
factor on the radial coordinate is usually systematically assumed. From the proper 
use of the relativistic transformations and by the direct numerical and analytic 
integration of the special relativistic equations of motion we demonstrate (see next 
sections) that no simple power-law relation can be derived for the equations of 
motion of the system. This situation is not new for workers in relativistic theories: 
scaling laws exist in the extreme ultrarelativistic regimes and in the Newtonian ones 
but not in the intermediate fully relativistic regimes (see e.g. Ref. 96). 
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8. The second paradigm: The Interpretation of the Burst 
Structure (IBS) paradigm 

We turn now to the second paradigm, which is more complex since it deals with 
all the different phases of the GRB phenomenon. We first address the dynamical 
phases following the dyadosphere formation. 

After the vacuum polarization process around a black hole, one of the topics of 
the greatest scientific interest is the analysis of the dynamics of the electron-positron 
plasma formed in the dyadosphere. This issue was addressed by us in a collaboration 
with Jim Wilson at Livermore. The numerical simulations of this problem were 
developed at Livermore, while the semi-analytic approach was developed in Rome 
(see Ruffini et al. 106 ' 107 and next sections). The corresponding treatment in the 
framework of the Cavallo, Rees et al. analysis was performed by Piran et al. 160 also 
using a numerical approach, by Bisnovatyi-Kogan & Murzina 161 using an analytic 
approach and by Meszaros et al. 162 using a numerical and semi-analytic approach. 

Although some similarities exist between these treatments, they are significantly 
different in the theoretical details and in the final results (see Ref. 163 and next 
sections). Since the final result of the GRB model is extremely sensitive to any 
departure from the correct treatment, it is indeed very important to detect at every 
step the appearance of possible fatal errors. 

8.1. The optically thick phase of the fireshell 

A conclusion common to all these treatments is that the electron-positron plasma is 
initially optically thick and expands till transparency reaching very high values of 
the Lorentz gamma factor. A second point, which is also common, is the discovery of 
a clearly new feature: the plasma shell expands but the Lorentz contraction is such 
that its width in the laboratory frame appears to be constant. This self acceleration 
of the thin shell is the distinguishing factor of GRBs, conceptually very different 
from the physics of a fireball developed by the inner pressure of an atomic bomb 
explosion in the Earth's atmosphere. In the case of GRBs the region interior to the 
shell is inert and with pressure totally negligible: the entire dynamics occurs on the 
shell itself. For this reason, we refer in the following to the self accelerating shell as 
the "fireshell." 

There is a major difference between our approach and those of Piran, Meszaros 
and Rees in that we assume the dyadosphere to be initially filled only with an 
electron-positron plasma. Such a plasma expands in substantial agreement with the 
results presented in the work of Ref. 161. In our model the fireshell of electron- 
positron pairs and photons (PEM pulse) 106 evolves and encounters the remnant of 
the star progenitor of the newly formed black hole. The fireshell is then loaded with 
baryons. A new fireshell is formed of electron-positron-photons and baryons (PEMB 
pulse) 107 which expands all the way until transparency is reached. At transparency 
the emitted photons give origin to what we define as the Proper-GRB (P-GRB, see 
Ref. 127 and Fig. 20). 
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Fig. 20. Above: The optically thick phase of the fireshcll evolution are qualitatively represented 
in this diagram. There are clearly recognizable 1) the PEM pulse phase, 2) the impact on the 
baryonic remnant, 3) the PEMB pulse phase and the final approach to transparency with the 
emission of the P-GRB. Details in Ref. 152. Below: The P-GRB emitted at the transparency point 
at a time of arrival which has been computed following the prescriptions of Eq. (116). Details 
in Ref. 127,152. 



In our approach, the baryon loading is measured by a dimensionless quantity 

5 = ^, (119) 



-E^dya 
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which gives direct information about the mass Mb = Nsnip of the remnant, where 
m p is the proton mass. The corresponding treatment done by Piran and collabora- 
tors 160,164 and by Ref. 162 differs in one important respect: the baryonic loading is 
assumed to occur from the beginning of the electron-positron pair formation and 
no relation to the mass of the remnant of the collapsed progenitor star is attributed 
to it. 

A further difference also exists between our description of the rate equation for 
the electron-positron pairs and the ones by those authors. While our results are 
comparable with the ones obtained by Piran under the same initial conditions, the 
set of approximations adopted by Ref. 162 appears to be too radical and leads 
to very different results violating energy and momentum conservation (see next 
sections and Ref. 163). 

From our analysis 107 it also becomes clear that such an expanding dynamical 
evolution can only occur for values of B < 10~ 2 (see Fig. 21). This prediction, as 
we will show shortly in the many GRB sources considered, is very satisfactorily 
confirmed by observations and is indeed essential in order to reach the high values 
of the Lorentz gamma factor observed in GRBs. 

From the value of the B parameter, related to the mass of the remnant, it 
therefore follows that the collapse to a black hole leading to a GRB is drastically 
different from the collapse to a neutron star. While in the case of a neutron star 
collapse a very large amount of matter is expelled, in many instances well above the 
mass of the neutron star itself, in the case of black holes leading to a GRB only a 
very small fraction of the initial mass (~ 10~ 2 or less) is expelled. The collapse to 
a black hole giving rise to a GRB appears to be much smoother than any collapse 
process considered until today: almost 99.9% of the star has to collapse at once to 
form the black hole! 

We summarize in Fig. 20 the optically thick phase of the fireshell evolution: we 
start from a given dyadosphere of energy Ed ya ; the fireshell self-accelerates outward; 
an abrupt decrease in the value of the Lorentz gamma factor occurs due to the 
cngulfment of the baryonic loading followed by a further self-acceleration until the 
fireshell becomes transparent. 

The photon emission at this transparency point is the P-GRB. An accelerated 
beam of baryons with an initial Lorentz gamma factor j starts to interact with the 
interstellar medium at typical distances from the black hole of r a ~ 10 14 cm and 
at a photon arrival time at the detector on the Earths surface of if ~ 0.1 s. These 
values determine the initial conditions of the afterglow. 

8.2. Hydrodynamics of the pair plasma 

We give a systematic derivation of the main equations, present a critical review 
of existing models for isotropic relativistic fireballs, and compare and contrast 
these models, following Ref. 163. In the next section, following Ref. 165 we de- 
rive basic equations and describe the approximations involved. Then we present the 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



48 



10 1£ 
10 1E 
10 17 I- 
10 16 
10 15 - 
10 14 h 

10 13 

10 12 

10 11 
10 io 

10 9 
10 8 I- 
10 7 
10 6 - 
10 5 - 
10 4 
1000 
100 
10 

1 








200 400 600 

Radius/R dyadosphere 



800 



Fig. 21. A sequence of snapshots of the coordinate baryon energy density is shown from the one 
dimensional hydrodynamic calculations of the Livermore code. The radial coordinate is given in 
units of dyadosphere radii (r^). At r ~ lOOr^ there is located a baryonic matter shell corre- 
sponding to a baryon loading B = 1.3 X 10~ 2 . For this baryon shell mass we see a significant 
departure from the constant thickness solution for the fireshcll dynamics and a clear instability 
occurs. Details in Ref. 107. As we will see, this result, peculiar to our treatment, will play a major 
role in the theoretical interpretation of GRBs. 



model 106,107 which differs from other models in the literature in that it describes the 
dynamics of the fireshell taking into account the rate equations for electron-positron 
pairs. Then we compare and contrast the above mentioned models. 
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8.2.1. Local, global and average conservation laws 

Particle number. The first relevant equation represents continuity of relativistic 
flux and reads c 

where n is the number density of relativistic fluid and [7 M is its velocity field. Defining 
the particle number by 

N= f V^gnU°dV, (121) 
Jv 

we see that 

^ = - f y/^nWdV = -<£ J~gnXJ l dS u (122) 
dt Jv j s 

where we have used the Ostrogradsky-Gauss theorem. Thus if particles do not cross 
the surface £ bounding the volume V considered under consideration, the total 
number of particles is constant during the system evolution. 

Now assume spherical symmetry 11 , which is usually done for fireball descriptions. 
With spherical spatial coordinates x l = {r, cp} the interval is 

ds 2 = -dt 2 + dr 2 + r 2 d{) 2 + r 2 sin 2 ddip 2 . (123) 

Assuming the absence of fluxes through the boundary S we rewrite (120) 



^ + ^|(--W^T)=0. ,124, 

Integrating this equation over the volume from some initial r^(t) to a radius r e (t) 
which we assume to be comoving with the fluid 

^f=fe*), d ^=^t), (125) 



c Greek indices denote four-dimensional components and run from to 3 while Latin indices run 
from 1 to 3. General relativistic effects are neglected, which is a good approximation, but we use 
the general definition of energy-momentum conservation to take into account the most general 
coordinate system. 

d The only nonvanishing components of the energy-momentum tensor are T 00 , T 01 , T 10 , T 11 , T 22 
, T 33 . The factor \/—g = r 2 sin 1? in all expressions above becomes simply a volume measure and 
the differentials are dV = drdi3d(p, dS = di3d(p, so the differential laboratory volume can be written 
as dV = yJ—gdV = r 2 sin tfdrdtfdip. 
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and ignoring a factor An we have 



d (n-f ) 
dt 



r z dr + / — r^nV7 2 - l) dr = (126) 



— I (n 7 ) r 2 dr - -^-n(r e , t)j(r e , t)r e + ~^ n (n, t)j(ri,t)r 2 + 
+r 2 e n(r e ,t)^-y 2 (r e ,t) - 1 - r 2 n(n,t)^ 2 (ri,t) - 1 = 

= dt J ^ r2dr = °' 

Since we deal with arbitrary comoving boundaries, this means that the number 
of particles in each shell between the boundaries is conserved as well as the total 
number of particles integrated over all shells, in other words, 

R(t) 

„2, 



N = 4ir J n~/r 2 dr = const, (127) 
o 

where R(t) is the external radius of the fircshell. 

Following Ref. 160 one can transform (124) from the variables (t, r) to the new 
variables (s = t — r, r) and then show that 



r 2 dr 



(^»v^T) - -I (^-) ■ 

Now assume the expansion velocity is ultrarelativistic, 

7>1- (129) 

In this approximation, therefore, 

dN = Airnjr 2 dr w const. (130) 
Relations (130) and (127) then imply 

47T J (rryr 2 ) dr = 4tt [n(r, t)f(r, t)r 2 ] J dr = 4ir (wyr 2 ) A w const, (131) 

where the first argument of the functions n(r,t) and j(r,t) is restricted to the 
interval r, < r < r e and 

A = r e — Ti i=a const. (132) 

This means that the fluid shell does not broaden, but rather has a constant thickness. 
This fact proves the constant thickness approximation, adopted in 106,107. 
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The volume element measured by the observer outside the hrcshcll (to be referred 
to as the lab frame in what follows), for which it appears to be moving with velocity 
(3 is just 

dV = ATrr 2 dr, (133) 

while the volume element comoving with the fireshell, for which the fluid is at rest, 
is 

dV = 47r 7 r 2 dr, (134) 
with the conversion of the volumes 

dV = jdV. (135) 
Then the average value of the Lorentz factor is defined as follows 

47T / jr 2 dr v 

<7> = ^ S = y • (136) 

47T / r 2 dr 

Now we can formulate the conservation law for the average value of the number 
density in the lab frame 



Kab= V = * !i:;7 ' 



47T J njr 2 dr 

r 

4tt J r 2 dr 



Assuming r » A we then obtain 



(">lab - 4 3A A = * )7(r ' t] " (138) 



Therefore, the average number density in the lab frame scales like r 2 . 

At the same time, recalling the expression for the divergence of the four-velocity 

^ - if, im 

where r is the proper time, and remembering that U^-^ = from (120) we get 



, t t 1 1 \ tt,, tt„ dn n dV 

(nU^ = + nEP ; „ = — + -— = 0, 

dlnn + dlnV = 0. (140) 

This means that the number of particles is conserved along the flow lines of the 
fluid. The solution of this equation provides the definition for the comoving average 
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number density 



47r / rvyr dr 



47r y 7r 2 dr 

Clearly, the condition (141) gives a link between the description of the fireshell 
evolution in terms of local functions entering (124) on one side of the equation, 
and global quantities (136) and (138) on the other. The presence of the global 
conservation (127) in both these cases ensures equivalence of the local (124) and the 
average (137) descriptions for the fireshell, unless its detailed structure is considered. 

Energy-momentum conservation. The basis of the description of a relativistic 
fireshell is the energy-momentum principle. It allows one to obtain the relativistic 
hydrodynamic equations, or equations of motion for the fireshell, the energy and 
momentum conservation equations which are used extensively to describe inter- 
action of relativistic baryons of the fireshell with the interstellar matter, and the 
boundary conditions which are used to understand shock wave propagation in the 
decelerating baryons and in the outer medium. 

Consider energy-momentum conservation in the most general form: 

(T„ ") iV - V) + V~gW x = 0, (142) 

where T^ x are the Christoffel symbols and g is the determinant of the metric tensor. 
Integrating over the entire three-dimensional volume we obtain 

/ T M v , v dV = 0. (143) 
Jv 

Integrating over the entire four-dimensional volume and applying the divergence 
theorem we get 166 

/ / T„ v . v dVdt = I T„ v KdV = 0, (144) 
Jt Jv Jv 

where \ a are covariant components of the outer normal to the three-dimensional 
hypersurface (volume V) enclosing the spacetime region. 

Now suppose that there is a discontinuity on the fluid flow. Taking the volume to 
be a spherical shell and choosing the coordinate system in which the discontinuity 
is at rest so that in (144) for normal vectors to the discontinuity hypersurface X a , 
we have 



A Q A a — 1, Aq — 0. 



(145) 
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Let the radius of the shell R s be very large and the shell thickness A very small. 
With R s — > oo and A — » from (144) wc get 

[T m ] = 0, (146) 

where the brackets mean that the quantity inside is the same on both sides of the dis- 
continuity surface. This equation together with continuity condition for particle den- 
sity flux [nU l ] = was used by Ref. 166 to obtain the relativistic Rankine-Hugoniot 
equations. These equations govern shock wave dynamics which are supposed to ap- 
pear during the collision of the baryonic material left from the fireshell with the 
interstellar medium. 167 The origin of the afterglow could be connected 168-170 to the 
conversion of kinetic energy into radiative energy in these shocks. 

Consider now the energy-momentum tensor of the perfect fluid in the lab frame 
(where the fluid was initially at rest) 

=pg^ +luU^U v , (147) 

where u> = p + p is the proper enthalpy, p is the proper pressure and p is the proper 
internal energy density. 

Rewrite (142) in the spherically symmetric case 

2 | ! + ^|:(^')-;(T' 3 » + r 3 =)=0, (149) 
arriving at the equations of motion of a relativistic fireshell 106,160,167 



9(7 w) dp Id 
dt dt r 2 dr 



(r 2 7 2 /3w) = 0, (150) 



where the four-velocity and the relativistic Lorentz factor are defined as follows 6 

E7" = (7,70,0,0), 7 ^(l-/3 2 )- 1 / 2 , (152) 

and (3 is the radial velocity. 

The total momentum of the spherically symmetric expanding shell vanishes. 
However, from the local conservation equations (150) one finds that the radial com- 
ponent of the four-momentum vector does not vanish. In analogy with the continuity 
equation (124) we integrate the first equation in (150) over the volume starting from 
some internal radius j-j(t) up to some external radius r e (t), and ignoring a factor 



c Throughout this chapter we set the speed of light equal to 1. 
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4tt we obtain 

-| J " j 2 Lur 2 dr + r?7 2 (riMrO/3(ri) - r 2 7 2 (r e )w(r e )/?(r e )- (153) 

— / pr 2 dr + rfp(ri) - r 2 p(r e ) + r 2 7 2 (r e )w(r e )/3(r e )- (154) 
r^ 2 ^)^^)/?^) = 0. 

If the boundaries rj(i) and r e (t) are comoving with the fluid we have 

■ y (j 2 oj - p) r 2 dr = r 2 e p(r e ) - r 2 p(ri). (155) 
Further, if one assumes (129), one gets the following result 

rR(t) 

E = 4tt / -f 2 ujr 2 dr = const. (156) 
Jo 

The differential conservation law follows from the same arguments which lead 
to (130), so we also have 

dE = 4ir~f 2 u>r 2 dr w const. (157) 
Analogously to (137) we introduce the average energy density in the lab frame 



d_ 
~dt 



At: j (7 2 w)r 2 dr 

'I' 



(PK b - f = . (158) 

47T / r 2 dr 



Taking the polytropic equation of state with the thermal index 

r = i + jj, (159) 

and requiring also r » A and (129) we find from (158) 

(Aab^PM7 2 Mcxr- 2 . (160) 
The radial momentum equation follows from (151) 
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dt 



9(7 ^ r 2 dr +TM V 2 l) «] ^ + f ^ 



^ ' 7 2 /?wr 2 dr + rf 7 2 (r 4 V(r,)/3 2 (r 2 ) - r 2 el 2 (r e )cu(r e )f3 2 (r e ) + (161) 



+r 2 frVe) - 1) uj(r e ) - r 2 ( 7 2 (r t ) - l) co(r 4 ) + ^ " ||r 2 dr : 
9 ( 7 2 /3 W ) r 2 dr + r 6 ^r 2 dr - 0. 

J n J ri 



■ ( 7 2 /3w) r 2 dr = 2 / prdr + r 2 p(r t ) - r 2 e p(r e ). (162) 



<9i ./,,_. .L- dr 

This leads to 

d 
~dt 

For the radial momentum we have 

dp. , d r R(t) r R{t) 

^J2L = fL / 47r ( 7 2^) r 2 dr = 8jT / prdr (163) 

The left hand side of this equation is the time derivative of the radial momentum, 
i.e., the radial "force". The right hand side is the integral of the pressure over all 
the shells, so it is clear that unless the pressure in the fireshcll is zero, it experiences 
self-acceleration due to internal pressure. 

Entropy conservation. Yet another relevant equation is entropy conservation 
which may be obtained from (142) by projection along the flow line 

U» (T„ ") .„ = (U^T, ") ;„ - T„ " (Un , v = (164) 
= - ( P Un ;M - uU v (U^. v ) - p£/% = 0. 

The second term on the last line vanishes since U^U^ = — 1, so we have another 
conservation equation 

-U" (7jt ") ,„ - ( P Un ,„ + P Ur^ = 0. (165) 

This conservation law corresponds to another conserved quantity, the entropy. 
In fact, Eq. (165) can be rewritten as 

(pUn iM + pUr^ = (ujUn ,„ - = 0. (166) 

Now using the continuity equation (120) and the identity uiU^ = nU^ (^) we find 



0. (167) 



But the functions inside the brackets are scalars, and therefore covariant deriva- 
tives can be replaced by ordinary derivatives. Then we recall the second law of 
thermodynamics 171 

d(-) =Td(-) + -dp, (168) 
\nJ \nJ n 
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and finally obtain 



nTU" (-) = 0, (169) 



which can be rewritten using (120) as 

(aUn ; , = 0. (170) 

This is the continuity equation for the entropy. Since it has exactly the same form 
as (120), all the conservation equations such as (130) and (127) hold for the entropy 
as well, 

da = 47T (0-7) r 2 dr w const, (171) 

R(t) 

S = J da = const. (172) 


Assuming (159) we find from (165) and (139) the following result 

+ = 

d\np + Td\nV = 0, (173) 
(P)com U r = const. 

Finally, due to the similarity of Eqs. (130) and (171), the average entropy can 
be defined in the same manner as (137). 

Analogy with a Friedmann Universe. It is easy to show that the conservation 
equations (130), (157) and (171) imply an analogy between the fireball and the 
Friedmann Universe, noticed first by Shemi and Piran. 160,164 In fact, this analogy 
is valid for a polytropic equation of state (159) and the ultrarelativistic expansion 
condition (129). First, using (159) and the integral form of (168) we obtain 

a = r 

which leads to 

poca r , (174) 
we rewrite the above mentioned conservation equations using (159) 

njr 2 — const, 
pTj r 2 = const, (175) 
/?7 2 r 2 = const. 
From these equations we then easily find 

2(r-i) 

7 oc r 2 - r , 

notr", (176) 
p oc r 2 - r . 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



57 

Taking the ultrarelativistic equation of state with r = 4/3 we immediately obtain 

7 cx r, 



-3 



-4 

p cx r , 



(177) 



as opposed to the nonrelativistic equation of state with T = 1 with different scaling 

7 = const, 

ncxr~ 2 , (178) 
p cx r~ 2 . 

Actually, scaling laws (177) take place for the homogeneous isotropic radiation- 
dominated Universe. 160,164 This fact allowed the authors of Ref. 160 to speak about 
the frozen-pulse profile for 7 » 1 where number, energy and entropy density is 
conserved within each differential shell with thickness dr, although the radial dis- 
tribution of matter and energy can be inhomogeneous. 

Although for observer inside the radiation-dominated fircshcll it looks indistin- 
guishable from a portion of radiation-dominated Universe, for the observer outside 
it looks drastically different. In fact, the validity of differential conservation laws 
(130), (157) and (171) together with integral ones (127), (156) and (172) implies 
constant thickness approximation assumed in 106,107. 

Clearly, if the condition (129) is satisfied, also 

A < R(t) (179) 

is valid. Given the scalings (177) we then find 

f R(t) f R(t) 
V = 4rr / jr 2 dr ~ 4tt / Sr 3 dr ~ 4iri? 3 , (180) 

JR(t)-A JR(t)-A 

where we put 7 = Sr, S is a constant. At the same time, 



V = 4tt 



/ r 2 dr = -Ji? 3 . (181) 
Jo 3 



Equality of (180) and (181) up to a numerical factor suggests that the initially 
homogeneous energy and particle number distribution looks highly compressed in 
the lab frame expanding with ultrarelativistic velocity and compression factor 7. 



8.2.2. Self- acceleration of the fireshell 

For a fireshell which is initially optically thick the total energy is conserved. Assume 
that the fireshell consists of relativistic electrons, positrons and photons, and also 
that some admixture of a plasma in the form of photons and electrons is present 
such that the total charge is zero. While electrons are relativistic, the protons are 
not. The equation of state for pairs and electrons in such a case is given to a good 
approximation by that of an ultrarelativistic fluid: p e ±.-y = p e ± t7 /3. At the same 
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time for protons we have p p ~ 0. Therefore positrons and electrons together with 
photons can be considered to be one fluid with p r = p r /3 since they are strongly 
coupled since the medium is optically thick. Instead protons have low pressure and 
little internal energy compared to their rest mass energy. 
According to (156) we find 

R(t) R(t) R(t) 

J (7 2 w - p)r 2 dr = J ~/ 2 p p r 2 dr + ^ J ^ 2 p r r 2 dr. (182) 



These two terms are the rest mass energy of protons Mb and the energy of the 
ultrarelativistic fluid E correspondingly, so we arrive at a simple result, expressing 
the total energy of the expanding relativistic shell in the lab frame 

-y{E + M) = const, (183) 

which reads simply as E + M = const in the comoving frame taking into account 
the conversion of volumes (135). 

For homogeneous distributions of matter, energy density and pressure the inte- 
grals (126), (155) and (162) reduce to 

717 V = const, 
[l 2 (P + P)] V = const, (184) 

while in the comoving frame instead we would have 

nV = const, (185) 
pV = const, (186) 

which means the energy and number of particles do not change. 
From the above we have 



nU° om V = nV = const = nlI? ab V = njV, (187) 
T™ m V = P V = const = TZV = [ 7 2 (p + p)} V, (188) 

remembering that all quantities n, p,p are always defined as comoving ones. 
Energy conservation (182) for (129) implies 




,„ + r„v < 189 > 



Clearly all the equations given above can be written for the average values of 
the number and energy densities. 
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8.2.3. Quasi-analytic model of GRBs 

The first detailed models for the expansion of a relativistic fireball were suggested 
in the beginning of nineties. 160,162,164 Independent calculations were performed in 
Ref. 106 and 107. The main difference of these last two articles from the other models 
in the literature is that initially not photons but pairs are created by an overcritical 
electric field, and these pairs produce photons later. This plasma referred to as the 
pair-electro-magnetic (PEM) pulse expands initially into the vacuum surrounding 
the black hole reaching relativistic velocities very quickly. Then the collision with 
the baryonic remnant of the collapsed star takes place and the PEM pulse becomes 
a pair-electro-magnetic-baryonic (PEMB) pulse, see Ref. 152 for details. This dif- 
ference is not large, since it was shown that the final gamma factor does not depend 
on the distance from the baryonic remnant or on the parameters of the black hole. 
The only crucial parameters are again the initial energy Eq and baryonic admixture 
B. 

The model is based on the numerical integration of the relativistic energy- 
momentum conservation equations (150,151) together with the baryonic number 
conservation equation (120). However, the most important point distinct from all 
previous models is that the rate equation for electron-positron pairs is added to the 
model and integrated simultaneously in order to give a more detailed description of 
the transparency. This latter fact leads to quantitative differences in predictions of 
the model with respect to the simplified models in the literature. 

Here we concentrate on the simple quasi-analytical treatment presented in 
106,107, see also 152. The PEMB pulse is supposed to contain a finite number 
of shells each with a flat density profile. Their dynamics is governed by the set 
of equations derived in the previous subsections. We collect (141), (173) and (189) 
together (omitting brackets for brevity) 

(190) 
(191) 

\ vo / \ VQ J \70/ 

2- = >° +r ^ V0 , (192) 

where subscript "0" denotes initial values, and all quantities are assumed to be av- 
eraged over a finite distribution of shells with constant width and density profiles. 
All components such as photons, electrons, positrons and plasma ions give a con- 
tribution to the energy density and pressure. This set of equations is equivalent to 
(200) and (202) (see below). The next step is to take into account the rate equation 
for positrons and elections, accounting for non- instant transparency: 




(n e ±U^=^(n 2 e± (T)-n 2 e± ), 



(193) 
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or, integrating over the volume 

l Ne± =-N e± ^+^(N 2 ± (T)-N 2 ± ), (194) 

where a is the mean pair annihilation-creation cross section and is the thermal 
velocity of e ± -pairs. The coordinate number density of e ± -pairs in equilibrium is 
N e ± (T) = 7n e ± (T) and the coordinate number density of e ± -pairs is N e ± — 7n e ± . 
For T > m e c 2 we have n e ±(T) ~ n 7 (T), i.e., the number densities of pairs and 
photons are nearly equal. The pair number densities are given by appropriate Fermi 
integrals with zero chemical potential at the equilibrium temperature T. For an 
infinitesimal expansion of the coordinate volume from Vo to V in the coordinate 
time interval t — to one can discretize the last differential equation for numerical 
computations. 

The most important results of the analysis performed in 107 are the following: 

• the appropriate model for the geometry of an expanding fircshcll (PEM-pulse) 
is given by the constant width approximation (this conclusion is achieved by 
comparing results obtained using (150), (151) and the simplified treatment de- 
scribed above), 

• there is a bound on parameter B which comes from violation of the constant 
width approximation, B < 10~ 2 (rj > 10 2 ). 

In a previous subsection we proved the applicability of the constant thickness 
approximation for the fireshell. The second conclusion appears to be crucial, since 
it shows that there is a critical loading of baryons at which their presence pro- 
duce a turbulence in the outflow from the fireshell, and its motion becomes very 
complicated and the fireshell evolution does not lead in general to a GRB. 

Exactly because of this reason, the optically thick fircshcll never reaches a radius 
as large as r& = r r] 2 which is discussed in 162, see section (8.2.4), since to do this 
the baryonic fraction should exceed the critical value B c = 10~ 2 . For larger values 
of B c the theory reviewed here does not apply. This means in particular that all 
the conclusions of Ref. 162 obtained for r > r& are invalid. In fact, for B < B c the 
gamma factor even does not reach saturation. 

Notice that another way to obtain the constraint B < B c is to require the optical 
depth of the emitting region to be smaller than 1, leading to the requirement that 
the Lorentz factor be greater than 7 > 10 2 , see the introduction. At the same 
time, there is a simple relation between the Lorentz factor and the baryonic loading 
parameter B = 7 _1 in the region 10~ 2 < B < 10~ 4 , see Fig. 23, which leads to 
B < 10~ 2 . 

The fundamental result coming from this model are the diagrams shown in 
Fig. 22 and Fig. 23. The first one shows basically which portion of the initial energy 
is emitted in the form of gamma rays E 1 when the fireshell reaches the transparency 
condition r ~ 1 and how much energy gets converted into the kinetic energy Ek of 
the baryons left after pair annihilation and the photons escape. 
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B 



10 8 10 7 10 6 10 5 10 4 10 3 10 2 

Fig. 22. Relative energy release in the form of photons emitted at the transparency point (solid 
line) and the kinetic energy of the plasma (dashed line) of the baryons in terms of the initial energy 
of the fireball depending on parameter B obtained on the basis of quasi-analytic model. The thick 
line denotes the total energy of the system in terms of the initial energy. 

The second diagram gives the value of gamma factor at the moment when the 
systems reaches transparency. 

Energy conservation holds, namely 

E = E y + E k . (195) 

Clearly when the baryon abundance is low, most energy is emitted when the fireshell 
becomes transparent. It is remarkable that almost all initial energy is converted into 
the kinetic energy of the baryons already in the region of validity of the constant 
thickness approximation B < 10~ 2 , so the region 10~ 8 < B < 10 -2 is the most 
interesting from this point of view. 

8.2.4. Alternative models 

Shemi and Piran model. In this section we discuss the model proposed by 
Ref. 164. This quantitative model gives a rather good general picture of relativistic 
fireballs. Shemi and Piran found that the temperature at which the fireball becomes 
optically thin is determined to be 



T esc = mm(Tg,T p ), 



(196) 
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B 



Fig. 23. Relativistic gamma factor of the fireball when it reaches transparency depending on the 
value of parameter B. The dashed line gives the asymptotic value. 

where T g and T p is the temperature when it reaches transparency with respect to 
the gas (plasma) or the pairs: 



2 45 m p 1 1 

3 ~ Sir* m e a2g i TJKO 111 



(197) 
(198) 



T p ~ 0.032, 

where m p is the proton mass, go = a — ™, and the dimensionless temperature 
T and radius 1Z of the fireball are measured in units of m , cC and A e = — — cor- 
respondingly, and the subscript "0" denotes initial values. The temperature at the 
transparency point in the case when the plasma admixture is unimportant is nearly 
a constant for a range of parameters of interest and it almost equals 



T p = 15 keV. 



(199) 



Adiabatic expansion of the fireball implies 



E 

Eo 



T 
To 



no 



(200) 
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where E — m E S is a radiative energy. From the energy conservation (142), assuming 
the fluid to be pressureless and its energy density profile to be constant, we have in 
the coordinate frame 



/ 



T °dV = jE tot = const. (201) 



Assuming at initial moment 70 = 1 and remembering that E tot — E + Mc 2 we 
arrive at the following fundamental expression for the relativistic gamma factor 7 
at the transparency point: 



E + Mc 2 77+I 

^-^W-(^TT' (202) 

where M = ^. 

One can use this relation to get such important characteristics of the GRB as 
the observed temperature and observed energy. In fact, they can be expressed as 
follows: 



T obs = 7 T esc , (203) 

E obs =E^. (204) 

These results are shown in Fig. 24. In the limit of small 77 we have 7 = (1 + 77), 
while for very large 77 the value of the gamma factor at the transparency point is 
7 = 7~o/7~ esc , and it has a maximum at intermediate values of 77. We denote by a 
dashed thick line the limiting value of the 77 parameter r\ c = B" 1 . For 77 < r\ c the 
approximations used to construct the model do not hold. It is clear that because 
of the presence of bound r\ c the value 7 = 77 can be reached only asymptotically. In 
effect, the value r/ c cuts the region where saturation of the gamma factor happens 
before the moment when the fireball becomes transparent. 

It was found that for relatively large 77 > 10 5 the photons emitted when the 
fireball becomes transparent carry most of the initial energy. However, since the 
observed temperature in GRBs is smaller than the initial temperature of the fireball, 
one may assume that a large part of the initial energy is converted to kinetic energy 
of the plasma. 

Shemi, Piran and Narayan model. 160 presents a generalization of this model 
to arbitrary initial density profile of the fireball. These authors performed numeri- 
cal integrations of energy-momentum relativistic conservation equations (150), (151) 
and baryon number conservation equation (122). They were mainly interested in the 
evolution of the observed temperature, gamma factor and other quantities as the the 
radius increases. Their study results in a number of important conclusions, namely: 
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Log 




Logr) 



Fig. 24. The relativistic gamma factor (upper dashed line), the observed temperature (solid line), 
and the ratio of observed energy to the initial energy of the fireball (lower dashed line) as a function 
of baryonic loading parameter, see Ref. 164. The values of parameters are the same as in the cited 
paper. The thick dashed line denotes the limiting value of the baryonic loading. Its values when 
the gamma factor reaches a maximum and becomes constant are also shown. 

• The expanding fireball has two basic phases: a radiation dominated phase and 
a matter-dominated phase. In the former, the gamma factor grows linearly with 
the radius of the fireball: 7 oc r, while in the latter the gamma factor reaches 
the asymptotic value 7 ~ 77 + 1. 

• The numerical solutions are reproduced with good accuracy by the frozen-pulse 
approximation when the pulse width is given by initial radius of the fireball. 

The last conclusion is important, since the fireball becomes a fireshell, and the 
volume V can be calculated as 



V = 4irR 2 A, (205) 

where A ~ i?o is the width of the leading shell with constant energy density profile 
and R is the radius of the fireball. 

They also present the following scaling solution: 

/ \ 1/2 

R = Ro( —D 3 ) , (206) 



7 

D 7 4777 477' ^ ^ 
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where the subscript "0" denotes some initial time when 7 > 2, which can be inverted 
to give 7(.R). 

Meszaros, Laguna and Rees model. The next step in developing this model 
was made by Ref. 162. In order to reconcile the model with observations, these 
authors proposed a generalization to the anisotropic (jet) case. Nevertheless, their 
analytic results apply to the case of homogeneous isotropic fireballs and we will 
follow their analytical isotropic model in this section. 

Starting from the same point as Shemi and Piran, consider (200) and (202). 
The analytic part of the paper describes the geometry of the fireball, the gamma 
factor behavior and the final energy balance between radiation and kinetic energy 
Magnetic field effects are also considered, but we are not interested in this part here. 

Three basic regimes are found in Ref. 162 for the evolution of the fireball. In 
two first regimes there is a correspondence between the analysis in the paper and 
results of Ref. 160, so the constant thickness approximation holds. It is claimed in 
162 that when the radius of the fireball reaches very large values comparable to 
Rb = Rorj 2 , a noticeable departure from the constant width of the fireball occurs. 
However, it is important to note that the fireball becomes transparent much earlier 
and this effect never becomes important (see section 8.2.3). 

The crucial quantity presented in the paper is T m — the maximum possible bulk 
Lorentz factor achievable for a given initial radiation energy E deposited within a 
given initial radius R : 



r m ee Vm = (Tor/) 173 = (£oK77) 1/3 , (208) 

M <7t 



S ° = r^2' «=— > (209) 



where S is initial baryon (plasma) mass surface density. 

All subsequent calculations in the paper Ref. 162 involve this quantity. It is 
evident from (208) that a linear dependence between the gamma factor T and 
parameter r\ is assumed. However, this is certainly not true as can be seen from 
Fig. 23. We will come back to this point in the following section. 

Another important quantity is given in this paper, namely 




(210) 



This is just the asymptotic behavior of the gamma factor at Fig. 24 for very large 
r\. Using it, the authors calculate the value of the ij parameter above which the pair 
dominated regime occurs: 



(211) 
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This means that above r\ v the presence of baryons in the fireball is insufficient to 
keep the fireball opaque after pairs are annihilated and almost all initial energy 
deposited in the fireball is emitted immediately. 

The estimate of the final radiation to kinetic energy ratio 162 is incorrect, because 
kinetic and radiation energies do not sum up to the initial energy of the fireball thus 
violating energy conservation. This is illustrated in Fig. 25. The correct analytic 
diagram is instead presented in Fig. 22. 



E/E 




Fig. 25. The ratios of radiation and kinetic energy to the initial energy of the fireball predicted 
by the Meszaros, Laguna and Rees model. The thick line denotes the total energy of the system 
in terms of the initial energy. Energy conservation does not hold. 



Approximate results. All models for isotropic fireballs are based on the following 
points: 

(1) flat space-time, 

(2) relativistic energy-momentum principle, 

(3) baryonic number conservation. 

Although the model 106 ' 107 starts with Reissner-Nordstrom geometry, the numer- 
ical code is written for the case of flat space-time simply because curved space-time 
effects become insignificant soon after the fireshell reaches relativistic expansion ve- 
locities. The presence of the rate equation in the model 106,107 has a strong physical 
basis and its absence in the other treatments means incompleteness of their mod- 
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els. Indeed, the number density of pairs influences the speed of expansion of the 
fireshell. However, in this section we neglect the rate equation and discuss those 
points common among all the models being considered. 

First of all, let us return to Fig. 23. For almost all values of the r\ parameter 
the gamma factor is determined by the gas (i.e., plasma or baryons) admixture 
according to (197); consider this case in what follows. For given initial energy and 
radius this temperature depends only on 77, so one can write: 

-- 77+1 - 71+1 (212) 



{%)V + 1 arji+1 



where 



a = 2.1-10 3 V 2 7e -°- 5 . (213) 
From this formula one has immediately the two asymptotic regimes, namely: 



7 |.M_1,,<W, (214) 

[ a^/rj' 'I ^ '' max 

Notice, that the constant a is an extremely small number, so that after obtaining 
a precise value of r? ma x by equating to zero the derivative of function (212), one can 
expand the result in a Taylor series and get to the lowest order in a the result: 



?7ma X - ( - ) - 2, (215) 



7max — T(^?max) — ^ 



1+1^" 

a 



(216) 



In particular, in the case shown in Fig. 24 one has ?7max = 2.8-10 5 , 7max = 9.3-10 4 
while according to (208), r m = rj max = 1.75 • 10 5 . Clearly, our result is much more 
accurate. Actually, the value T m in (208) is obtained from equating asymptotes in 
(214) and there exists the following relation: 



a=Mr 1/2 - (217) 

Now we are ready to explain why the observed temperature (and consequently 
the observed energy) does not depend on rj in the region rj max < 77 < rj p . From the 
second line in (214) it follows that the gamma factor in this region behaves like 
7 oc rj^ 1 / 2 , while T esc oc 77 1 / 2 . These two exactly compensate each other leading to 
the independence of the observed quantities on rj in this region. This remains the 
same for rj > r\ p also, since here T esc = T p =const and from (212) 7 =const. 
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Nakar, Piran and Sari revision. Recently a revision of the fireball model was 
made in Ref. 172. These authors presented a new diagram for the final Lorentz 
gamma factor and for the energy budget of the fireball. Their work was motivated 
by the observation of giant flares with the subsequent afterglow spreading up to the 
radio region with a thermal spectrum. They concluded that the fireball has to be 
loaded by either baryons or magnetic field, and cannot be only a pure e ± , 7 plasma 
in order to have 10 -3 of the total energy radiated in the giant flare. 

In analogy with cosmology the authors define the number density of pairs which 
survive because the expansion rate becomes larger than the annihilation rate f which 
gives the condition 

n± « — ^. (218) 
Then, recalling (200), if we want to estimate the number of pairs it turn out to be 

where we identify A = ct in (205). In Ref. 172 the authors obtained the third power 
of the ratio of temperatures which influences all of their subsequent results. 

Having reached the conclusion that the afterglow cannot be obtained as the 
result of interaction of the e ± ,7 plasma with the CircumBurst Medium (CBM), 
the authors turn to baryonic loading considerations. They attempt to define critical 
values of the loading parameter r] finding in general 4 such values 8 , in particular: 

m = T ^- s (^)\ (220) 



4TtR ctm e c 2 \ T 

E a T ( T± 
ATrRoctrripC 2 \ Tq 



m= ._:ZL -A ^) . (221) 



1/4 



*-Usss?J - ,222) 

We recall that the first two quantities are based on the formula for N± and 
should contain factors of (-5^) instead. 

The first 'critical' value, 771, comes from the condition N p m p = N±m e , where 
N p = m E °i is just the number of protons in the plasma admixture. It does not 
correspond to any critical change in the physics of the phenomena; for instance it 
cannot be interpreted as the equality of masses (equal inertia) of pairs and baryons 
since the former is mainly due to their total energy E±, while the latter to their 
rest mass N p m p . This value is, however, close to the one defined above r\\ ss r\ p . 



This effect is accounted for automatically in our approach where rate equations for pairs include 
an expansion term. 

g The last value r?4 corresponds to the case of heavy loading where spreading of the expanding shell 
is observed, and is not considered here. 
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The second 'critical' value, 772, corresponds to the condition N p = N±, namely 
equality of the numbers of protons and pairs. It is also incorrectly interpreted 
as equal contributions to the Thompson scattering. In fact, the cross-section for 

Thompson scattering of protons contains an additional factor ( — ) with respect 

\ m p J 

to the usual formula for electrons. 

The definition of the third 'critical' value, 773, is not clear, but an important 
feature is its closeness to the critical value r\ c quoted above. 

Assuming adiabatic conditions (200) the authors present the new diagram for 
the final gamma factor and energy budget of the pair-baryonic plasma at trans- 
parency. In fact, this diagram, shown in our Fig. 27 by the dashed curve for the 
parameter B, is very similar to the one obtained in Ref. 173, which considered the 
hydrodynamics of relativistic , 7 winds. That problem is very different from ours, 
because of different boundary conditions' 1 . In particular, in the wind energy conser- 
vation, (142) does not hold; the reason is that constant energy (mass) supply takes 
place parametrized in Ref. 173 by E (M). In that paper in fact authors present the 
diagram for the asymptotic value of the Lorentz gamma factor depending on the 
ratio which is very different from the quantity r\. 

Surprisingly, the fundamental result about the presence of a maximum in the 
diagram for the gamma factor on Fig. 27 which was found by the same authors 
previously in Ref. 164 (see Fig. 24) that comes from the energy conservation (202) 
is ignored in 172. It can be understood in the following way. For small loading (small 
B) the more baryons are present in the plasma the larger the number density of 
corresponding electrons becomes, and the larger the optical depth is. Therefore, 
transparency is reached later, which gives a larger gamma factor at transparency. 
On the other hand, for heavy baryon loading (relatively large B) the more baryons 
are present, the more inertia the plasma has, and by energy conservation, the less 
the final gamma factor has to be. 

8.2.5. Significance of the rate equation 

The rate equation describes the evolution of the number densities for electrons 
and positrons. In analytic models it is assumed that pairs are annihilated instantly 
when the transparency condition is fulfilled. Moreover, the dynamics of expansion is 
influenced by the electron-positron energy density as can be seen from (190)-(194). 
Therefore, it is important to clarify whether neglecting the rate equation is a crude 
approximation or not. 

Using Eq. (202) one can obtain the analytic dependence of the energy emitted 
at the transparency point on parameter B and we compare them in Fig. 26. 

We also show the difference between numerical results based on integration of 
Eqs. (190-194) and analytic results from the Shemi and Piran model. The values of 



Note that the authors of Ref. 173 also use rate equations describing decoupling plasma from 
photons. 
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Fig. 26. Relative energy release in the form of photons emitted at the transparency point of the 
GRB in terms of the initial energy of the fireball depending on parameter B. Thick line represents 
numerical results and it is the same as in Fig. 22. The normal line shows the results for the 
analytic model of Ref. 164. The dashed line shows the difference between the exact numerical and 
approximate analytical results. 

the parameters are: /x = 10 3 and £ = 0.1 (which correspond to Eg = 2.87 • 10 54 ergs 
and i?o = 1-08 • 10 9 cm). One can see that the difference peaks at intermediate 
values of B. The crucial deviations, however, appear for large B, where analytical 
predictions for the observed energy are about two orders of magnitude smaller than 
the numerical ones. This is due to the difference in predictions of the radius of the 
fireshell at the transparency moment. In fact, the analytical model overestimates this 
value by about two orders of magnitude for B = 10~ 2 . So for large B with correct 
treatment of pair dynamics the fireshell becomes transparent earlier comparing to 
the analytical treatment. 

At the same time, the difference between numerical and analytical results for 
gamma factor is significant for small B as illustrated at Fig. 27. While both results 
coincide for B > 10~ 4 there is a constant difference for the range of values 10~ 8 < 
B < 10~ 4 and asymptotic constant values for the gamma factor are also different. 
Besides, this asymptotic behavior takes place for larger values of B in disagreement 
with analytical expectations. Thus the acceleration of the fireshell for small B is 
larger if one accounts for pair dynamics. 

It is clear that the error coming from neglecting the rate equation is significant. 
This implies that the simple analytic model of Shemi and Piran gives only a qual- 
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Fig. 27. Relativistic gamma factor when transparency is reached. The thick line denotes exact 
numerical results, the normal line corresponds to the analytical estimate from the Shemi and Piran 
model, while the dotted line denotes the asymptotic value of the baryonic loading parameter. The 
dashed line shows the results of Nakar, Piran and Sari. 

itative picture of the fireshell evolution and in order to get the correct description 
of the fireshell one cannot neglect the rate equation. 

Moreover, the difference between the exact numerical model 106 ' 107 and the 
approximate analytical models 164 becomes apparent in various physical aspects, 
namely in predictions of the radius of the shell when it reaches transparency, the 
gamma factor at transparency and the ratio between the energy released in the 
form of photons and that converted into kinetic form. The last point is crucial. It is 
assumed in the literature that the entire initial energy of the fireshell gets converted 
into kinetic energy of the shell during adiabatic expansion. Indeed, taking a typical 
value of the parameter B like 10 -3 we find that according to the Shemi and Piran 
model wc have only 0.2% of the initial energy left in the form of photons. However, 
exact numerical computations 106,107 give 3.7% for the energy of photons radiated 
when the fireshell reaches transparency, which is a significant value and it cannot 
be neglected. 

To summarize the above discussion, we present the results of this survey in the 
Table 1. It is important to notice again that comparing to the simplified analytic 
treatment, accounting for the rate of change of electron-positron pairs densities 
gives quantitatively different results for the ratio of kinetic versus photon energies 
produced in the GRB and the gamma factor at the transparency moment, which in 
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Table 1. Comparison of different models for fireballs. 





Rufnni et al. 


Shemi, Piran 


Piran, et al. 


Meszaros, et al. 


conservation: 










EM, 


yes 


yes 


yes 


yes 


bar. number, 


yes 


not consider 


yes 


yes 


rate equation 


yes 


no 


no 


no 


const width 


justify 


not consider 


justify 


in part 


approx. 










model 7(r) 


num./analyt. 


no 


num./analyt. 


num./analyt. 


model 7(77) 


num. 


analyt. 


not consider 


analytic 



turn leads to different afterglow properties. Therefore, although analytical models 
presented in sections 8.2.4 and 8.2.4 agree and give the correct qualitative descrip- 
tion of the fireshell, one should use the numerical approach described in section 
8.2.3 in order to compare the theory and observations. 




■iqO I Z I 1 I 1 I 1 I 1 I 1 I 1 

1Q 8 1Q 9 1(J 10 10 11 1Q 12 1Q 13 1Q 14 

Fireshell radius (cm) 

Fig. 28. The Lorentz factor as a function of the radius for selected values of the baryonic loading 
parameter. Filled circles show the radius at the moment of transparency. Squares show the moment 
of departure from thermal distributions of electron-positron pairs. The star denotes the moment 
when the temperature of the fireshell equals 511 keV. 



Reaching transparency. To demonstrate the richness of physical phenomena as- 
sociated with gamma-ray bursts we provide below several illustrations. We calculate 
the temperature in the comoving and laboratory frames of the plasma as well as 
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the Lorentz factor with the code described in Ref. 106,107 for different values of 
the parameter B. We show our results in Fig. 28 and 29. From Fig. 29 it is clear 
that during early phases of the expansion the temperature decreases down to 0.511 
MeV, as shown by the star in Fig. 28, and the ratio of electron-positron pairs to 
photons becomes exponentially suppressed. However, because of the accelerated ex- 
pansion the apparent temperature in the observer's frame remains almost constant, 
see Fig. 29. Then, after collision with a baryonic remnant after which the Lorentz 
factor decreases, the plasma continues to expand. At a certain moment, shown by 
squares in Fig. 28, a departure from a thermal distribution of electron-positron pairs 
occurs, due to the fact that the rate of the reaction 77 — ► e + e~ becomes smaller 
than the expansion rate. From that moment electron-positron pairs freeze out, anal- 
ogous to what happens in the early Universe. Finally, transparency is reached, as 
denoted by circles in Fig. 28. For high values of the baryonic loading B > 10~ 4 the 
comoving temperature decreases in the late stages of expansion in the same way as 
the observed temperature. 

Total opacity due to pair production and due to Compton scattering is shown 
in Fig. 30, while the ratio of separate contributions, i.e., opacity due to pairs and 
baryons is shown in Fig. 31. The expansion starts with t» 1 and the optical depth 
starts to decrease. Then, after collision with the baryonic remnant containing also 
associated electrons, the opacity may increase, but only for large baryonic loading 
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Fig. 30. Total opacity of the fireshell. 



B > 1CP 4 . The change of exponential decrease into a power law decrease seen in 
Fig. 30 corresponds to the departure of distributions of electrons and positrons from 
thermal ones. Finally, as one can see from Fig. 31 at the late stage of expansion 
of the fireshell the opacity is dominated by pairs for B < 10~ 7 and by Compton 
scattering for B > 10~ 7 . 



8.3. The afterglow 

After reaching transparency and the emission of the P-GRB, the accelerated bary- 
onic matter (the ABM pulse) interacts with the CBM and gives rise to the afterglow 
(see Fig. 32). Also in the descriptions of this last phase many differences exist be- 
tween our treatment and the other ones in the current literature (see next sections) . 

We first look to the initial value problem. The initial conditions of the afterglow 
era are determined at the end of the optically thick era when the P-GRB is emitted. 
As recalled in the last section, the transparency condition is determined by a time 
of arrival tf~ , a value of the gamma Lorentz factor 7 , a value of the radial coordinate 
r , and the amount of baryonic matter Mb , all of which are only functions of the 
two parameters Ed ya and B (see Eq. (119)). 

This connection to the optically thick era is missing in the current approach in 
the literature which attributes the origin of the "prompt radiation" to an unspecified 
inner engine activity (see Ref. 153 and references therein). The initial conditions at 
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Fireshell radius (cm) 

Fig. 31. The ratio of the opacity due to electron-positron pairs to the opacity due to electrons 
associated with baryons in the fireshell. 



tttt 




Fig. 32. The GRB afterglow phase is represented here together with the optically thick phase 
(see Fig. 20). The value of the Lorentz gamma factor is given here from the transparency point all 
the way to the ultrarelativisitc, relativistic and nonrelativistic regimes. Details in Ref. 152. 

the beginning of the afterglow era are obtained by a best fit of the later parts of 
the afterglow. This approach is quite unsatisfactory since the theoretical treatments 
currently adopted in the description of the afterglow are not appropriate. 174-176 The 
fit which uses an inappropriate theoretical treatment leads necessarily to the wrong 
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conclusions as well as, in turn, to the determination of incorrect initial conditions. 

The order of magnitude estimate usually quoted for the characteristic time scale 
to be expected for a burst emitted by a GRB at the moment of transparency at the 
end of the optically thick expansion phase is given by r ~ GM/c 3 . For a 10M Q black 
hole this will give ~ 10~ 3 s. There are reasons today not to take seriously such an 
order of magnitude estimate. 177 In any case this time is much shorter than the ones 
typically observed in "prompt radiation" of the long bursts, from a few seconds all 
the way to 10 2 s. In the current literature (see e.g. Ref. 153 and references therein), 
in order to explain the "prompt radiation" and overcome the above difficulty it 
has been generally assumed that its origin should be related to a prolonged "inner 
engine" activity preceding the afterglow which is not well identified. 

To us this explanation has always appeared logically inconsistent since there 
remain to be explained not one but two very different mechanisms, independent of 
each other, of similar and extremely large energetics. This approach has generated 
an additional very negative result: it has distracted everybody working in the field 
from the earlier very interesting work on the optically thick phase of GRBs. 

The way out of this dichotomy in our model is different: 1) indeed the optically 
thick phase exists and is crucial to the GRB phenomenon and terminates with a 
burst: the P-GRB; 2) the "prompt radiation" follows the P-GRB; 3) the "prompt 
radiation" is not a burst: it is actually the temporally extended peak emission of 
the afterglow (E-APE). The observed structures of the prompt radiation can all be 
traced back to inhomogeneities in the interstellar medium (see Fig. 33 and Ref. 129). 

This approach was first tested on GRB 991216. Both the relative intensity and 
time separation of the P-GRB and the afterglow were duly explained (see Fig. 33) 
choosing a total energy of the plasma E to ± — Ed ya — 4.83 x 10 53 erg and a baryon 
loading B = 3.0 x 1CT 3 (see Ref. Ref. 127,129,152,178). Similarly, the temporal 
substructure in the prompt emission was explicitly shown to be related to the CBM 
inhomogeneities (see the following sections). 

Following this early analysis and the subsequent ones on additional sources, it 
became clear that the CBM structure revealed by our analysis is quite different from 
the traditional description in the current literature. Far from considering analogies 
with shock wave processes developed within a fluidodynamic approach, it appears 
to us that the correct CBM description is a discrete one, composed of overdense 
"blobs" of typical size Ai? ~ 10 14 cm widely spaced in underdense and inert regions. 

We can then formulate the second paradigm, the interpretation of the burst 
structure (IBS) paradigm, 127 which covers three fundamental issues leading to the 
unequivocal identification of the canonical GRB structure: 

a) the existence of two different components: the P-GRB and the afterglow related 
by precise equations determining their relative amplitude and temporal sequence 
(see Fig. 34, Ref. 152 and next section); 

b) what in the literature has been addressed as the "prompt emission" and consid- 
ered as a burst, in our model is not a burst at all — instead it is just the emission 
from the peak of the afterglow (see the clear confirmation of this result by the Swift 
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Detector arrival time (t B a ) (s) 



Fig. 33. The detailed features of GRB 991216 shown by our theoretical models are reproduced 
here: the P-GRB, the "prompt radiation" and what is generally called the afterglow. It is clear that 
the prompt emission observed by BATSE coincides with the extended afterglow peak emission (E- 
APE) and has been considered as a burst only as a consequence of the high noise threshold in the 
observations. The small precursor is identified with the P-GRB. Details in Ref. 127,129,152,178. 

data of e.g. GRB 050315 in the next sections and in Ref. 180,181); 
c) the crucial role of the parameter B in determining the relative amplitude of the 
P-GRB to the afterglow and discriminating between the short and the long bursts 
(see Fig. 35). Both short and long bursts arise from the same physical phenomena: 
the gravitational collapse to a black hole endowed with electromagnetic structure 
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Fig. 34. Bolometric luminosity of P-GRB and afterglow as a function of the arrival time. Details 
in Ref. 152. Reproduced and adapted from Ref. 179 with the kind permission of the publisher. 

and the formation of its dyadospherc. 

The fundamental diagram determining the relative intensity of the P-GRB and 
the afterglow as a function of the dimensionless parameter B is shown in Fig. 35. The 
main difference relates to the amount of baryonic matter engulfed by the electron- 
positron plasma in their optically thick phase prior to transparency. For B < 10~ 5 
the intensity of the P-GRB is larger and dominates the afterglow. This corresponds 
to the "genuine" short bursts. 182 For 10~ 5 < B < 10~ 2 the afterglow dominates the 
GRB. For B > 10~ 2 we may observe a third class of "bursts", eventually related 
to a turbulent process occurring prior to transparency. 107 This third family should 
be characterized by smaller values of the Lorentz gamma factors than in the case 
of the short or long bursts. 

Particularly enlightening for the gradual transition to the short bursts as a 
function of the B parameter is the diagram showing how the GRB 991216 bolometric 
light curve would scale changing the only the value of B (see Fig. 36). 

Moving from these two paradigms, and the prototypical case of GRB 991216, 
we have extended our analysis to a larger number of sources, such as GRB 
9702 28, 182 GRB 98 4 25, 183 - 184 GRB 03 3 29, 185 GRB 03 1 203, 186 GRB 050315, 180 
GRB 060218 187 which have led to a confirmation of the validity of our canonical 
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Fig. 35. Above: The energy radiated in the P-GRB (the red line) and in the afterglow (the green 
line), in units of the total energy of the dyadosphere {E^ ya ), are plotted as functions of the B 
parameter. Below: The arrival time delay between the P-GRB and the peak of the afterglow is 
plotted as a function of the B parameter for three selected values of Ed ya - 



GRB structure (see Fig. 37). In addition, progress has been made in our theoretical 
comprehension, which will be presented in the following sections. 
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Detector arrival time (s) 

Fig. 36. The bolometric luminosity of a source with the same total energy and CBM distribution 
of GRB 991216 is represented here for selected values of the B parameter, ranging from B = 10~ 2 
to B = 10~ 4 . The actual value for the GRB 991216 is B = 3.0 X 10 -3 . As expected, for smaller 
values of the B parameter the intensity of the P-GRB increases and the total energy of the 
afterglow decreases. What is most remarkable is that the luminosity in the early part of the 
afterglow becomes very spiky and the peak luminosity actually increases. 

8.3.1. The simple model for the afterglow 

When electron-positron pairs are annihilated and the transparency condition r ~ 1 
is satisfied, photons escape and the only ingredient of the fireshell which is left 
is a relativistically expanding shell of protons and electrons. The latter are not 
important kinematically and one can say that the shell consists of baryons. This 
shell propagates in the interstellar medium sweeping up the cold gas. The constant 
width approximation is good enough to describe this process, however, we do not 
restrict ourselves to this case. The only requirement is that the collision is inelastic 
and the energy released in the collision is shared within the whole shell for a short 
time. 

Consider the collision process in the lab reference frame where the CBM is 
initially at rest. Assuming that the expanding baryons as well as CBM are cold, 
the total energy of the shell is E — Mbc 2 j = Mj with 7 being the Lorentz factor 
of the shell and M is the rest mass of the shell (together with its thermal energy) 
in units of energy. The mass-energy of the CBM swept up within an infinitesimal 
time interval is dm. The energy released during this process is dE. In our definitions 
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Fig. 37. Same as Fig. 35 with the values determined for selected GRBs. In order to determine 
the value of the B parameter and the total energy we have performed the complete fit of each 
source. In particular, for each source we have fit the observed luminosities in selected energy bands 
of the entire afterglow including the prompt emission. We have verified that in each source the 
hard-to-soft spectral evolution is correctly fit and we have compared the theoretically computed 
spectral lag with the observations. Where applicable, we have also computed the relative intensity 
and temporal separation between the P-GRB and the peak of the afterglow and compared these 
values with the observed ones. The absence of spectral lag in the P-GRB is automatically verified 
by our model. 



7 /3 = ^ - 1. Finally, the gamma factor after the collision becomes j+d-y. Rewrite 
energy and momentum conservation: 



M7 + dm = (M + dm + dE){~f + cfry), 



M V7 2 - 1 = (M + dm + dE) y / (j + dj) 2 - 1. 
This set of equations is equivalent to the one used in 152. From (224) we get: 



dE = -dm - M + Ma 



7 2 - 1 



(7 + dj) 2 - 1 ' 
Substituting the last equality into (223) we arrive at 

where 



(223) 
(224) 

(225) 

(226) 
(227) 
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Substituting these expressions into (225) we find 



dE = -dm- M(l-b). 



(228) 



Since dm is infinitcsimally small we can expand b in a series in dm/M 



6~ 1 



dm 
~M 



(229) 



so we finally obtain: 



dE = (j~ l)dm, 



(230) 




(231) 



These equations are used to describe collision of the baryonic remnant of the fireshell 
with the CBM. 107 ' 152 - 153 Note that it is a mistake to assume that dE = and use 
only (224) to obtain a differential equation for 7 as was done in Ref. 170. The total 
mass-energy swept up during this process is given by 



where 5 denotes the portion of the energy that is radiated. There are basically two 
approximations which follow from this expression, namely the fully radiative con- 
dition 5=1 and the adiabatic condition 5 = 0. With the adiabatic condition no 
energy is radiated and one must assume that the kinetic energy of the decelerating 
baryons is converted into radiation in dissipation processes in shocks which are sup- 
posed to form during the collision, 168,169 see also 153. With the radiative condition 
no additional mechanisms are required to describe the afterglow since it results from 
the emission of part of the energy released during inelastic collision of accelerated 
baryonic pulse with the CBM. 152 

8.3.2. Blandford and McKee radiative solution 

Eqs. (230), (231) can be found already in Ref. 167. However, the meaning of the 
quantities used there as well as the derivation are incorrect. 

Following Ref. 167 consider the relativistic shock (blast) wave resulting from 
an impulsive energy release, propagating into the homogeneous external medium 
with the Lorentz factor 7. Authors assume that the shock wave sweeps up external 
matter. They treat this swept-up matter as being in a thin, cold shell. 

Consider the reference frame where the shock is at rest. For an observer in this 
frame the external medium has relativistic Lorentz factor 7. The shock front is a 
spherical massless surface. For our observer there is a flow of matter, energy and 
momentum through this surface. We can calculate these quantities using (142) and 



dM = M + dm + (1 - 5)dE — M = dm + (1 — 5)dE, 



(232) 
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<f T 01 dS = -<f P1 2 l3dS = -4nR 2 p(3j 2 , (233) 
Js Js 

<j> T n dS = - <j> p( 7 2 - l)dS = -4nR 2 p(-/ 2 - 1), (234) 

p^dS = -AnR 2 p(3 1 , (235) 



(122). Recall that the external medium is cold so p <C p. Thus we get 

dEtot 
dt 
dP 
~dt 
dm 
~~dt ' 

where P is a radial momentum of external medium, E tot is its total energy. The 
kinetic energy Ek change is simply connected to its flux: 
dE k dE tot dm 2 

~df = ~d~T ~H±^ ^ 7(7 - 1} - (236) 

This formula coincides with Eq. (84) in Ref. 167. The next equation (85) is of the 
form of our (231). This implies that in order to derive this result these authors 
considered the full derivative of kinetic energy of the external medium with respect 
to the expanding shock wave, namely: 

dE k dm d'j 

^T = (7 - 1) ^ + m *' (237) 

and equate these two expressions. However, if one proceeds in the same manner 
with the radial momentum equation, one arrives at a strange result: they become 
inconsistent. The reason is the following. 

Eqs. (233), (234), (235) and (236) are of course correct. The next step is in doubt. 
To write down the total derivative means to assume that the Lorentz factor changes 
with time. This means that the properties of the external medium change. However, 
the problem considered by Blandford and McKee is an idealized one, namely the 
hydrodynamics of rclativistic shocks without consideration of any microphysics 1 . In 
this context the damping and dissipation effect cannot be incorporated directly and 
the shock wave generally speaking will propagate with constant velocity. Apart from 
difficulties in definition of the quantity m, the mass of external medium which can 
be thought of as infinite, the last term in (237) cannot exist and therefore there is 
no way to get (85) from (84) in 167. 

From the physical point of view, instead of a massless shock (even if it is present) 
a massive shell must be considered^. The only natural way to deal with this prob- 
lem is to consider interaction of this shell with a small shell of external medium. 
The problem is consequently equivalent to the interaction of two massive particles 
which is obvious to treat with the help of conservation equations instead of shock 
equations. 

Thus, equations (230), (231) can be derived only from energy and radial momen- 
tum conservation equations as was done in 188 and independently in Ref. 107. 



'Microphysical interactions responsible for changes in the blast wave can be Coulumb interactions 
within the plasma, ionization of external medium losses, interactions with magnetic field and 
various radiative processes. 

J In some approximation it can be considered as the relativistic piston problem. 
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9. The "canonical GRB" bolometric light curve 

We assume that the internal energy due to kinetic collision is instantly radiated away 
and that the corresponding emission is isotropic. Let Ae be the internal energy 
density developed in the collision. In the comoving frame the energy per unit of 
volume and per solid angle is simply 

( dE \ Ae 



\dVdn) o 4tt 



(238) 



due to the fact that the emission is isotropic in this frame. The total number of 
photons emitted is an invariant quantity independent of the frame used. Thus we 
can compute this quantity as seen by an observer in the comoving frame (which we 
denote with the subscript "o") and by an observer in the laboratory frame (which 
we denote with no subscripts). Doing this we find: 

dN~ ( dN~ \ , , n /nnns 
= A- 3 costf, (239) 



where d is the angle between the radial expansion velocity of a point on the fireshell 
surface and the line of sight, cosi? comes from the projection of the elementary 
surface of the shell on the direction of propagation and A = 7(1 — /3cos-#) is the 
Doppler factor introduced in the two following differential transformation 

dfl =dflx A- 2 (240) 

for the solid angle transformation and 

dt Q =dtx A- 1 (241) 

for the time transformation. An extra A factor comes from the energy transforma- 
tion: 

E = E x A (242) 
(see also Ref. 189). Thus finally we obtain: 

A -4 cos ,. (243) 



dtdndY, \dtdndEj G " 

Doing this we clearly identify ( dtdndu ) as ^ ne ener §y density in the comoving frame 
up to a factor ^ (see Eq. (238)). Then we have: 

4^ = I —vcosd A" 4 , (244) 
dtdn J shell An y ' 

where the integration in dT, is performed over the visible area of the ABM pulse 
at laboratory time t, namely with < $ < i? mal and dmax is the boundary of the 
visible region defined by: 

v 

cos d max = - . (245) 
c 
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Eq. (244) gives us the energy emitted toward the observer per unit solid angle and 
per unit laboratory time t in the laboratory frame. 

What we really need is the energy emitted per unit solid angle and per unit 
detector arrival time t d , so we must use the complete relation between t d a and t 
given by: 



f d a = (l + z) 



i - - (i) cos + — 

c c 



(246) 



where r* is the initial size of the fireshell. First we have to multiply the integrand 
in Eq. (244) by the factor (dt/dtf.) to transform the energy density generated per 
unit of laboratory time t into the energy density generated per unit arrival time 
t d a . Then we have to integrate with respect to e?£ over the EquiTemporal Surfaces 
(EQTS) corresponding to arrival time t d instead of the ABM pulse visible area at 
laboratory time t. The analog of Eq. (244) for the source luminosity in detector 
arrival time is then: 

dEL f Ae . a dt . . 

-r-r^-r = —v costf A~ 4 — T dZ . (247 

dtidSl J eqts 4tt dt* 1 ; 

It is important to note that, in the present case of GRB 991216, the Doppler factor 
A~ 4 in Eq. (247) enhances the apparent luminosity of the burst compared to the 
intrinsic luminosity by a factor which at the peak of the afterglow is in the range 
between 10 10 and 10 12 ! 

We are now able to reproduce in Fig. 34 the general behavior of the luminosity 
starting from the P-GRB to the latest phases of the afterglow as a function of 
the arrival time. It is generally agreed that the GRB afterglow originates from an 
ultrarelativistic shell of baryons with an initial Lorentz factor j Q ~ 200-300 with 
respect to the CBM (see e.g. Ref. 152,174 and references therein). Using GRB 991216 
as a prototype, in Ref. 126,127 we have shown how from the time varying bolometric 
intensity of the afterglow it is possible to infer the average density (n c b m ) = 1 
particle/cm 3 of the CBM in a region of approximately 10 17 cm surrounding the 
black hole giving rise to the GRB phenomenon. 

The summary of these general results are shown in Fig. 33, where the P-GRB, 
the emission at the peak of the afterglow in relation to the "prompt emission" and 
the latest part of the afterglow are clearly identified for the source GRB 991216. 
Details in Ref. 152. 

Summarizing, unlike treatments in the current literature (see e.g. Ref. 190,191 
and references therein), we define a "canonical GRB" light curve with two sharply 
different components (see Fig. 33): 1,127 ' 182 

(1) The P-GRB: it has the imprint of the black hole formation, a harder spectrum 
and no spectral lag. 177,192 

(2) The afterglow: it presents a clear hard-to-soft behavior; 180,186 ' 193 the peak 
of the afterglow contributes to what is usually called the "prompt emis- 

gion » 127,180,187 
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The ratio between the total time-integrated luminosity of the P-GRB (namely, its 
total energy) and the corresponding one for the afterglow is the crucial quantity for 
the identification of a GRBs' nature. Such a ratio, as well as the temporal separation 
between the corresponding peaks, is a function of the B parameter. 127 

When the P-GRB is the leading contribution to the emission and the afterglow 
is negligible we have a "genuine" short GRB. 127 This is the case where B < 10~ 5 
(see Fig. 37): in the limit B — > the afterglow vanishes (see Fig. 37). In the other 
GRBs, with 10~ 4 < B < 10~ 2 , the afterglow contribution is generally predominant 
(see Fig. 37; for the existence of the upper limit B < 1CT 2 see Ref. 107,187). Still, 
this case presents two distinct possibilities: 

• The afterglow peak luminosity is larger than the P-GRB peak luminosity. A 
clear example of this situation is GRB 991216, represented in Fig. 33. 

• The afterglow peak luminosity is smaller than the P-GRB one. A clear example 
of this situation is GRB 970228, represented in Fig. 53. 

The simultaneous occurrence of an afterglow with total time-integrated lumi- 
nosity larger than the P-GRB one, but with a smaller peak luminosity, is indeed 
explainable in terms of a peculiarly small average value of the CBM density, com- 
patible with a galactic halo environment, and not due to the intrinsic nature of the 
source (see Fig. 53). 182 Such a small average CBM density deflates the afterglow 
peak luminosity. Of course, such a deflated afterglow lasts much longer, since the 
total time-integrated luminosity in the afterglow is fixed by the value of the B pa- 
rameter (see above and Fig. 55). In this sense, GRBs belonging to this class are 
only "fake" short GRBs. This is GRB class identified by Ref. 194, which also GRB 
060614 belongs to, and which has GRB 970228 as a prototype. 182 

Our "canonical GRB" scenario, therefore, especially points out the need to dis- 
tinguish between "genuine" and "fake" short GRBs: 

• The "genuine" short GRBs inherit their features from an intrinsic property 
of their sources. The very small fircshell baryon loading, in fact, implies that 
the afterglow time-integrated luminosity is negligible with respect to the P-GRB 
one. 

• The "fake" short GRBs instead inherit their features from the environment. 
The very small CBM density in fact implies that the afterglow peak luminosity 
is lower than the P-GRB one, even if the afterglow total time-integrated lumi- 
nosity is higher. This deflated afterglow peak can be observed as a "soft bump" 
following the P-GRB spike, as in GRB 97 2 28, 182 GRB 060614 (Caito et al., in 
preparation), and the sources analyzed by Ref. 194. 

A sketch of the different possibilities depending on the fireshell baryon loading B 
and the average CBM density (n c b m ) is given in Fig. 38. 
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B < 10 s 

The total time-integrated P-GRB 
luminosity is larger than the afterglow 
one. 


► 


"Genuine" short GRBs: 

the P-GRB is the leading contribution to 
the emission and the afterglow is 
negligible 






n dim ~ 10 3 l/tm 3 
(compatible with a galactic 
halo environment) 


— ► 


"Fake" short GRBs: 

the P-GRB appears to be the leading 
contribution to the prompt emission 
because the afterglow is deflated by the 
small CBM density 
(e.g. GRB 970228, GRB 060614, etc.). 


10-"<B<10> 

The total time-integrated P-GRB 
luminosity is smaller than the 
afterglow one. 
















Normal ("long") GRBs: 

(e.g. GRB 991216, GRB 050315, etc.). 




„ ta -l#/cm' 



Fig. 38. A sketch summarizing the different possibilities predicted by the "canonical GRB" sce- 
nario depending on the fireshell baryon loading B and the average CBM density (n c f, m ). 



10. The spectra of the afterglow 

In our approach we focus uniquely on the X and gamma ray radiation, which ap- 
pears to be conceptually more predictable in terms of fundamental processes than 
the optical and radio emission. It is perfectly predictable by a set of constitutive 
equations, which leads to directly verifiable and very stable features in the spectral 
distribution of the observed GRB afterglows. In line with the observations of GRB 
991216 and other GRB sources, we assume in the following that the X and gamma 
ray luminosity represents approximately 90% of the energy of the afterglow, while 
the optical and radio emission represents only the remaining 10%. 

This approach differs significantly from the other ones in the current literature, 
where attempts are made to explain at once all the multi-wavelength emission in the 
radio, optical, X and gamma ray as coming from a common origin which is linked to 
boosted synchrotron emission. Such an approach has been shown to have a variety 
of difficulties 195 ' 196 and cannot anyway have the instantaneous variability needed 
to explain the structure in the "prompt radiation" in an external shock scenario, 
which is indeed confirmed by our model. 

Here the fundamental new assumption is adopted (see also Ref. 193) that the X 
and gamma ray radiation during the entire afterglow phase has a thermal spectrum 
in the comoving frame. The temperature is then given by: 

T s = [AE int / (47rr 2 Arcr^)] 1/4 , (248) 

where A_E; nt is the internal energy developed in the collision with the CBM in a 
time interval At in the co-moving frame, a is the Stcfan-Boltzmann constant and 

K = A eff /A ms , (249) 

is the ratio between the "effective emitting area" of the ABM pulse of radius r and 
its total visible area, which accounts for the CBM filamentary structure. 197 Due 
to the CBM inhomogencities the ABM emitting region is in fact far from being 
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homogeneous. In GRB 991216 such a factor is observed to be decreasing during the 
afterglow between: 3.01 x 10~ 8 > K > 5.01 x 10~ 12 . 193 

The temperature in the comoving frame corresponding to the density distribu- 
tion described in Ref. 129 is shown in Fig. 39. 



10"' 



10 




r (cm) 



Fig. 39. The temperature in the comoving frame of the shock front corresponding to the density 
distribution with the six spikes A,B,C,D,E,F presented in Ref. 197. The green line corresponds to 
an homogeneous distribution with n c j, m = 1. Details in Ref. 197. 



We are now ready to evaluate the source luminosity in a given energy band. The 
source luminosity at a detector arrival time t^, per unit solid angle dil and in the 
energy band [^1,^2] is given by: 152,193 

dE^ V2] f Ae „ dt 



-J EQTS 4. V C ° S * A ^W(^,T arr )dX, (250) 

where Ae = AEi n t/V is the energy density released in the interaction of the ABM 
pulse with the CBM inhomogeneities measured in the comoving frame, A = 7(1 — 
(v/c) cost?) is the Doppler factor, W {u\, f2,T arr ) is an "effective weight" required to 
evaluate only the contributions in the energy band 1^2], dS is the surface element 
of the EQTS at detector arrival time i„ on which the integration is performed (see 
also Ref. 197) and T arr is the observed temperature of the radiation emitted from 
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ci£: 

T arr = T s / [ 7 (1 - (v/c) COS 0) (1 + z)} . (251) 

The "effective weight" W {u\, V2,T arr ) is given by the ratio of the integral over 
the given energy band of a Planckian distribution at a temperature T arr to the total 
integral aT^ rr : 

W(u 1 ,u 2 ,T arr ) = -^- J p(T arr ,v)d{^—) , (252) 

where p (T arr , v) is the Planckian distribution at temperature T arr : 

P {Tarr, v) = (2 / ' h 3 ) hv / (e hv ^ kT ^ - l) (253) 

11. Application to GRB 031203: the first spectral analysis 

The first verification of the above theoretical framework came from the analysis 
of GRB 031203. GRB 031203 was observed by IBIS, on board the INTEGRAL 
satellite, 198 as well as by XMM 199 and Chandra 200 in the 2 - -10 kcV band, and 
by VLT 200 in the radio band. It appears as a typical long burst, 201 with a simple 

profile and a duration of w 40 s. The burst ffuencc in the 20 200 kcV band is 

(2.0±0.4) x 10~ 6 erg/cm 2 , 201 and the measured rcdshift is z = 0.106. 202 We analyze 
in the following the gamma ray signal received by INTEGRAL. The observations 
in other wavelengths, in analogy with the case of GRB 98 4 2 5, 184 ' 203 ' 204 could be 
related to the supernova event, as also suggested by Ref. 200, and they will be 
examined elsewhere. 

The INTEGRAL observations find a direct explanation in our theoretical model. 
We reproduce correctly the observed time variability of the prompt emission (see 
Fig. 40). 186 The radiation produced by the interaction with the CBM of the baryonic 
matter shell, accelerated by the fireshell, agrees with observations both for intensity 
and time structure. 

The progress in reproducing the X and 7— ray emission as originating from a 
thermal spectrum in the comoving frame of the burst 193 leads to the characterization 
of the instantaneous spectral properties which are shown to drift from hard to soft 
during the evolution of the system. The convolution of these instantaneous spectra 
over the observational time scale is in very good agreement with the observed power- 
law spectral shape. 

11.1. The initial conditions 

The best fit of the observational data leads to a total energy of the electron-positron 
plasma i?*± = 1.85 x 10 50 erg. Assuming a black hole mass M = 10M Q , we then 
have a black hole charge to mass ratio £ = 6.8 x 10~ 3 ; the plasma is created between 
the radii r\ = 2.95 x 10 6 cm and r2 = 2.81 x 10 7 cm with an initial temperature 
T = 1.52 MeV and a total number of pairs iV e ± = 2.98 x 10 55 . The amount of 
baryonic matter in the remnant is B = 7.4 x 10~ 3 . 
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After the transparency point and the P-GRB emission, the initial Lorentz 
gamma factor of the accelerated baryons is 70 = 132.8 at an arrival time at the 
detector t d a = 8.14 x 10~ 3 s and a distance from the Black Hole r D = 6.02 x 10 12 cm. 
The CBM parameters are: < n cbm >= 0.3 particle/cm 3 and < TZ >= 7.81 x 10~ 9 . 

11.2. The GRB luminosity in fixed energy bands 

The aim of our model is to derive from first principles both the luminosity in 
selected energy bands and the time resolved/integrated spectra. We recall that the 
luminosity in selected energy bands is evaluated integrating over the EQTSs 193 the 
energy density released in the interaction of the accelerated baryons with the CBM 
measured in the co-moving frame, duly boosted in the observer frame. The radiation 
viewed in the comoving frame of the accelerated baryonic matter is assumed to have 
a thermal spectrum and to be produced by the interaction of the CBM with the 
front of the expanding baryonic shell. 

In order to evaluate the contributions in the band [v\ , v 2 ] we have to multiply 
the bolometric luminosity by an "effective weight" W(vi,V2,T arr ), where T arr is 
the observed temperature. W(vi, v^, T arr ) is given by the ratio of the integral over 
the given energy band of a Planckian distribution at temperature T arr to the total 
integral aT^ rr . 193 The resulting expression for the emitted luminosity is Eq. (250). 

11.3. The "prompt emission" 

In order to compare our theoretical prediction with the observations, it is important 
to notice that there is a shift between the initial time of the GRB event and the 
moment in which the satellite instrument has been triggered. In fact, in our model 
the GRB emission starts at the transparency point when the P-GRB is emitted. If 
the P-GRB is under the threshold of the instrument, the trigger starts a few seconds 
later with respect to the real beginning of the event. Therefore it is crucial, in the 
theoretical analysis, to estimate and take this time delay into proper account. In 
the present case it results in At% — 3.5 s (see the bold red line in Fig. 40). In what 
follows, the detector arrival time is referred to the onset of the instrument. 

The structure of the prompt emission of GRB 031203, which is a single peak 
with a slow decay, is reproduced assuming an CBM which does not have a constant 
density but instead several density spikes with < n c t, m >= 0.16 particle/cm 3 . Such 
density spikes corresponding to the main peak are modeled as three spherical shells 
with width A and density contrast An/n: we adopted for the first peak A = 3.0 x 
10 15 cm and An/n = 8, for the second peak A = 1.0 x 10 15 cm and An/n = 1.5 and 
for the third one A = 7.0 x 10 14 cm and An/n = 1. To describe the details of the 
CBM filamentary structure we would require intensity vs. time information with 
an arbitrarily high resolving power. With the finite resolution of the INTEGRAL 
instrument, we can only describe the average density distribution compatible with 
the given accuracy. Only structures at scales of 10 15 cm can be identified. Smaller 
structures would need a stronger signal and/or a smaller time resolution of the 
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Fig. 40. Theoretically simulated light curve of the GRB 031203 prompt emission in the 20 — 200 
keV energy band (solid red line) is compared with the observed data (green points) from Ref. 201. 
The vertical bold red line indicates the time position of the P-GRB. 



detector. The three clouds here considered are necessary and sufficient to reproduce 
the observed light curve: a smaller number would not fit the data, while a larger 
number is unnecessary and would be indeterminable. 

The result (see Fig. 40) shows a good agreement with the light curve reported 
by Ref. 201, and it provides further evidence for the possibility of reproducing light 
curves with a complex time variability through CBM inhomogeneities. 129 ' 152 ' 178 



11.4. The instantaneous spectrum 

In addition to the the luminosity in fixed energy bands we can also derive the 
instantaneous photon number spectrum N(E). In Fig. 41 are shown samples of 
time-resolved spectra for five different values of the arrival time which cover the 
whole duration of the event. 

It is clear from this picture that, although the spectrum in the comoving frame 
of the expanding pulse is thermal, the shape of the final spectrum in the laboratory 
frame is clearly not thermal. In fact, as explained in Ref. 193, each single instanta- 
neous spectrum is the result of an integration of hundreds of thermal spectra over 
the corresponding EQTS. This calculation produces a nonthermal instantaneous 
spectrum in the observer frame (see Fig. 41). 

Another distinguishing feature of the GRBs spectra which is also present in these 
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Fig. 41. Five different theoretically predicted instantaneous photon number spectra N(E) for 
tJJ = 2, 6, 10, 14, 18 s are represented here (colored curves) together with their own temporal 
convolution (black bold curve). The shapes of the instantaneous spectra are not blackbodies due 
to the spatial convolution over the EQTS (see text). 

instantaneous spectra, as shown in Fig. 41, is the hard to soft transition during the 
evolution of the eve nt. 153 ' 195 ' 205 ' 206 In fact the peak of the energy distributions E p 
drift monotonically to softer frequencies with time (see Fig. 42). This feature ex- 
plains the change in the power-law low energy spectral index a 207 which at the 
beginning of the prompt emission of the burst (t° = 2 s) is a — 0.75, and progres- 
sively decreases for later times (see Fig. 41). In this way the link between E p and a 
identified by Ref. 205 is explicitly shown. This theoretically predicted evolution of 
the spectral index during the event unfortunately cannot be detected in this par- 
ticular burst by INTEGRAL because of the insufficient quality of the data (poor 
photon statistics, see Ref. 201). 

11.5. The time-integrated spectrum: comparison with the observed 
data 

The time-integrated observed GRB spectra show a clear power-law behavior. Within 
a different framework Shakura, Sunyaev and Zel'dovich (see e.g. Ref. 208 and ref- 
erences therein) argued that it is possible to obtain such power-law spectra from 
a convolution of many non-power-law instantaneous spectra evolving in time. This 
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Fig. 42. The energy of the peak of the instantaneous photon number spectrum N(E) is here 
represented as a function of the arrival time during the "prompt emission" phase. The clear hard 
to soft behavior is shown. 

result was recalled and applied to GRBs by Blinnikov et al. 209 assuming for the in- 
stantaneous spectra a thermal shape with a temperature changing with time. They 
showed that the integration of such energy distributions over the observation time 
gives a typical power-law shape possibly consistent with GRB spectra. 

Our specific quantitative model is more complicated than the one considered by 
Blinnikov et al.: 209 the instantaneous spectrum here is not a black body. Each in- 
stantaneous spectrum is obtained by an integration over the corresponding EQTS: 
it is itself a convolution, weighted by appropriate Lorentz and Doppler factors, of 
~ 10 6 thermal spectra with variable temperature. Therefore, the time-integrated 
spectra are not plain convolutions of thermal spectra: they are convolutions of con- 
volutions of thermal spectra (see Fig. 41). 

The simple power-law shape of the integrated spectrum is more evident if we 
sum tens of instantaneous spectra, as in Fig. 43. In this case we divided the prompt 
emission in three different time interval, and for each one we integrated on time 
the energy distribution. The resulting three time-integrated spectra have a clear 
nonthermal behavior, and still show the characteristic hard to soft transition. 

Finally, we integrated the photon number spectrum N(E) over the entire du- 
ration of the prompt event (see again Fig. 43): in this way we obtain a typical 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



91 



convolution 
INTEGRAL data 

s < tg < 5 s 
5 s < tg < 1 s 
10s<t^<20s 




10 



100 
Energy (keV) 



1000 



Fig. 43. Three theoretically predicted time-integrated photon number spectra N(E) are here 
represented for < tjj < 5 s, 5 < < 10 s and 10 < ij[ < 20 s (colored curves). The hard to soft 
behavior presented in Fig. 42 is confirmed. Moreover, the theoretically predicted time-integrated 
photon number spectrum N(E) corresponding to the first 20 s of the "prompt emission" (black 
bold curve) is compared with the data observed by INTEGRAL (green points). 201 ' 210 This curve is 
obtained as a convolution of 108 instantaneous spectra, which are enough to get a good agreement 
with the observed data. 



nonthermal power-law spectrum which turns out to be in good agreement with the 
INTEGRAL data 201 ' 210 and gives clear evidence of the possibility that the observed 
GRBs spectra originate from a thermal emission. 

The precise knowledge we have here acquired on GRB 031203 helps in clarifying 
the overall astrophysical system GRB 031203 - SN 20031w - the 2 - 10 keV XMM 

and Chandra data (see sections 17.3 and 17.5, where the late 2 10 keV XMM 

and Chandra data are also discussed) . 



12. Application to GRB 050315: the first complete light curve 
fitting 

The fit of GRB 991216 was specially important in showing a good agreement be- 
tween the bolomctric luminosity predicted by our theory and the observations, ev- 
idencing the clear separation between the P-GRB and the afterglow. The INTE- 
GRAL observations of GRB 031203 have been crucial in leading to the confirmation 
of our theoretical approach for the spectral shape in the prompt emission of GRBs. 
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It has been the welcome consequence of the Swift satellite to have the possibility 
of receiving high quality data continuously spanning from the early part of the af- 
terglow (the prompt emission) to the late phases of the afterglow, leading to the 
emergence of a standard GRB structure. Our first analysis of the Swift data came 
with GRB 050315. 

GRB 050315 211 had been triggered and located by the BAT instrument 212 - 213 
on board the Swift satellite 214 at 2005-March-15 20:59:42 UT. 215 The narrow field 
instrument XRT 216 ' 217 began observations <~ 80 s after the BAT trigger, one of the 
earliest XRT observations yet made, and continued to detect the source for <~ 10 
days. 211 The spectroscopic redshift has been found to be z = 1.949. 218 

We present here the results of the fit of the Swift data of this source in 5 energy 
bands in the framework of our theoretical model, pointing out a new step toward 
the uniqueness of the explanation of the overall GRB structure. We first recall the 
essential features of our theoretical model; then we fit the GRB 050315 observations 
by both the BAT and XRT instruments; we also present the instantaneous spectra 
for selected values of the detector arrival time ranging from 60 s (i.e., during the 
so called "prompt emission") all the way to 3.0 x 10 4 s (i.e., the latest afterglow 
phases). 

12.1. The fit of the observations 

The best fit of the observational data leads to a total energy of the black hole 
dyadosphere, generating the e ± plasma, E*± = 1.46 x 10 53 erg (the observational 
Swift Ei SO is > 2.62 x 10 52 erg), 211 so that the plasma is created between the 
radii r\ — 5.88 x 10 6 cm and = 1.74 x 10 8 cm with an initial temperature 
T = 2.05MeV and a total number of pairs N e+e - = 7.93 x 10 57 . The second 
parameter of the theory, the amount Mb of baryonic matter in the plasma, is found 
to be such that B = Mbc 2 / Ed ya — 4.55 x 10 -3 . The transparency point and the P- 
GRB emission occurs then with an initial Lorentz gamma factor of the accelerated 
baryons "f Q = 217.81 at a distance r = 1.32 x 10 14 cm from the black hole. 

12.1.1. The BAT data 

In Fig. 44 we represent our theoretical fit of the BAT observations in the three 
energy channels 15-25 keV, 25-50 keV and 50-100 keV and in the whole 15-350 
keV energy band. 

We have already recalled how in our model the GRB emission starts at the 
transparency point when the P-GRB is emitted; this instant of time is often different 
from the moment in which the satellite instrument triggers, due to the fact that 
sometimes the P-GRB is under the instrumental noise threshold or comparable 
with it. In order to compare our theoretical predictions with the observations, it 
is important to estimate and take into account this time shift. In the present case 
of GRB 050315 there has been observed 211 a possible precursor before the trigger. 
Such a precursor is indeed in agreement with our theoretically predicted P-GRB, 
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both in its isotropic energy emitted (which we theoretically predict to be Ep-qrb = 
1.98 x 10 51 erg) and its temporal separation from the peak of the afterglow (which 
we theoretically predicted to be Aif = 51 s). In Fig. 44a the blue line shows our 
theoretical prediction for the P-GRB in agreement with the observations. 

After the P-GRB emission, all the observed radiation is produced by the inter- 
action of the expanding baryonic shell with the interstellar medium. In order to 
reproduce the complex time variability of the light curve of the prompt emission 
as well as of the afterglow, we describe the CBM filamentary structure, for sim- 
plicity, as a sequence of overdense spherical regions separated by much less dense 
regions. Such overdense regions are nonhomogeneously filled, leading to an effective 
emitting area A e ff determined by the dimensionless parameter 7?.. 193,197 Clearly, in 
order to describe any detailed structure of the time variability an authentic three- 
dimensional representation of the CBM structure would be needed. However, this 
finer description would not change the substantial agreement of the model with the 
observational data. Anyway, in the "prompt emission" phase, the small angular size 
of the source visible area due to the relativistic beaming makes such a spherical 
approximation an excellent one (see also for details Ref. 129). 

The structure of the "prompt emission" has been reproduced assuming three 
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overdense spherical CBM regions with width A and density contrast An/(n): we 
chose for the first region, at r = 4.15 x 10 16 cm, A = f .5 x 10 15 cm and An/ (n) = 
5.17, for the second region, at r = 4.53 x 10 16 cm, A = 7.0 x 10 14 cm and An/ (n) = 
36.0 and for the third region, at r — 5.62 x 10 16 cm, A = 5.0 x 10 14 cm and An/ (n) = 
85.4. The CBM mean density during this phase is (n c b m ) = 0.81 particles/cm 3 and 
(TZ) = 1.4 x 10 -7 . With this choice of the density mask we obtain agreement with the 
observed light curve, as shown in Fig. 44. A small discrepancy occurs in coincidence 
with the last peak: this is due to the fact that at this stage the source visible area due 
to the relativistic beaming is comparable with the size of the clouds, therefore the 
spherical shell approximation should be properly modified by a detailed analysis of 
a full three-dimensional treatment of the CBM filamentary structure. Such a topic 
is currently under investigation (see also for details Ref. 129). Fig. 44 also shows 
the theoretical fit of the light curves in the three BAT energy channels in which the 
GRB has been detected (15-25 keV in Fig. 44b, 25-50 keV in Fig. 44c, 50-100 keV 
in Fig. 44d). 



12.1.2. The XRT data 




1CT 1 10° 10 1 10 2 10 3 10 4 10 5 10 6 

Detector arrival time (s) 



Fig. 45. Our theoretical fit (blue line) of the XRT observations (green points) of GRB 050315 in 
the 0.2—10 keV energy band. 211 The theoretical fit of the BAT observations (see Fig. 44a) in the 
15—350 keV energy band is also represented (red line). 
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The same analysis can be applied to explain the features of the XRT light curve 
in the afterglow phase. It has been recently pointed out 219 that almost all the GRBs 
observed by Swift show a "canonical behavior" : an initial very steep decay followed 
by a shallow decay and finally a steeper decay. In order to explain these features 
many different approaches have been proposed. 191 ' 219 ~ 221 In our treatment these 
behaviors are automatically described by the same mechanism responsible for the 
prompt emission described above: the baryonic shell expands in a CBM region, 
between r = 9.00 x 10 16 cm and r = 5.50 x 10 18 cm, which is at significantly 
lower density {(n cbm ) = 4.76 x 10~ 4 particles/cm 3 , (K) = 7.0 x 10~ 6 ) then the 
one corresponding to the prompt emission, and this produces a slower decrease of 
the velocity of the baryons with a consequently longer duration of the afterglow 
emission. The initial steep decay of the observed flux is due to the smaller number 
of collisions with the CBM. In Fig. 45 is represented our theoretical fit of the XRT 
data, together with the theoretically computed 15-350 keV light curve of Fig. 44a 
(without the BAT observational data in order not to overwhelm the picture too 
much). 

What is impressive is that no different scenarios need to be advocated in order to 
explain the features of the light curves: both the prompt and the afterglow emission 
are just due to the thermal radiation in the comoving frame produced by inelastic 
collisions with the CBM properly boosted by the relativistic transformations over 
the EQTSs. 

12.2. The instantaneous spectrum 

In addition to the the luminosity in fixed energy bands we can derive also the in- 
stantaneous photon number spectrum N(E) starting from the same assumptions. 
In Fig. 46 are shown samples of time-resolved spectra for eight different values of 
the arrival time which cover the entire duration of the event. It is clear from this 
picture that although the spectrum in the co-moving frame of the expanding pulse is 
thermal, the shape of the final spectrum in the laboratory frame is clearly nonther- 
mal. In fact, as we have recalled and explained in Ref. 193, each single instantaneous 
spectrum is the result of an integration of thousands of thermal spectra over the cor- 
responding EQTS. This calculation produces a nonthermal instantaneous spectrum 
in the observer frame (see Fig. 46). 

A distinguishing feature of the GRB spectra which is also present in these 
instantaneous spectra is the hard to soft transition during the evolution of the 
event. 153 ' 195,205 ' 206 In fact the peak of the energy distribution E p drifts monoton- 
ically to softer frequencies with time. This feature is linked to the change in the 
power-law low energy spectral index a, 207 so the correlation between a and E p 205 
is explicitly shown. 

It is important to stress that there is no difference in the nature of the spec- 
trum during the prompt and the afterglow phases: the observed energy distribution 
changes from hard to soft, with continuity, from the "prompt emission" all the way 
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Fig. 46. Eight theoretically predicted instantaneous photon number spectra N(E) are represented 
here for different values of the arrival time (colored curves). The hard to soft behavior is confirmed. 

to the latest phases of the afterglow. 

13. Problems with the definition of "long" GRBs evidenced by 
GRB 05015 

The confirmation by Swift of our prediction of the overall afterglow structure, and 
especially the coincidence of the "prompt emission" with the peak of the afterglow, 
opens a new problematic in the definition of the long GRBs. It is clear in fact 
that the identification of the "prompt emission" in the current GRB literature is 
not at all intrinsic to the phenomenon but is merely due to the threshold of the 
instruments used in the observations (e.g. BATSE in the 50-300 keV energy range, 
or BeppoSAX GRBM in 40-700 keV, or Swift BAT in 15-350 keV). As it is clear 
from Fig. 47, there is no natural way to identify in the source a special extension 
of the peak of the afterglow that is not the one purely defined by the experimental 
threshold. It is clear, therefore, that long GRBs, as defined until today, are just 
the peak of the afterglow and there is no way, as explained above, to define their 
"prompt emission" duration as a characteristic signature of the source. As the Swift 
observations show, the duration of the long GRBs has to coincide with the duration 
of the entire afterglow. A Kouveliotou-Tavani plot of the long GRBs, done following 
our interpretation which is clearly supported by the recent Swift data (see Fig. 47), 
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will present enormous dispersion on the temporal axis. 




We recall that in our theory both "short" and "long" GRBs originate from 
the same process of black hole formation. The major difference between the two 
is the value of the baryon loading parameter B (see Fig. 35). In the limit of small 
baryon loading, all the plasma energy is emitted at the transparency in the P-GRB, 
with negligible afterglow observed flux. For higher values of the baryon loading, the 
relative energy content of the P-GRB with respect to the afterglow diminishes (see 
Ref. 178 and references therein). 

14. Application to GRB 060218: the first evidence of the critical 
value of the baryon loading B ~ 10 -2 

GRB 060218 triggered the BAT instrument of Swift on 18 February 2006 at 03:36:02 
UT and has a T 90 = (2 1 00 ± 100) s. 222 The XRT instrument 222 ' 223 began observa- 
tions ~ 153 s after the BAT trigger and continued for ~ 12.3 days. 224 The source 
is characterized by a flat gamma ray light curve and a soft spectrum. 225 It has an 
X-ray light curve with a long, slow rise and gradual decline and it is considered 
an X-Ray Flash (XRF) since its peak energy occurs at E p — 4.9^03 keV. 226 It 
has been observed by the Chandra satellite on February 26.78 and March 7.55 UT 
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(t ~ 8.8 and 17.4 days) for 20 and 30 ks respectively. 227 The spectroscopic redshift 
has been found to be z — 0.033. 228,229 The corresponding isotropic equivalent en- 
ergy is Ei SO = (1.9 ± 0.1) x 10 49 erg 224 which sets this GRB as a low luminous one, 
consistent with most of the GRBs associated with supernovas. 230-232 

GRB 060218 is associated with SN2006aj whose expansion velocity is v ~ 
0.1c. 231,233-235 The host galaxy of SN2006aj is a low luminosity, metal poor star- 
forming dwarf galaxy 236 with an irregular morphology 237 similar to that of other 
GRBs associated with supernovas. 228,238 



14.1. The fit of the observed data 
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Fig. 48. GRB 060218 prompt emission: a) our theoretical fit (blue line) of the BAT observations in 
the 15-150 keV energy band (pink points); b) our theoretical fit (red line) of the XRT observations 
in the 0.3—10 keV energy band (green points). 226 



In this section we present the fit of our fireshell model to the observed data 
(see Figs. 48, 51). The fit leads to a total energy of the plasma i?*± = 2.32 x 
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10 50 erg, with an initial temperature T = 1.86 MeV and a total number of pairs 
A e ± = 1.79 x 10 55 . The second parameter of the theory, B = 1.0 x 10~ 2 , is the 
highest value ever observed and is close to the limit for the stability of the adiabatic 
optically thick acceleration phase of the fircshell (for further details see Ref. 107). 
The Lorentz gamma factor obtained solving the fireshell equations of motion 175 ' 176 
is 70 = 99.2 at the beginning of the afterglow phase at a distance from the progenitor 
r Q = 7.82 x 10 12 cm. It is much larger than 7 ~ 5 estimated by Ref. 239 and Ref. 240. 

In Fig. 48 we show the afterglow light curves fitting the prompt emission both 
in the BAT (15-150 keV) and in the XRT (0.3-10 keV) energy ranges, as expected 
in our "canonical GRB" scenario. 187 Initially the two luminosities are comparable 
to each other, but for a detector arrival time t d > 1000 s the XRT curves becomes 
dominant. The displacement between the peaks of these two light curves leads to 
a theoretically estimated spectral lag greater than 500 s in perfect agreement with 
the observations. 241 We obtain that the bolometric luminosity in this early part 
coincides with the sum of the BAT and XRT light curves (see Fig. 51) and the 
luminosity in the other energy ranges is negligible. 

We recall that at t d <~ 10 4 s there is a sudden enhancement in the radio luminos- 
ity and there is an optical luminosity dominated by the SN2006aj emission. 226 ' 227 ' 242 
Although our analysis addresses only the BAT and XRT observations, for r > 10 18 
cm corresponding to t d > 10 4 s the fit of the XRT data implies two new features: 
1) a sudden increase of the ft factor from ft = 1.0 x lO" 11 to ft = 1.6 x 10~ 6 , 
corresponding to a significantly more homogeneous effective CBM distribution (see 
Fig. 52b); 2) an XRT luminosity much smaller than the bolometric one (see Fig. 51). 
These theoretical predictions may account for the energetics of the enhancement of 
the radio and possibly optical and UV luminosities. Therefore, we identify two dif- 
ferent regimes in the afterglow, one for t d < 10 4 s and the other for t d a > 10 4 s. 
Nevertheless, there is a unifying feature: the determined effective CBM density de- 
creases with the distance r monotonically and continuously through both these two 
regimes from n c b m = 1 particle/cm 3 at r = r to n c \, m = 10~ 6 particle/cm 3 at 
r = 6.0 x 10 18 cm: n cbm oc r~ a , with 1.0 < a < 1.7 (sec Fig. 52a). 

Our assumption of spherical symmetry is supported by the observations which 
set for GRB 060218 an opening beaming angle larger than - 370.226,227,230,232 



14.2. The procedure of the fit 

The arrival time of each photon at the detector depends on the entire previous 
history of the fireshell. 126 Moreover, all the observables depends on the EQTS 174 ' 175 
which in turn depend crucially on the equations of motion of the fireshell. The 
CBM cngulfment has to be computed self-consistently through the entire dynamical 
evolution of the fireshell and not separately at each point. Any change in the CBM 
distribution strongly influences the entire dynamical evolution of the fireshell and, 
due to the EQTS structure, produces observable effects up to a much later time. 
For example if we change the density mask at a certain distance from the black 
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hole we modify the shape of the lightcurve and consequently the evolution changes 
at larger radii corresponding to later times. Anyway the change of the density is 
not the only problem to face in the fitting of the source, in fact first of all we have 
to choose the energy in order to have a Lorentz gamma factor sufficiently high to 
fit the entire GRB. In order to show the sensitivity of the fitting procedure I also 
present two examples of fits with the same value of B and a different value of E*± . 

The first example has an S*± = 1.36 x 10 50 erg . This fit was unsuccessful as 
we see from the Fig. 49, because the bolometric lightcurve is under the XRT peak 
of the afterglow. This means that the value of the energy chosen is too small to fit 
any data points after the peak of the afterglow. So we have to increase the value of 
the energy to a have a better fit. In fact the parameter values have been found with 
various attempts in order to obtain the best fit. 

The second example is characterized by £?'±* = 1.61 x 10 50 erg and the all the 
data are fit except for the last point from 2.0 x 10 2 s to the end (see Fig. 50). We 
attempt to fit these last points trying to diminish the R values in order to enhance 
the energy emission, but again the low value of the Lorentz gamma factor that in 
this case is 3 prevents the fireshell from expanding. So again in this case the value 
of the energy chosen is too small, but it is better than the previous attempt. In this 
case we increased the energy value by 24%, but it is not enough so we decide to 
increase 16%. 

So the final fit is characterized by the B = 1.0 x 10~ 2 and by the E*± — 
2.32 x 10 50 erg. With this value of the energy we are able to fit all the experimental 
points. 

14.3. The fireshell fragmentation 

GRB 060218 presents different peculiarities: the extremely long Tg , the very low 
effective CBM density decreasing with the distance and the largest possible value 
of B = 10~ 2 . These peculiarities appear to be correlated. Following Ref. 184, we 
propose that in the present case the fireshell is fragmented. This implies that the 
surface of the fireshell does not increase any longer like r 2 but like r@ with (i < 2. 
Consequently, the effective CBM density n c b m is linked to the actual one by: 

n cbm =TZ shell n^ m , with n she ii = (r*/r) a , (254) 

where r* is the starting radius at which the fragmentation occurs and a = 2 — (3 (see 
Fig. 52a). For r* — r Q we have n"^ = 1 particles/cm 3 , as expected for a "canonical 
GRB" 1 and in agreement with the apparent absence of a massive stellar wind in 
the CBM. 227 ' 242,243 

The TZ parameter defined in Eq. (255) has to take into account both the effect 
of the fireshell fragmentation (TZ s heii) and the effective CBM porosity (TZ c b m ): 

TZ = U s hell x TZ cbm ■ (255) 

The phenomenon of the dumpiness of the ejecta, whose measure is the filling 
factor, is an aspect well known in astrophysics. For example, in the case of Novae 
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Fig. 49. GRB 060218 light curves with E l °l = 1.36 X 10 5U erg: our theoretical fit (blue line) 
of the 15-150 keV BAT observations (pink points), our theoretical fit (red line) of the 0.3—10 
keV XRT observations (green points) and the 0.3-10 keV Chandra observations (black points) 
are represented together with our theoretically computed bolometric luminosity (black line) (Data 
from: Ref. 226,227). 



the filling factor has been measured to be in the range 10 2 -10 5 . 244 Such a filling 
factor coincides, in our case, with TZ s heii- 



14.4. Binaries as progenitors of GRB-supernova systems 

The majority of the existing models in the literature appeal to a single astrophysical 
phenomenon to explain both the GRB and the SN ( "collapsar" , see e.g. Ref. 245). 
On the contrary, a distinguishing feature of our theoretical approach is to distinguish 
between the supernova and the GRB process. The GRB is assumed to occur during 
the formation process of a black hole. The supernova is assumed to lead to the 
formation of a neutron star (NS) or to a complete disruptive explosion without 
remnants and in no way to the formation of a black hole. In the case of SN2006aj 
the formation of such a NS has been actually inferred by Ref. 246 because of the 
large amount of 58 Ni (0.05M©). Moreover the significantly small initial mass of the 
supernova progenitor star M ss 20M Q is expected to form a NS rather than a black 
hole when its core collapses. 236 ' 246 ~ 248 In order to fulfill both the above requirement, 
we assume that the progenitor of the GRB and the supernova consists of a binary 
system formed by a NS close to its critical mass collapsing to a black hole, and a 
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Fig. 50. GRB 060218 light curves with E*^ = 1.61 X 10 50 erg: our theoretical fit (blue line) 
of the 15-150 keV BAT observations (pink points), our theoretical fit (red line) of the 0.3—10 
keV XRT observations (green points) and the 0.3-10 keV Chandra observations (black points) 
are represented together with our theoretically computed bolometric luminosity (black line). Data 
from: Ref. 226,227. 



companion star evolved out of the main sequence originating the supernova. The 
temporal coincidence between the GRB and the supernova phenomenon is explained 
in terms of the concept of "induced" gravitational collapse. 128,184 There is also the 
distinct possibility of observing the young born NS out of the supernova (see e.g., 
Ref. 184 and references therein). 

It has been often proposed that GRBs associated with Ib/c supernovas, at a 
smaller redshift 0.0085 < z < 0.168 (see e.g. Ref. 249 and references therein), 
form a different class, less luminous and possibly much more numerous than the 
high luminosity GRBs at higher redshift. 200,233,246 ' 249 Therefore they have been 
proposed to originate from a separate class of progenitors. 230 ' 231 In our model this 
is explained by the nature of the progenitor system leading to the formation of the 
black hole with the smallest possible mass: the one formed by the collapse of a just 
overcritical NS. 184 

The recent observation of GRB 060614 at z = 0.125 without an associated 
supernova 250,251 gives strong support to our scenario, alternative to the collapsar 
model. Also in this case the progenitor of the GRB appears to be a binary system 
composed of two NSs or a NS and a white dwarf. 



April 17, 2008 15:40 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



106 
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Fig. 51. GRB 060218 complete light curves: our theoretical fit (blue line) of the 15-150 keV 
BAT observations (pink points), our theoretical fit (red line) of the 0.3-10 keV XRT observations 
(green points) and the 0.3—10 keV Chandra observations (black points) 226,227 are represented 
together with our theoretically computed bolometric luminosity (black line). In this case we have 
= 2.32 X 10 50 ergs. 

14.5. Conclusions on GRB 060218 

GRB 060218 presents a variety of peculiarities, including its extremely large T90 
and its classification as an XRF. Nevertheless, a crucial point of our analysis is that 
we have successfully applied to this source our "canonical GRB" scenario. 

Within our model there is no need for inserting GRB 060218 in a new class 
of GRBs, such as the XRFs, alternative to the "canonical" ones. This same point 
recently received strong observational support in the case of GRB 06 02 1 8 241 and a 
consensus by other models in the literature. 239 

The anomalously long T 90 led us to infer a monotonic decrease in the CBM 
effective density giving the first clear evidence for the fragmentation of the fireshell, 
which indeed was predicted for values of the baryon loading B > 10~ 2 . For GRB 
060218 there is no need within our model for a new or unidentified source such as 
a magnetar or a collapsar. 

GRB 060218 is the first GRB associated with a supernova with complete cover- 
age of data from the onset all the way up to ~ 10 6 s. This fact offers an unprece- 
dented opportunity to verify theoretical models on such a GRB class. For example, 
GRB 060218 fulfills the Ref. 252 relation unlike other sources in its same class. This 
is particularly significant, since GRB 060218 is the only source in such a class to 
have an excellent data coverage without gaps. We are currently examining if the 
missing data in the other sources of such a class may have a prominent role in their 
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Fig. 52. The CBM distribution parameters: a) the effective CBM number density (red line) 
monotonically decreases with the distance r following Eq. (254) (green line); b) the 1Z parameter 
vs. distance. 

non-fulfillment of the Ref. 252 relation. 253 

15. Application to GRB 970228: the appearance of "fake" short 
GRBs 

GRB 970228 was detected by the Gamma-Ray Burst Monitor (GRBM, 40-700 
keV) and Wide Field Cameras (WFC, 2-26 keV) on board BeppoSAX on Febru- 
ary 28.123620 UT. 254 The burst prompt emission is characterized by an initial 5 s 
strong pulse followed, after 30 s, by a set of three additional pulses of decreasing 
intensity. 254 Eight hours after the initial detection, the NFIs on board BeppoSAX 
were pointed at the burst location for a first target of opportunity observation and 
a new X-ray source was detected in the GRB error box: this is the first "after- 
glow" ever detected. 102 A fading optical transient has been identified in a position 
consistent with the X-ray transient, 255 coincident with a faint galaxy with redshift 
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z = 0.695. 256 Further observations by the Hubble Space Telescope clearly showed 
that the optical counterpart was located in the outskirts of a late-type galaxy with 
an irregular morphology. 257 

The BeppoSAX observations of GRB 970228 prompt emission revealed a dis- 
continuity in the spectral index between the end of the first pulse and the beginning 
of the three additional ones. 102 ' 206 ' 254 The spectrum during the first 3 s of the sec- 
ond pulse is significantly harder than during the last part of the first pulse, 206 ' 254 
while the spectrum of the last three pulses appear to be consistent with the late 
X-ray afterglow. 206 ' 254 This was soon recognized by Ref. 206,254 as pointing to an 
emission mechanism producing the X-ray afterglow already taking place after the 
first pulse. 

The simultaneous occurrence of an afterglow with total time-integrated lumi- 
nosity larger than the P-GRB one, but with a smaller peak luminosity, is indeed 
explainable in terms of a peculiarly small average value of the CBM density and not 
due to the intrinsic nature of the source. In this sense, GRBs belonging to this class 
are only "fake" short GRBs. We show that GRB 970228 is a very clear example 
of this situation. We identify the initial spikelike emission with the P-GRB, and 
the late soft bump with the peak of the afterglow. GRB 970228 shares the same 
morphology and observational features with the sources analyzed by Ref. 194 as 
well as with e.g. GRB 05 7 09, 258 GRB 05 7 24 259 and GRB 06 06 1 4. 260 Therefore, 
we propose GRB 970228 as a prototype for this new GRB class. 



15.1. The analysis of GRB 970228 prompt emission 

In Fig. 53 we present the theoretical fit of BeppoSAX GRBM (40-700 keV) and 
WFC (2-26 keV) light curves of GRB 970228 prompt emission. 254 Within our 
"canonical GRB" scenario we identify the first main pulse with the P-GRB and 
the three additional pulses with the afterglow peak emission, consistent with the 
above mentioned observations by Ref. 102 and Ref. 254. The last three such pulses 
have been reproduced assuming three overdense spherical CBM regions (see Fig. 54) 
with very good agreement (see Fig. 53). 

We therefore obtain for the two parameters characterizing the source in our 
model E t( ± = 1.45 x 10 54 erg and B = 5.0x 10~ 3 . This implies an initial e ± plasma 
created between the radii ri = 3.52 x 10 7 cm and r 2 = 4.87 x 10 8 cm with a 
total number of e ± pairs N e ± = 1.6 x 10 59 and an initial temperature T = 1.7 
MeV. The theoretically estimated total isotropic energy emitted in the P-GRB is 
Ep-GRB = 1.1%£ I *±* = 1.54x 10 52 erg, in excellent agreement with the one observed 
in the first main pulse (Epl GRB <~ 1.5 x 10 52 erg in 2 — 700 keV energy band, see 
Fig. 53), as expected due to their identification. After the transparency point at 
tq = 4.37 x 10 14 cm from the progenitor, the initial Lorentz gamma factor of the 
fircshcll is 70 = 199. On average, during the afterglow peak emission phase we have 
for the CBM (K) = 1.5 x 10~ 7 and (n cbm ) = 9.5 x 10~ 4 particles/cm 3 . This very 
low average value for the CBM density is compatible with the observed occurrence 
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Fig. 53. The "canonical GRB" light curve theoretically computed for the prompt emission of 
GRB 970228. BeppoSAX GRBM (40-700 keV, above) and WFC (2-26 keV, below) light curves 
(data points) are compared with the afterglow peak theoretical ones (solid lines). The onset of the 
afterglow coincides with the end of the P-GRB (represented qualitatively by the dotted lines). For 
this source we have B ~ 5.0 X 10~ 3 and (n c (, m ) ~ 10~ 3 particles/cm 3 . 



of GRB 970228 in its host galaxy's halo 255,257,261 and it is crucial in explaining the 
light curve behavior. 

The values of E l °± and B we determined are univocally fixed by two tight con- 
straints. The first one is the total energy emitted by the source all the way up to 
the latest afterglow phases (i.e., up to ~ 10 6 s). The second one is the ratio be- 
tween the total time-integrated luminosity of the P-GRB and the corresponding 
one of the entire afterglow (i.e., up to - 10 6 s). In particular, in GRB 970228 such 
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Fig. 54. The CBM density profile we assumed to reproduce the last three pulses of the 
GRB 970228 prompt emission (red line), together with its average value (n cbm ) = 9.5 X 10 — 4 
particles/cm 3 (green line). 



a ratio turns out to be ~ 1.1% (see Fig. 37). However, the P-GRB peak luminosity 
actually turns out to be much more intense than the afterglow one (see Fig. 53). 
This is due to the very low average value of the CBM density (n cbm ) = 9.5 x 10~ 4 
particles/cm 3 , which produces a less intense afterglow emission. Since the afterglow 
total time-integrated luminosity is fixed, such a less intense emission lasts longer 
than what we would expect for an average density (n c b m ) ~ 1 particles/cm 3 . 



15.2. Reseating the CBM density 

We present now an explicit example in order to probe the crucial role of the average 
CBM density in explaining the relative intensities of the P-GRB and of the afterglow 
peak in GRB 970228. We keep fixed the basic parameters of the source, namely the 
total energy 2?*°' and the baryon loading B, therefore keeping fixed the P-GRB 
and the afterglow total time-integrated luminosities. Then we rescale the CBM 
density profile given in Fig. 54 by a constant numerical factor in order to raise its 
average value to the standard one (n c b m ) = 1 particle/cm 3 . We then compute the 
corresponding light curve, shown in Fig. 55. 

We notice a clear enhancement of the afterglow peak luminosity with respect 
to the P-GRB one in comparison with the fit of the observational data presented 
in Fig. 53. The two light curves actually cross at t% ~ 1.8 x 10 4 s since their total 
time-integrated luminosities must be the same. The GRB "rescaled" to (n c f, m ) — 
1 particle/cm 3 appears to be totally similar to, e.g., GRB 05 03 1 5 180 and GRB 
99 1 2 1 6. 152 ' 178 ' 193 
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Fig. 55. The theoretical fit of the BeppoSAX GRBM observations (red line, see Fig. 53) is com- 
pared with the afterglow light curve in the 40—700 keV energy band obtained rescaling the CBM 
density to (n c j, m ) = 1 particle/cm 3 keeping constant its shape and the values of the fundamental 
parameters of the theory and B (black line). The P-GRB duration and luminosity (blue line), 
depending only on -E^°* and B, are not affected by this process of rescaling the CBM density. 

It is appropriate to emphasize that, although the two underlying CBM density 
profiles differ by a constant numerical factor, the two afterglow light curves in Fig. 55 
do not. This is because the absolute value of the CBM density at each point affects 
in a nonlinear way all the following evolution of the fireshell due to the feedback 
on its dynamics. 176 Moreover, the shape of the surfaces of equal arrival time of the 
photons at the detector (EQTS) is strongly elongated along the line of sight. 175 
Therefore photons coming from the same CBM density region are observed over a 
very long arrival time interval. 

15.3. GRB 970228 and the Amati relation 

We turn now to the "Amati relation" 252,262 between the isotropic equivalent energy 
emitted in the prompt emission E iso and the peak energy of the corresponding time- 
integrated spectrum E P: i in the source rest frame. It has been shown by Ref. 252,262 
that this correlation holds for almost all the "long" GRBs which have a redshift and 
an E p ^ measured, but not for the ones classified as "short" . 262 If we focus on the 
"fake" short GRBs, namely the GRBs belonging to this new class, at least in one 
case (GRB 050 724) 259 it has been shown that the correlation is recovered if also the 
extended emission is considered. 263 

It clearly follows from our treatment that for the "canonical GRBs" with large 
values of the baryon loading and high (n c b m ), which presumably are most of the 
GRBs for which the correlation holds, the leading contribution to the prompt emis- 
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sion is the afterglow peak emission. The case of the "fake" short GRBs is completely 
different: it is crucial to consider separately the two components since the P-GRB 
contribution to the prompt emission in this case is significant. 

To test this scenario, we evaluated from our fit of GRB 970228 Ei SO and E p i only 
for the afterglow peak emission component, i.e., from = 37 s to if = 81.6 s. We 
found an isotropic energy emitted in the 2-400 keV energy band E iso — 1.5 x 10 52 
erg, and E p j = 90.3 keV. As it is clearly shown in Fig. 56, the sole afterglow 
component of GRB 970228 prompt emission is in perfect agreement with the Amati 
relation. If this behavior is confirmed for other GRBs belonging to this new class, 
this will reinforce our identification of the "fake" short GRBs. This result will also 
provide a theoretical explanation for the apparent absence of such a correlation for 
the initial spikclikc component in the different nature of the P-GRB. 



10000 



1000 



100 




Fig. 56. The estimated values for E„j and Ei ao obtained by our analysis (black dot) compared 
with the "Amati relation": 252 the solid line is the best fitting power law 262 and the dashed lines 
delimit the region corresponding to a vertical logarithmic deviation of 0.4. 262 The uncertainty in 
the theoretical estimated value for E p ^ has been assumed conservatively as 20%. 



15.4. Conclusions on GRB 970228 

We conclude that GRB 970228 is a "canonical GRB" with a large value of the 
baryon loading quite near to the maximum B ~ 10~ 2 (see Fig. 37). The difference 
with e.g. GRB 05 03 1 5 180 or GRB 991216 152 - 178 ' 193 is the low average value of the 
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CBM density (n c b m ) ~ 10 -3 particles/cm 3 which deflates the afterglow peak lumi- 
nosity. Hence, the predominance of the P-GRB, coincident with the initial spikelike 
emission, over the afterglow is just apparent: 98.9% of the total time-integrated 
luminosity is indeed in the afterglow component. Such a low average CBM density 
is consistent with the occurrence of GRB 970228 in the galactic halo of its host 
galaxy, 255,257 where lower CBM densities have to be expected. 261 

We propose GRB 970228 as the prototype for the new class of GRBs comprising 
GRB 060614 and the GRBs analyzed by Norris & Bonnell. 194 We naturally explain 
the hardness and the absence of spectral lag in the initial spikelike emission with the 
physics of the P-GRB originating from the gravitational collapse leading to the black 
hole formation. The hard-to-soft behavior in the afterglow is also naturally explained 
by the physics of the relativistic fireshell interacting with the CBM, clearly evidenced 
in GRB 03 1 203 186 and in GRB 05 03 1 5. 180 Also justified is the applicability of the 
Amati relation to the sole afterglow component. 262 ' 263 

This class of GRBs with z <~ 0.4 appears to be nearer than the other GRBs 
detected by Swift (z ~ 2. 3). 264 This may be explained by the afterglow peak lumi- 
nosity deflation. The absence of a jet break in those afterglows has been pointed 
out, 259 ' 265 consistently with our spherically symmetric approach. Their association 
with non-star-forming host galaxies appears to be consistent with the merging of a 
compact object binary. 266 ' 267 It is here appropriate, however, to caution on this con- 
clusion, since the association of GRB 060614 and GRB 970228 with the explosion 
of massive stars is not excluded. 250 ' 268 

Most of the sources of this class appear indeed not to be related to bright "Hy- 
pernovae" , to be in the outskirts of their host galaxies 267 and a consistent fraction of 
them are in galaxy clusters with CBM densities (n c b m ) ~ 10~ 3 particles/cm 3 . 269,270 
This suggests a spiraling out binary nature of their progenitor systems 271 made of 
neutron stars and/or white dwarfs leading to a black hole formation. 

Moreover, we verified the applicability of the Amati relation to the sole afterglow 
component in GRB 970228 prompt emission, in analogy with what happens for some 
of the GRBs belonging to this new class. In fact it has been shown by Ref. 262,263 
that the "fake" short GRBs do not fulfill the E p ^-Ei SO correlation when the sole 
spikelike emission is considered, while they do if the long soft bump is included. 
Since the spikclike emission and the soft bump contributions are comparable, it is 
natural to expect that the soft bump alone will fulfill the correlation as well. 

Within our "canonical GRB" scenario the sharp distinction between the P-GRB 
and the afterglow provide a natural explanation for the observational features of the 
two contributions. We naturally explain the hardness and the absence of spectral 
lag in the initial spikelike emission with the physics of the P-GRB originating from 
the gravitational collapse leading to the black hole formation. The hard-to-soft 
behavior in the afterglow is also naturally explained by the physics of the relativistic 
fireshell interacting with the CBM, clearly evidenced in GRB 03 1 203 186 and in GRB 
05 03 1 5. 180 Therefore, we expect naturally that the E p ^-E iso correlation holds only 
for the afterglow component and not for the P-GRB. Actually we find that the 
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correlation is recovered for the afterglow peak emission of GRB 970228. 

In the original work of Ref. 252,262 only the prompt emission is considered and 
not the late afterglow one. In our theoretical approach the afterglow peak emission 
contributes to the prompt emission and continues up to the latest GRB emission. 
Hence, the meaningful procedure within our model to recover the Amati relation 
is to look at a correlation between the total isotropic energy and the peak of the 
time-integrated spectrum of the whole afterglow. A first attempt to obtain such a 
correlation has already been performed using GRB 050315 as a template, giving 
very satisfactory results (...). 

16. The GRB-Supernova Time Sequence (GSTS) paradigm: the 
concept of induced gravitational collapse 



Tabic 2. 
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Note: see: a) Ref. 184; b) Ref. 186; c) Ref. 272; d) Ref. 180; e) Ref. 182; f) Ref. 178; g) see Ref. 239; h) Mazzali, P., private communication 
at MG11 meeting in Berlin, July 2006; i) evaluated fitting the URCAs with a power law followed by an exponentially decaying part; 
j) respectively Ref. 229, Ref. 273, Ref. 202, Ref. 274, Ref. 218, Ref. 275, Ref. 256, Ref. 276; k) respectively Ref. 223, Ref. 224, XRR is 
considered in Ref. 223, while XRF as suggested by Ref. 199, Ref. 223, Ref. 211, Ref. 256. 



Following the result of Ref. 273 who discovered the temporal coincidence of 
GRB 980425 and SN 1998bw, the association of other nearby GRBs with Type 
Ib/c SNe has been spectroscopically confirmed (see Tab. 2). The approaches in 
the current literature have attempted to explain both the supernova and the GRB 
as two aspects of the same astrophysical phenomenon. Hence, GRBs have been 
assumed to originate from a specially strong supernova process, a hypernova or a 
collapsar (see e.g. Ref. 245,277-279 and references therein). Both these possibilities 
imply very dense and strongly wind-like CBM structure. 

In our model we assumed that the GRB consistently originates from the gravita- 
tional collapse to a black hole. The supernova follows instead the complex pattern 
of the final evolution of a massive star, possibly leading to a neutron star or to 
a complete explosion but never to a black hole. The temporal coincidence of the 
two phenomena, the supernova explosion and the GRB, have then to be explained 
by the novel concept of "induced gravitational collapse", introduced in Ref. 128. 
We have to recognize that still today we do not have a precise description of how 
this process of "induced gravitational collapse" occurs. At this stage, it is more a 
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framework to be implemented by additional theoretical work and observations. Two 
different possible scenarios have been outlined. In the first version 128 we have con- 
sidered the possibility that the GRBs may have caused the trigger of the supernova 
event. For the occurrence of this scenario, the companion star had to be in a very 
special phase of its thermonuclear evolution and three different possibilities were 
considered: 

(1) A white dwarf, close to its critical mass. In this case, the GRB may implode 
the star enough to ignite thermonuclear burning. 

(2) The GRB enhances in an iron-silicon core the capture of the electrons on the 
iron nuclei and consequently decreases the Fermi energy of the core, leading to 
the onset of gravitational instability. 

(3) The pressure waves of the GRB may trigger a massive and instantaneous nuclear 
burning process leading to the collapse. 

More recently, 184 a quite different possibility has been envisaged: the supernova, 
originating from a very evolved core, undergoes explosion in presence of a companion 
neutron star with a mass close to its critical one. The supernova blast wave may 
then trigger the collapse of the companion neutron star to a black hole and the 
emission of the GRB (see Fig. 57). It is clear that, in both scenarios, the GRB and 
the supernova occur in a binary system. 

There are many reasons to propose this concept of "induced gravitational col- 
lapse" : 

(1) The fact that GRBs occur from the gravitational collapse to a black hole. 

(2) The fact that CBM density for the occurrence of GRBs is inferred from the 
analysis of the afterglow to be on the order of 1 particle/cm 3 (see Tab. 2) 
except for few cases (see e.g. sections 17 and 15). This implies that the process 
of collapse has occurred in a region of space filled with a very little amount 
of baryonic matter. The only significant contribution to the baryonic matter 
component in this process is the one represented by the fircshcll baryon loading, 
which is anyway constrained by the inequality B < 1CP 2 . 

(3) The fact that the energetics of the GRBs associated with supernovas appears 
to be particularly weak is consistent with the energy originating from the grav- 
itational collapse to the smallest possible black hole: the one with mass M just 
over the neutron star critical mass. 

There are also at work very clearly selection effects among the association between 
supernovas and GRBs: 

(1) Many type Ib/c supernovas exist without an associated GRB. 232 

(2) Some GRBs do not show the presence of an associated supernova, although 
they are close enough for the supernova to be observed. 250 

(3) The presence in all observed GRB-supernova systems of an URCA source, a 
peculiar late time X-ray emission. These URCA sources have been identified 
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Fig. 57. A possible process of gravitational collapse to a black hole "induced" by the Ib/c super- 
nova on a companion neutron star in a close binary system. 



and presented for the first time at the Tenth Marcel Grossmann meeting held 
in Rio de Janeiro (Brazil) in the Village of Urea, and named consequently. They 
appear to be one of the most novel issues still to be understood about GRBs. 
We will return on these aspects in the section 17.5. 

The issue of triggering the gravitational collapse instability induced by the GRB 
on the progenitor star of the supernova or, vice versa, by the supernova on the pro- 
genitor star of the GRB needs accurate timing. The occurrence of new nuclear 
physics and/or relativistic phenomena is very likely. The general relativistic insta- 
bility induced on a nearby star by the formation of a black hole needs some very 
basic new developments. 

Only a very preliminary work exists on this subject, by Jim Wilson and his 
collaborators. 280 The reason for the complexity in answering such a question is 
simply stated: unlike the majority of theoretical work on black holes and binary X- 
ray sources, which deals mainly with one-body black hole solutions in the Newtonian 
field of a companion star, we now have to address a many-body problem in general 
relativity. We are starting in these days to reconsider, in this framework, some 
classic work by Ref. 281-287 which may lead to a new understanding of general 
relativistic effects in these many-body systems. This is a welcome effect of GRBs 
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on the conceptual development of general relativity. 

17. Some unexplained features in the GRB-supernova association: 
the URCA sources 

Models of GRBs based on a single source (the "collapsar") generating both the 
supernova and the GRB abounds in the literature. 245 Since the two phenomena 
arc qualitatively very different, in our approach we have emphasized the concept 
of induced gravitational collapse, which occurs strictly in a binary system. The 
supernova originates from a star evolved out of the main sequence and the GRB 
from the collapse to a black hole. The concept of induced collapse implies at least 
two alternative scenarios. In the first, the GRB triggers a supernova explosion in 
the very last phase of the thermonuclear evolution of a companion star. 128 In the 
second, the early phases of the SN induce gravitational collapse of a companion 
neutron star to a black hole. Of course, in absence of a supernova, there is also the 
possibility that the collapse to a black hole, generating the GRB, occurs in a single 
star system or in the final collapse of a binary neutron star system. Still, in such a 
case there is also the possibility that the black hole progenitor is represented by a 
binary system composed by a white dwarf and/or a neutron star and/or a black hole 
in various combinations. What is most remarkable is that, following the "uniqueness 
of the black hole", 288 all these collapses lead to a common GRB independently of 
the nature of their progenitors. 

Having obtained success in the fit of GRB 991216, as well as of GRB 031203 and 
GRB 050315 (see sections 11 and 12), we turn to the application of our theoretical 
analysis to the GRBs associated with supernovas. We start with GRB 980425 / 
SN 1998bw. We must, however, be cautious about the validity of this fit. From the 
available data of BeppoSAX, BATSE, XMM and Chandra, only the data of the 
prompt emission (if < 10 2 s) and of the latest afterglow phases (if > 10 5 s all the 
way to more than 10 8 s!) were available. Our fit refers only to the prompt emission, 
as usually interpreted as the peak of the afterglow. The fit, therefore, represents an 
underestimate of the GRB 980425 total energy and in this sense it is not surprising 
that it does not fit the Ref. 252 relation. The latest afterglow emission, the URCA-1 
emission, presents a different problematic which we will shortly address (see below). 

17.1. GRB 980425 / SN 1998bw / URCA-1 

The best fit of the observational data of GRB 98 4 25 203 ' 206 leads to E^f = 1.2 x 10 48 
erg and B = 7.7 x 10~ 3 . This implies an initial e ± plasma with N e + e - — 3.6 x 10 53 
and with an initial temperature T = 1.2 MeV. After the transparency point, the 
initial Lorentz gamma factor of the accelerated baryons is 70 = 124. The variability 
of the luminosity, due to the inhomogencitics of the CBM, is characterized by a 
density contrast Sn/n ~ 10 _1 on a length scale of A ~ 10 14 cm. We determine 
the effective CBM parameters to be: (n c {, TO ) = 2.5 x 10~ 2 particle/cm 3 and (1Z) = 
1.2 x 10~ 8 . 
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Fig. 58. Theoretical light curves of GRB 980425 prompt emission in the 40-700 keV and 2-26 keV 
energy bands (red line), compared with the observed data respectively from Beppo-SAX GRBM 
and WFC. 203 ' 206 



In Fig. 58 we test our specific theoretical assumptions comparing our theoreti- 
cally computed light curves in the 40-700 and 2-26 keV energy bands with the ob- 
servations by the BeppoSAX GRBM and WFC during the first 60 s of data. 203 ' 206 
The agreement between observations and theoretical predictions in Fig. 58 is very 
satisfactory. 

In Fig. 59 we summarize some of the problematics implicit in the old pie-Swift 
era: data are missing in the crucial time interval between 60 s and 10 5 s, when the 
BeppoSAX NFI starts to point the GRB 980425 location. In this region we have 
assumed, for the effective CBM parameters, constant values inferred by the last 
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Fig. 59. Theoretical light curves of GRB 980425 in the 40-700 keV (red line), 2-26 keV (green 
line), 2—10 keV (blue line) energy bands, represented together with URCA-1 observational data. 
All observations are by BeppoSAX, 203 with the exception of the last two URCA-1 points, which 
is observed by XMM and Chandra. 289 - 290 



observational data. Currently we are relaxing this condition, also in view of the 
interesting paper by Ref. 253. 

The follow-up of GRB 980425 with BcppoSAX NFI 10 hours, one week and 6 
months after the event revealed the presence of an X-ray source consistent with 
SN1998bw, 203 confirmed also by observations by XMM 289 and Chandra. 290 The SI 
X-ray light curve shows a decay much slower than usual X-ray GRB afterglows. 203 
We then address to this peculiar X-ray emission as "URCA-1" (see section 16 and 
the following sections). In Fig. 60A we represent the URCA-1 observations. 203,289,290 
The separation between the light curves of GRB 980425 in the 2-700 keV energy 
band, of SN 1998bw in the optical band, 233,248 and of the above mentioned URCA-1 
observations is evident. 



17.2. GRB 030329 / SN 2003dh / URCA-2 

For GRB 030329 we have obtained 184,185,272 a total energy E^± = 2.12 x 10 52 
erg and a baryon loading B = 4.8 x 10~ 3 . This implies an initial e ± plasma with 
N e + e - = 1.1 x 10 57 and with an initial temperature T — 2.1 MeV. After the 
transparency point, the initial Lorentz gamma factor of the accelerated baryons 
is 70 = 206. The effective CBM parameters are (n c b m ) = 2.0 particle/cm 3 and 
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(TV) = 2.8 x 10 -9 , with a density contrast Sn/n ~ 10 on a length scale of A ~ 10 14 
cm. The resulting fit of the observations, both of the prompt phase and of the 
afterglow have been presented in Ref. 185,272. We compare in Fig. 60B the light 
curves of GRB 030329 in the 2-400 keV energy band, of SN 2003dh in the optical 
band 233 ' 248 and of the possible URCA-2 emission observed by XMM-EPIC in 2-10 
keV energy band. 291 ' 292 

17.3. GRB 031203 / SN 2003lw / URCA-3 




Detector arrival iime (t a ) (s) 



Fig. 60. Theoretically computed light curves of GRB 980425 in the 2-700 keV band (A), of 
GRB 030329 in the 2-400 keV band (B) and of GRB 031203 in the 2-200 keV band (C) are 
represented, together with the URCA observational data and qualitative representative curves for 
their emission, fit with a power law followed by an exponentially decaying part. The luminosity of 
the supernovas in the 3000 — 24000 A are also represented. 233,248 

We will show in section 11 the detailed analysis of GRB 031203 which leads to 
a total energy E l °± — 1.85 x 10 50 erg and to a baryon loading B = 7.4 x 10~ 3 . 
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This implies an initial plasma with N e + e - — 3.0 x 10 55 and with an initial 
temperature T = 1.5 MeV. After the transparency point, the initial Lorentz gamma 
factor of the accelerated baryons is 70 = 132. The effective CBM parameters arc 
{n c bm) = 1-6 x 10 _1 particle/cm 3 and (TV) = 3.7 x 10~ 9 , with a density contrast 
Sn/n ~ 10 on a length scale of A ~ 10 15 cm. In Fig. 60C we compare the light 
curves of GRB 031203 in the 2-200 keV energy band, of SN 20031w in the optical 
band 233 - 248 and of the possible URCA-3 emission observed by XMM-EPIC in the 
0.2-10 keV energy band 199 and by Chandra in the 2-10 keV energy band. 200 

17.4. The GRB / SN / URCA connection 



Table 3. 
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with a power law followed by an exponentially decaying part; d) evaluated assuming a mass of the neutron star M = 1.5M© and 
T ~ 5-7 keV in the source rest frame; e) see Ref. 202,229,273,274. 



In Tab. 3 we summarize the representative parameters of the above three GRB- 
supernova systems together with GRB 060218-SN 2006aj, including the very large 
kinetic energy observed in all supernovas. 293 Some general conclusions on these weak 
GRBs at low rcdshift, associated to SN Ib/c, can be established on the grounds of 
our analysis: 

1) From the detailed fit of their light curves, as well as their accurate spectral 
analysis, it follows that all the above GRB sources originate consistently from the 
formation of a black hole. This result extends to this low-energy GRB class at small 
cosmological redshift the applicability of our model, which now spans over a range 
of energy of six orders of magnitude from 10 48 to 10 54 ergs. 152 ' 180 ' 184 " 186 ' 204 ' 272 Dis- 
tinctive of this class is the very high value of the baryon loading which in one case 
(GRB 060218) 187 is very close to the maximum limit compatible with the dynam- 
ical stability of the adiabatic optically thick acceleration phase of the GRBs. 107 
Correspondingly, the maximum Lorentz gamma factors are systematically smaller 
than the ones of the more energetic GRBs at large cosmological distances. This in 
turn implies the smoothness of the observed light curves in the so-called "prompt 
phase". The only exception to this is the case of GRB 030329. 

2) The accurate fits of the GRBs allow us to also infer some general properties of 
the CBM. While the size of the clumps of the inhomogeneitics is A 10 14 cm, 
the effective CBM average density is consistently smaller than in the case of more 
energetic GRBs: we have in fact (n c b m ) m the range between <~ 10~ 6 particle/cm 3 
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(GRB 060218) and - lO" 1 particle/cm 3 (GRB 031203), while only in the case of 
GRB 030329 it is - 2 particle/cm 3 . 

3) Still within their weakness these four GRB sources present a large variability in 
their total energy: a factor 10 4 between GRB 980425 and GRB 030329. Remarkably, 
the supernova emissions both in their very high kinetic energy and in their bolo- 
metric energy appear to be almost constant respectively 10 52 erg and 10 49 erg. The 
URCAs present also a remarkably steady behavior around a "standard luminosity" 
and a typical temporal evolution. The weakness in the energetics of GRB 980425 
and GRB 031203, and the sizes of their dyadospheres, suggest that they originate 
from the formation of the smallest possible black hole, just over the critical mass of 
the neutron star (see Fig. 57). 

17.5. URCA-1, URCA-2 and URCAS 

We turn to the search for the nature of URCA-1, URCA-2 and URCA-3. These 
systems are not yet understood and may have an important role in the compre- 
hension of the astrophysical scenario of GRB sources. It is important to perform 
additional observations in order to verify if the URCA sources are related to the 
black hole originating the GRB phenomenon or to the supernova. Even a single ob- 
servation of an URCA source with a GRB in absence of a would prove their relation 
with the black hole formation. Such a result is today theoretically unexpected and 
would open new problematics in relativistic astrophysics and in the physics of black 
holes. Alternatively, even a single observation of an URCA source during the early 
expansion phase of a Type Ib/c supernova in absence of a GRB would prove the 
early expansion phases of the supernova remnants. In the case that none of such 
two conditions are fulfilled, then the URCA sources must be related to the GRBs 
occurring in presence of a supernova. In such a case, one of the possibilities would 
be that for the first time we are observing a newly born neutron star out of the 
supernova phenomenon unveiled by the GRB. This last possibility would offer new 
fundamental information about the outcome of the gravitational collapse, and es- 
pecially about the equations of state at supranuclear densities and about a variety 
of fundamental issues of relativistic astrophysics of neutron stars. 

The names of "URCA-1" and "URCA-2" for the peculiar late X-ray emission 
of GRB 980425 and GRB 030329 were given in the occasion of the Tenth Marcel 
Grossmann meeting held in Rio dc Janeiro (Brazil) in the Village of Urea. 294 Their 
identification was made at that time and presented at that meeting. However, there 
are additional reasons for the choice of these names. Another important physical 
phenomenon was indeed introduced in 1941 in the same Village of Urea by George 
Gamow and Mario Schoenberg. 295 The need for a rapid cooling process due to 
neutrino anti-neutrino emission in the process of gravitational collapse leading to 
the formation of a neutron star was there considered for the first time. It was 
Gamow who named this cooling as "Urea process". 296 Since then, a systematic 
analysis of the theory of neutron star cooling was advanced by Ref. 297-301. The 
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coming of age of X-ray observatories such as Einstein (1978-1981), EXOSAT (1983- 
1986), ROSAT (1990-1998), and the contemporary missions of Chandra and XMM- 
Newton since 1999 dramatically presented an observational situation establishing 
very embarrassing and stringent upper limits to the surface temperature of neutron 
stars in well known historical supernova remnants. 302 For some remnants, notably 
SN 1006 and the Tycho supernova, the upper limits to the surface temperatures 
were significantly lower than the temperatures given by standard cooling times. 302 
Much of the theoretical work has been mainly directed, therefore, to find theoretical 
arguments in order to explain such low surface temperature T s ~ 0.5-1.0 x 10 6 K — 
embarrassingly low, when compared to the initial hot (~ 10 11 K) birth of a neutron 
star in a supernova explosion. 302 Some important contributions in this research 
have been presented by Ref. 303-307. The youngest neutron star to be searched for 
thermal emission has been the pulsar PSR J0205+6449 in 3C 58, 307 which is 820 
years old! Ref. 308 reported evidence for the detection of thermal emission from the 
crab nebula pulsar which is, again, 951 years old. 

URCA-1, URCA-2 and URCA-3 may explore a totally different regime: the X- 
ray emission possibly from a recently born neutron star in the first days - months 
of its existence. The thermal emission from the young neutron star surface would 
in principle give information on the equations of state in the core at supranuclear 
densities and on the detailed mechanism of the formation of the neutron star itself 
with the related neutrino emission. It is also possible that the neutron star is initially 
fast rotating and its early emission could be dominated by the magnetospheric 
emission or by accretion processes from the remnant which would overshadow the 
thermal emission. A periodic signal related to the neutron star rotational period 
should in principle be observable in a close enough GRB-supernova system. In order 
to attract attention to this problematic, we have given in Tab. 3 an estimate of the 
corresponding neutron star radius for URCA-1, URCA-2 and URCA-3. It has been 
pointed out 203 the different spectral properties between the GRBs and the URCAs. 
It would be also interesting to compare and contrast the spectra of all URCAs 
in order to evidence any analogy among them. Observations of a powerful URCA 
source on time scales of 0.1-10 seconds would be highly desirable. 

18. Conclusions 

GRBs are giving the first clear evidence for the extraction of energy from black 
holes during the last phases of their formation process. This new form of energy is 
unprecedented in the Universe, both for its magnitude and its very high efficiency in 
transforming matter into radiation, which reaches the 50% limit while the nuclear 
energy reaches efficiency of 2 — 3% only. These sources, with their energy of 10 54 
ergs/pulse, dwarf the corresponding nuclear energy events with their energy of <~ 
10 22 ergs/pulse. 

We have shown how the quest for understanding the initial conditions leading 
to this new gravitational electromagnetic phenomenon has originated a new inquiry 
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into the properties of nuclear matter in bulk which sheds new light on two classical 
physical problems: the heavy nuclei and neutron stars. Such an analysis follows a 
new conceptual unified approach of important heuristic content. It may well lead 
to the solution of yet unsolved problems in physics and astrophysics and on much 
different scales, such as the emission of the remnant in supernova phenomena on a 
macroscopic scale or the limits on the dimension of the constituent of elementary 
particles on a microscopic scale. 

The richness of the experimental data obtained from GRBs, especially thanks to 
the Swift and INTEGRAL satellites, has been exemplified in a selected number of 
sources. This allowed the development and for the first time the testing of the theory 
of ultrarelativistic collisions of baryonic matter with a Lorcntz gamma factor up to 
7 <~ 300 and involving scales extending up to a few light years. The interpretation 
of these observations needed the development of a highly self-consistent theoretical 
framework, which has been tested with high accuracy analyzing both spectra and 
luminosities in selected energy bands on unprecedented time scales ranging from 
few milliseconds all the way to 10 6 seconds. In turn, the CBM structure has been 
analyzed using GRB light curves as its tomographic image. From all these analyses a 
"canonical" GRB scenario is emerging, quite independent on their enormous energy 
ranging from 10 48 erg to 10 54 erg. Such a "canonical" GRB scenario promises to be 
relevant for the future use of GRBs as cosmological standard candles. 

As the comprehension of the GRB phenomenon progresses, so the astrophysical 
setting where they originate is further clarified. We have given evidence for the GRB 
observation originating in active star formation region, with n c b m ~ 1 particle/cm 3 , 
as well as in galactic halos, with n c & m ~ 10~ 3 particlcs/cm 3 . There is also mount- 
ing evidence that all observed GRBs originate from a variety of binary systems. 
There are, in this respect, at least three different systems: 1) GRBs associated with 
supernova originating from an initial binary system formed by a massive star and 
a neutron star and evolving to a neutron star, 2) originating from the supernova 
event, and 3) a black hole, originating from the gravitational collapse of the neu- 
tron star. The occurrence of these GRBs is strictly linked to the process of induced 
gravitational collapse. These supernova-related GRBs are therefore the less ener- 
getic ones, being formed by the smallest possible black holes just over the critical 
mass of a neutron star. There are also much more energetic phenomena occurring 
both in the galactic halos (e.g. GRB 970228) and in star-forming regions (e.g. GRB 
050315), which must definitely originate from the collapse of binary systems formed 
by two neutron stars or a white dwarf and a neutron star or two intermediate mass 
black holes. 

What would appear to be a priori very surprising is that such a large variety 
of initial conditions leads to a "canonical" GRB scenario consistent over a range 
of energies spanning 6 orders of magnitude. Again, the crucial explanation for this 
is due to the uniqueness of the black hole, which was represented in a thought- 
provoking form in the paper "Introducing the black hole" 16 which we reproduce 
here in Fig. 61 and we also represent in its updated version for the GRB connection 
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Fig. 61. The uniqueness of the black hole. Reproduced from the paper "Introducing the black 
hole". 16 

in Fig. 62. 
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